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abstract 


Structural analyses are developed to determine the linear elastic and the geo- 
metncally nonlinear elastic response of an internally pressurised, orthogonally stiff- 
ened, composite material cylindrical shell. The configuration is a long circular cylin- 
drical shell stiffened on the inside by a regular arrangement of identical stringers 
and identical rings. Periodicity permits the analysis of a unit cell model consisting 
of a portion of the shell wall centered over one stringer-ring joint. The stringer- 
ring-shell joint is modeled in an idealised manner; the stiffeners are mathematically 
permitted to pass through one another without contact, but do interact indirectly 
through their mutual contact with the shell at the joint. Discrete beams models of 
the stiffeners include a stringer with a symmetrical cross section and a ring with 
either a symmetrical or an asymmetrical open section. Mathematical formulations 
presented for the linear response include the effect of transverse shear deformations 
and the effect of warping of the ring’s cross section due to torsion. These effects 
are important when the ring has an asymmetrical cross section because the loss of 
symmetry in the problem results in torsion and out-of-plane bending of the ring, 
and a concomitant rotation of the joint at the stiffener intersection about the cir- 
cumferential axis. Data from a composite material crown panel typical of a large 
transport fuselage structure are used for two numerical examples. Although the 
inclusion of geometric nonlinearity reduces the “pillowing” of the shell, it is found 
that bending is localised to a narrow region near the stiffener. Including warping 
deformation of the ring into the analysis changes the sense of the joint rotation. 
Transverse shear deformation models result in increased joint flexibility. 
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CHAPTER 1 
INTRODUCTION 

1.1 COMPOSITE MATERIALS IN PRIMARY STRUCTURES 

Composite materials are being used increasingly for variety of structural applica- 
tions in aerospace engineering and other related weight sensitive applications where high 
strength- to- weight and stiffness-to- weight ratios are required. The success of composite 
materials results from the ability to make use of the outstanding strength, stiffness and 
low specific gravity of fibres such as glass, graphite or Kevlar. When superior specific 
mechanical properties are combined with the unique flexibility in design and the ease of 
fabrication that composites offer, it is no wonder that their growth rate has far surpassed 
that of other materials. 

Development of the state-of-the-art manufacturing techniques has made it possible 
to replace complicated structural components/assemblies by single co-cured or adhesively 
bonded composite parts, thereby minimizing the number of fasteners to be used in a 
structure, and hence, enhancing the structural integrity. While the use of bonded com- 
posite structures as secondary and tertiary load carrying members has been widespread 
in aerospace industry, their use as primary load carrying members is still very limited. 
Most of the applications of composites as primary structural components have been in the 
area of fabrication of empennage or control surfaces of an aircraft. Thus, the potential of 
composite materials as a primary load carrying structure, such as fuselage of an aircraft, 
has not been fully realized yet. One of the main reasons for this could be the lack of 
confidence of aerospace industry in utilizing composite materials for fuselage manufactur- 
ing, which, in turn, could be due to the lack of a full scale analysis, design, and testing 
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to qualify composite materials for use in the fuselage of both civil and military transport 
aircraft. 

1.2 FUSELAGE LOADS AND DESIGN 

As described in the text by Niu 1 , the loads affecting fuselage design of a transport 
aircraft can result from flight maneuvers, landings, cabin pressurization and ground han- 
dling, etc. Fuselage (or cabin) pressurization of a transport aircraft induces hoop and 
longitudinal stresses in the fuselage. The fuselage internal pressure depends on the cruise 
altitude and the comfort desired for the flight crew and/or passengers, and can cause a 
pressure differential of up to 10 psi across the fuselage skin. An unstiffened, or monocoque, 
fuselage would carry this internal pressure load as a shell in membrane response, like a 
pressure vessel. However, internal longitudinal and transverse stiffeners are necessary to 
carry the loads resulting from flight maneuvers, landings, and ground handling, etc. The 
longitudinal stiffeners, called stringers or longerons, carry the major portion of the fuse- 
lage bending moment. The transverse stiffeners, called frames or rings, are spaced at 
regular intervals along the length of the fuselage to prevent buckling of the longitudinals 
and maintain cross-sectional shape of the fuselage. The presence of these internal stiffen- 
ers introduces the following two important aspects in the fuselage design of a transport 

aircraft: 

1. The stiffeners, i.e. stringers and rings, are attached to the fuselage skin by some 
kind of fasteners, or perhaps bonded to it. Thus, there is a transmission of loads between 
the skin and the stiffeners all along their attachment lines, and at the stiffeners’ intersection 
a local concentration of the interacting loads due to joint stiffness occurs. Understanding of 
the load transfer mechanism in the stiffener-to-skin joints under pressurization is necessary 
for determining the load capacity of these joints. 
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2. The presence of internal stiffeners, particularly the presence of frames or rings, 
prevents expansion of the fuselage skin as a membrane, and the skin bulges, or “pillows”, 
between the stiffeners under the action of the internal presure as shown in Fig. 1.1. Hence, 
where the skin is restrained against its expansion as a membrane along the stiffeners, a 
bending boundary layer is formed. 

1.3 STIFFENER-TO-SKIN JOINTS 

The design of stiffener-to-skin joints was cited by Jackson et al. 2 as one of the major 
technology issues in utilizing graphite/epoxy composites in the fuselage of a large transport 
aircraft. In order to realize the full potential of advanced composites in lightweight aircraft 
structure, it is particularly important to ensure that the joints, either adhesively bonded 
or mechanically fastened, do not impose a reduced efficiency on the structure and should 
be cost effective as well. The use of graphite/epoxy composites in conjunction with metal 
fasteners in conventional, mechanically fastened joints is a critical design factor. Improper 
coupling of joint materials can cause serious corrosion problems to metals because of the 
difference in electric potential between these metals and graphite. In other words, insuring 
the galvanic compatibility of fastener materials with graphite composites is essential to 
avoid corrosion problems in the structure 3,4 . 

It has been established that materials such as titanium, corrosion-resistant steels, 
nickel and cobalt alloys can be coupled to graphite composites without such corrosive 
effects. In contrast, aluminum, magnesium and stainless steel are most adversely affected 
because of the difference in electric potential between these materials and graphite, and 
their use would lead to serious corrosion problems in the structure. However, fasteners 
made of materials such as titanium, corrosion-resistant steels, nickel and cobalt alloys 
are much more expensive than the more conventional fastner materials like aluminum, 
magnesium and stainless steel, etc. With thousands of fasteners, e.g., rivets, bolts, nuts. 
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Fig.1.1 Orthogonally stiffened cylindrical shell subjected to internal pressure. 



etc. required to assemble stiffeners to the fuselage skin for a large transport aircraft, 
mechanically fastened joints using corrosion-resistant materials are costly and may offset 
the advantages of using high strength-to-weight composite materials in structures where 
assembly of two or more components is imperative. Niu 1 has pointed out that, in general, 
adhesively bonded joints are more cost efficient for lightly loaded joints, and mechanically 
fastened joints are more cost efficient for highly loaded joints. Thus, if the loads transferred 
between the stiffeners and the fuselage skin are small enough, the adhesively bonded joints 
can be used thereby eliminating all or most fasteners. Elimination of fasteners, or even 
a reduction in the number of fasteners, would enhance the use of advanced composite 
materials in fuselage of a, transport aircraft. As an example, a graphite-epoxy crown panel 
for the fuselage of a large transport aircraft was recently fabricated without fasteners 
by co-curing the stringers and co-bonding the rings, or frames, to the skin 5,6 . Also, the 
curved graphite-epoxy fuselage frames were manufactured by resin transfer molding into 
two-dimensional braided preforms of net structural shape 7 . Clearly, the strength of the 
bond line is a critical issue for these primary fuselage structures made from advanced 
composite materials. 

1.4 CONTACT PROBLEMS 

As described in the monograph by Grigoluk and Tolkachev 8 , contact (or load dif- 
fusion) problems occur in the theory of plates and shells when dealing with interaction 
of plates and shells with rigid and elastic bodies (stamps), stiffening ribs/stiffeners, and 
with plate-shell contacts. The class of contact problems can also include laminated plates 
and shells, if one introduces reactions of interaction between layers. Furthermore, they 
pointed out that the selection of the theory used to formulate the given contact problem 
may also influence the final results. 
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The study of load diffusion in the stiffener-to-skin joints of an orthogonally stiffened 
shell subjected to internal pressure is also a shell contact problem. The type of structural 
theory used to model the discrete elements, i.e„ the shell, the stringer, and the ring, 
influences the distribution of the interacting loads a. the shell-stiffener interface. A brief 
literature survey on the work done in the area of contact problems is presented in the 

remainder of this section. 

Contact problems have always attracted scientists, academicians and designers alike 
because of their inherent importance for any structure analysis involving an assembly of 
two or more components. The first work is by Melan", who considered a semi-infinite 
plate with an infinite stiffener attached to its edge. A concentrated longitudinal force is 
applied to the stiffener. In 1932, Melan obtained a closed form solution for tangential 
forces in the plate along the line of stiffener attachment and also for the axial force in the 
stiffener. Buell'", in 1946, analysed a semi-infinite plate to which a semi-inlinite stiffener 
is attached, loaded at the origin with a longitudinal force. An infinite series solution for 
the airy stress function reduced the problem to an infinite set of algebraic equations, and 
Buell obtained a numerical solution by reducing the set to six equations in six unknowns. 

A solution to Buell’s problem and an identical problem for an infinite plate were obtained 
by Koiter" in 1955. Using as a Green’s function the solution with a concentrated force. 
Roller obtained a singular integral equation for the interacting tangential force between 
the stiffener and the plate. Through a series of complex mathematical steps using Mellin 
transformation. Koiter found the longi.ud.nal force in the stiffener as an infinite series. 
Roller's solution can serve as a criterion of exactness of Buell’s numerical solution. 

The load diffusion problem for a finite stiffener attached to an infinite plate was first 
solved by Benscooter' 1 in 1949. He obtained an Integra-differential equation, oi the same 
form as the Prandtl equation for the distribution of aerodynamic forces in aircraft wing 
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(also known as monoplane equation), with stiffener axial force as an unknown variable. 
First, Benscooter expanded the variable into a series of Chebyshev polynomials of the 
second kind to obtain the discretized equations, and then solved them for unknown coef- 
ficients. Budiansky and Wu 13 extended Melan’s problem for the case where the stiffener 
is rivetted to the plate at discrete points with constant spacing. Subsequent to some of 
these landmark works, numerous authors have studied the load diffusion problem between 
sheet and stiffener. An extensive biblography on the subject is given in Chapter 3 of Ref. 
[ 8 ], 

As for circular cylindrical shells stiffened by longtudinal stiffeners, studies are few. 
Fischer 14 was the first to analyze an infinitely long circular cylindrical shell reinforced by 
equally spaced, continuously attached longitudinal stiffeners, each stiffener being loaded by 
a single concentrated longitudinal force (a counterpart of Melan’s plate problem). Fischer 
accounted for bending of stiffeners and obtained a solution for the membrane shearing 
stress transmitted by a loaded stringer to the shell, and the axial stress developed within 
the stringer. Grigoluk and Tolkachev 8 also analyzed this problem but did not take into 
account the bending of stiffeners. A detailed biblography on some other types of shell 
contact problems can be found in Chapter 8 of Ref. [8]. 

1.5 PRESSURIZED, STIFFENED SHELLS 

A literature survey on the work done in the area of stiffened shells under internal 
pressure suggests that in the past, only a few studies have been carried out in this area. In 
1952, Fliigge 15 studied the stress problems in pressurized cabins of high altitude aircraft 
by dividing it into two problems. First problem was concerned with curved walls of the 
cabin or pressure vessel, hence was called shell problem. The second problem, called the 
plate problem, was concerned with small rectangular panels of the cabin wall, framed 
by stiffeners. Of interest here are the former problems where Fliigge obtained analytical 
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expressions for stresses in the shell and the stiffeners (i.e., stringer and ring) for a single 
cylinder model, and a double cylinder model, using a smeared stiffness approach. In 1958, 
Houghton 16 computed the stresses occuring in stringer reinforced pressurized cylindrical 
shells due to restraining action of the frames. He presented results showing the effect of 
variation of frame pitch and stiffness on the bending moment and shear force in the skins, 
and the hoop stress in the skins between the frames. Houghton’s analysis wa, limited 
metallic components, and did not take into account the eccentricity of stiffeners with 


respect to the skin. Pressure-cabin problems are described in Chapter 9 of Williams" 
1960 text on aircraft structures. In the preface Williams justified the need for a chapter 
devoted to this subject on the importance of high speed civil air-transport. The effect of 
frames and bulkheads on the stresses in a cabin shell was considered in some detail, and 
it was shown how the presence of reinforcing stringers de.locaU.es the constricting effect 
of a frame or bulkhead. Williams analyses were also limited to metallic components, and 


did not take into account the eccentricity of stiffeners with respect to the skin. Wang , m 
1970, carried out a discrete analysis of a metallic, orthogonally stiffened cylindrical shell 
subjected to internal pressure. Stiffener eccentricity, the normal component of interacting 
load between shell and stiffeners, and closed-end pressure vessel effects were taken into 
account. In 1985, Wang and Hsu 1 " improved the earlier work by including in the analysis. 


a composite material shell wall, interacting shear forces between the skin and stiffeners, 
and a direct accounting of closed-end pressure vessel effects. In both of these works, the 
results were obtained for a linear elastic response and symmetric stiffeners. Skin-stiffener 
interactions were computed but results for them were not presented. In 1985, Boitnott- 
examined by experiment and analysis the pressure pillowing of a cylindrical composite 
panel clamped in a stiff fixture. Boitnott’s geometrically nonlinear analysis correlated 
well with the experiments when panel slip from the fixture was taken into account. The 
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analysis showed that the boundary layer decay length decreased with increasing pressure 
and decreasing panel thickness. 

1.6 OBJECTIVES 

In the light of the foregoing discussions, it seems pertinent and timely to extend 
the work in the area of pressurized, stiffened shells by including in the analysis some of 
the features which would improve the understanding of the subject and hence, further 
reinforce the support for the design of a composite material fuselage for a large transport 
aircraft. The objectives of the present research work are to develop analyses of an or- 
thogonally stiffened, laminated composite, cylindrical shell subjected to internal pressure. 

The stiffeners and shell are modeled as distinct elements in order to make available in the 
analyses 

• the distribution of the interacting loads between the shell and stiffeners, and 

• the stress concentration in the shell adjacent to the stiffeners due to “pillowing”. 
Other analysis issues to be addressed in support of these objectives include 

• geometrically nonlinear response versus linear response, and 

• the influence of a ring, or frame, with an asymmetrical open cross section on the 
linear elastic response. 

The intent is to develop analyses that could be used for the design of stiffener- 
to-skin joints and the design of laminated wall construction for the skin. A potential 
benefit of such an analysis/design capability is (i) to use fewer expensive fasteners in the 
graphite/epoxy fuselage, and (ii) to obtain an optimum structural geometery (e.g. shell 
wall thickness and lay-up, frame and stringer stiffnesses, and stiffener spacing etc.) for an 
optimum interface load distribution. 
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1.7 PROBLEM DEFINITION 


An idealized model is assumed for the semi-monocoque fuselage. This configurat.on 
is a dosed-end, stiffened, pressurized shell in which closure is mathematically presumed to 
occur a. infinity. The long circular cylindrical shell is stiffened on the inside by a regular 
arrangement of identical stringers and identical rings (frames). With respect to the applied 
internal pressure load, which is assumed spatially uniform, the model is periodic in the 
circumferential and longitudinal directions both in geometry and in material properties. 
Periodicity of this configuration permits the analysis of a portion of the shell wall centered 
over a generic stringer-ring joint as shown in Fig. 1.2; i.e., deforma., on of a structural 
unit cell (or repeating unit) determines the deformation of the entire stiffened shell. The 
radius of the middle surface of the undeformed cylindrical shell is denoted by R. and the 
thickness of the shell is denoted by t. Axial coordinate a and the circumferential angle S 
are lines of curvature on the middle surface, and the thickness coordinate is denoted by =. 
with -t/2 <z<t/2. The origin of the surface coordinates is centered over the stiffeners 
intersection so that -1 < x < ( and -0 < * < 0, where 21 is the axial length, and 2RB 
is the circumferential me length of the repeating unit. The stringer is assumed to have 
a symmetrical cross section, and the frame is assumed to have either an asymmetrical 
or a symmetrical open section. Asymmetrical open section frames are commonly used 
as transverse stiffeners in the fuselage structure. The stiffeners are modeled as discrete 
beams perfectly bonded to the inside shell wall, so that the interacting loads between 
the stiffeners and shell wall are line load intensities. These line load intensities represent 
resultants of the tractions integrated across the width of the attachment flanges of the 

stiffeners. 

Mathematical formulations for the linear elastic and a geometrically nonlinear elastic 
response are presented in this work. The formulations for the linear elastic response 
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include the effect of transverse shear deformations and the effect of warping deformation 
of the ring’s cross section due to torsion. These effects are important when the ring has 
an asymmetrical cross section, because the loss of symmetry in the problem results in 
torsion of the ring, as well as out-of-plane bending, and a concomitant rotation of the 
joint at the stiffener intersection about the circumferential axis. For symmetric section 
stiffeners, the response of the unit cell (see Fig. 1.2) is symmetric about the stringer axis 
and the ring axis, and there is no rotation of stringer-ring-shell joint. The formulations 
for a geometrically nonlinear response are presented for symmetric stiffeners only, and are 
based on classical theory. The stringer-ring-shell joint is modeled in an idealized manner; 
the stiffeners are mathematically permitted to pass through one another without contact, 
but do interact indirectly through their mutual contact with the shell at the joint. 

On the basis of the symmetry about the x-axis for the unit, only the interacting line 
load components tangent and normal to the stringer are included in the analysis. However, 
due to the ring’s asymmetrical cross section, the components of line loads between shell 
and the ring consist of three force intensities and two moment intensities. The shell- 
stringer interacting force components per unit length along the contact lines are denoted 
by A xs (x) for the component tangent to the stringer, and A„(z) for the component normal 
to the stringer. The three shell-ring interacting force components per unit length along 
the contact lines are denoted by X„(0) for the component acting in the axial direction. 
A , r( *) for the component tangent to the ring, and X„{9) for the component normal to 
the ring. The two shell-ring interacting moment components per unit length along the 
contact lines are denoted by A e r (0) for the component tangent to the ring, and A, r (0) for 
the component normal to the ring. These interacting loads acting in a positive sense on 

the inside surface of the shell are shown in Fig. 1.3. 
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ensities shown in the positive 
of the shell. 
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1.8 ANALYSIS APPROACH 

For both the linear elastic and geometrically nonlinear elastic response of the repeat- 
ing unit to internal pressure, the Ritz method is used. The principle of virtual work is 
applied separately to the shell, stringer and ring. Displacements are individually assumed 
for the shell, stringer, and the ring as Fourier Series expansions. The virtual work func- 
tionals are augmented by Lagrange multipliers to enforce kinematic constraints between 
the structural components of the repeating unit. As a result, point-wise displacement con- 
tinuity between structural elements is achieved. The Lagrange multipliers represent the 
interacting line loads between the stiffeners and the shell, and are also expanded in Founer 
Series. Closed-end pressure vessel effects are included. Data for the example problems are 
representative of the dimensions of large transport fuselage structure. 

The primary advantage of using the analysis approach discussed above results from 
the fact that a point-wise displacement continuity is achieved between the structural 
elements. In commercial finite element analysis codes viz., ABAQUS 21 , NASTRAN' 2 , 
etc., the interpolation functions used for the displacement fields of the structural elements 
(e.g., the shell and beam elements) are, in general, not the same. Thus, the continuity 
between the structural elements can only be satisfied at discrete points, i.e., at the nodes. 
Another shortcoming of these finite element codes is that they can not model the torsional 
warping deformation of an open section, laminated, curved beam. This is a disadvantage 
since the restraint of warping deformation in the ring due to continous contact with 
the shell results in significant circumferential normal stresses in the ring. To account 
for torsional warping deformation in the ring, the ring would have to be modeled as a 
branched shell with these finite element codes. Branched shell models of the stiffeners 
would significantly increase the degrees- of- freedom in the finite element model. It should 
be mentioned that beam models including torsional warping deformation require seven 
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nodal degrees of freedom between elements. The seventh degree of freedom is related to 
the rate of twist. However, it is standard in finite element codes to have only six nodal 
degrees of freedom (three displacements and three rotations) between one-dimensional 
elements. 
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CHAPTER 2 

GOVERNING EQUATIONS FOR LINEAR ANALYSES 


2.1 STRUCTURAL MODEL AND ASSUMPTIONS 

Linear elastic analyses are carried out for a unit cell model (Fig. 1.2) defined in 
Section 1.7 of an internally pressurized, orthogonally stiffened, long circular cylindrical 
shell. The interacting line loads between the shell and stiffeners acting in a positive 
sense on the inside surface of the shell are shown in Fig. 1.3. For a symmetrical section 
ring, the repeating unit (or unit cell model) is symmetric about 0-axis as well, which 
implies that there is no out-of-plane bending and torsion of the ring, and consequently, no 
rotation of the joint at the stiffener intersection about the circumferential axis. Thus, for 
the symmetrical section stiffeners only the interacting line load components tangent and 
normal to the stiffeners are non-zero. 

Mathematical formulations for the linear elastic response presented in this chapter 
include the effect of transverse shear deformations and the effect of warping deformation 
of the ring’s cross section due to torsion. These effects are important when the ring has 
an asymmetrical cross section, because the loss of symmetry in the problem results in 
torsion of the ring, as well as out-of-plane bending, and a concomitant rotation of the 
joint at the stiffener intersection about the circumferential axis. This stringer-ring-shell 
joint is modeled in an idealized manner; the stiffeners are mathematically permitted to 
pass through one another without contact, but do interact indirectly through their mutual 
contact with the shell at the joint. Restraint of cross-sectional warping, as occurs here in 
the ring due to contact with the shell, is an important contributor to the normal stresses in 
thin-walled open section bars, as was demonstrated by Hoff 23 . Based on transverse shear 
deformation and cross-sectional warping of the ring, four structural models are defined. 
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The simplest model uses non-transverse-shear-deformable theory, or classical theory, and 
neglects warping due to torsion. The most complex model includes both effects. Models 
of intermediate complexity occur for inclusion of one effect without the other. 

The purpose of linear elastic analyses is two fold. First, the linear elastic analy- 
sis developed in this chapter is compared with a geometrically nonlinear elastic analysis 
developed in Chapter 3 for the unit cell model with symmetrical cross section stiffeners. 
Second, the effect of skin-stringer-ring joint flexibility, and the effect of warping of the 
ring’s cross section due to torsion, on the response are quantified. The following gen- 
eral assumptions, which are valid for classical as well as transverse shear deformation 
formulations, are made for linear elastic analyses of the unit cell model: 

1. Normals to the undeformed reference surface remain straight and are 
inextensional. 

2. Material behavior is linearly elastic. 

3. The thickness normal stress is assumed to be small with respect to the normal 
stresses in the axial and circumferential directions, and hence it is neglected 
in the material law. 

2.2 TRANSVERSE SHEAR DEFORMATION FORMULATIONS 

2.2.1 SHELL 

A consistent first order transverse shear deformation theory is developed to model 
the shell. Based on the assumption that the shell thickness t is relatively small and hence, 
does not change during loading, the displacements at an arbitrary material point in the 
shell are approximated by 

U(x,6,z) = u(x,8) + z<fi x (x,8) (2.1) 
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V(x,6,z) = v{x,6) + z<f>e(x 7 0) 


( 2 . 2 ) 


W(x,0,z) = w(x,0) ( 2 . 3 ) 

where u(x,0), v(x,0) and w(x,0) are the displacements of the points of the reference 
surface, and <f> x (x,0) and <f>o(x,0) are the rotations of the normal to the reference surface 
as shown in Fig. 2 . 1 . Assuming small displacement gradients, the three-dimensional 
engineering strains are related to the displacements by 



(2.4) 

(2.5) 

( 2 . 6 ) 


in which the polar radius r in cylindrical coordinates is replaced by R + z. Substituting 
Eqs. ( 2 . 1 ) to (2.3) into Eqs. (2.4) to (2.6), and rearranging the terms results in the 


following expressions for the three-dimensional engineering strains: 


^ ^ ^ tee + zk$8 

€xx — €xx T ZK XX €$$ — /I , 7 \ 

(* + *) 

e Z z = 0 

(2.7) 

2 

7x61 T z{\ T -j- 9 ft 

6x6 " (1 + *) 


( 2 . 8 ) 

lOz 

€xz — Ixz €0z — (% 1 z\ 

i 1 + r) 


(2.9) 


in which e xx , n xx , Kee , 7 x $, KxB, KxQ, 7 xz, and 7 g z are the shell strains independent of z- 
coordinate. These shell strains are defined in the following sub-subsection. The transverse 
shear strains e xz and e$ z given in Eqs. (2.9) were obtained through differentiation of Eqs. 
(2.1) to (2.3) with respect to z. However, Eqs. (2.1) to (2.3) are approximate in the 
z-coordinate, so that differentiating with respect to z cannot capture the distribution 
of the transverse shear strains through the thickness of the shell. Since the material is 
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assumed rigid in the 2 -direction (e zz = 0 ), the distribution of the transverse shear strains, 
and consequently the distribution of the transverse shear stresses, does not influence the 
shell behavior. It is the integral of the transverse shear stresses through the thickness, 
or transverse shear resultants, that influences shell behavior. Thus, Eqs. (2.9) should be 
viewed as average values of the transverse shear strains, or as the transverse shear strains 

evaluated at the reference surface (* — 0). 


2.2.1. 1 STRAIN-DISPLACEMENT RELATIONS 

In Eqs. (2.7) to (2.9), the two-dimensional, or shell, strain measures, which are 
independent of the 2 -coordinate, are defined by 


du _ d$x_ (2.10) 



1 dv w 1 

'"^Rdo + R K66 'Rde 

(2.11) 


dv 1 du 

lx6 = !te + Rd6 

(2.12) 


dfo 1 dcj) x J_dv 

hxff ~ dx + R d0 + Rdx 

(2.13) 


dcp e 1 d<j> x l_dv 

Kx6 ~ ~dx R d9 R dx 

(2.14) 


dw , v , ^ 

Ixz = + -fa R + R qq 

(2.15) 

If we set 

the (average) transverse shear strains in Eqs. (2.9) to zero, 

then the 

rotations of the 

normal are 

dw 

<i>x = 

(2.16) 


v 1 dw 

^ = a ~ him 

(2.17) 

so that 

2 d 2 w 2 dv - n 

„ r. „ 4 n x 0 - 0 

h x e -k x0 - R dxde t R q x 

(2.18) 
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Hence, the thickness distribution of the shear strain reduces to 


^x9 


1x9 d~ ^(1 d - 2 f{)^x9 


(i + *) 


(2.19) 


which coincides with the results of Novozhilov’s 24 classical shell theory. 

It is evident from Eq. (2.8) that three shell strain measures are needed to represent 
the distribution of the in-plane shear strain through the thickness in the transverse shear 
deformation shell theory. Whereas, only two shell strain measures are required in classical 
shell theory to represent the shearing strain distribution through the thickness (refer to 
Eq. (2.19) ). Also it can be shown that under rigid body motions of the shell, the nine shell 
strain measures, given by Eqs. (2.10) through (2.15) vanish. (For Novozhilov’s classical 
shell theory, six shell strain measures given by Eqs. (2.10-2.12) and (2.18) vanish under 
rigid body motions.) 


2.2. 1.2 VIRTUAL WORK 


In the three-dimensional elasticity theory, the internal virtual work for the shell is 
given by 




~f I J x \° xx fie xx + veebeee + ^zz^ e zz + ^xe^xO d- ^xz^xz d- cr$z^9z] dV 


( 2 . 20 ) 


where V denotes the volume of shell and dV = (1 + ^) dxR d$ dz . Substitute the variation 
of Eqs. (2.7) to (2.9) into Eq, (2.20), and note that the virtual strains are explicit functions 
of £. Integrals of the stresses with respect to z give force and moment resultants conjugate 
to the shell strains. Hence, the volume integral in Eq. (2.20) reduces to an area integral, 
and the internal virtual work becomes 


«>Vfn 


shell 

t 


-II, 


fit shell ° shell dS, 


( 2 . 21 ) 



where S denotes the area of the reference surface with dS = dxRdd. The generalized 9x1 
stress vector for the shell in Eq. (2.21) is defined by 

5 shell = [IViii N$$ i N$xi M xx , Mg$ , M x $i M x $, Q x -, Q$] > (2.22) 


and the generalized strain vector for the shell is 


€ shell = [ € xxi 1x6, K 66, K x9, MiTizi l6z\ 


(2.23) 


The physical stress resultants and stress couples for the shell, some of which appear in 
Eq. (2.22), are defined in terms of stress components of the symmetric stress tensor in 
cylindrical coordinates by 


(1 + j)dz 


{N xx , M xx ) — J (1»^)^ 

(N»$, M e ») =j{l,z)a $ e dz 
(N xe , M x g) =J(l,z)a x6 {l+^)dz 
( Ng x , M ex ) = J(l,z)a 9x dz 


(2.24) 


Qx — J&xz{l T dz 

Qe =Jcrezdz 

In Eq. (2.22), M x g and M x e are the mathematical quantities conjugate to the modified 
twisting measures k x b and respectively, and are defined in terms of the physical stress 
couples by 

Mxe = \{Mx6 + Mex) M x6 = ^(M x6 - M Sx ) (2.25) 

The nine elements of the stress vector in Eq. (2.22) and the relations of Eq. (2.25) 
determine all the stress resultants and stress couples listed in Eq. (2.24) except for shear 
resultant N x $. The shear stress resultant N xS is determined from moment equilibrium 
about the normal for an element of the shell. This so-called sixth equilibrium equation is 

Mex 


N x $ = N 9x + 


R 


(2.26) 
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Written out in full, the internal virtual work for the shell is given by 

bWint = JJ [N xx he xr + N 9 oSee e + N 6x S'y xe + M xx Sk xx + Meebnee + M x$ 8k x() 
+ M x6 6k X 6 + Qxhxz + Qehez] dS 


The external virtual work for the shell is 

= 6W; he “ + m s x he “ (2.28) 

where 6W* he “ is the external virtual work due to the spatially uniform internal pressure 
load, and 6W s x heli is the external (or augmented) virtual work due to interacting loads. 
The external virtual work for a cylindrical shell under uniform internal pressure, including 
an axial load due to the closed-end effect, is written as 


© 

^ p he “ = J J s P dS + p I d6 [6u(l,0) - 


6u(—l,8)} 


The discussion on the augmented virtual work due to interacting loads is given in Section 


2.2.1.3 CONSTITUTIVE RELATIONS 

The material law for an orthotropic lamina with one material axis in the normal 
direction is given by 

j °xx | Qn Q 12 Q i6 e xx 1 

) r = Qu Q 22 Q 26 < eee > (2.30) 

I. a x6 ) _<5i6 Q 26 Qee _ . e x e J 

where Qij are the transformed reduced stiffnesses given in the text by Jones 25 . (The 
thickness normal stress is assumed to be zero in the material law.) Substitution of Eq. 
(2.30) into Eqs. (2.24) in conjunction with Eqs. (2.25), and use of Eqs. (2.7) to (2.9) for 
the three-dimensional engineering strains, results in the following linear elastic constitutive 
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law for a laminated composite shell wall: 
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in which stiffnesses Aij,Bij and Dij are given by 


(Au,Bn,Dn)=J(l,z,z 2 )Qn( l + ^ dz 

(An, #12, #12) =j(hz,z 2 )Qudz 

{Am, Bn, Dx) = J{ 1 , z,z 2 )Q 22 {^ + dz 

(A w ,Bn)=j(h*)Qi6** 

f - z ~~ ^ 

(A 2 6 ,# 62 ) = j (l,^)Q 26 (l + -g) ^ 

A 66 =^< 966(1 + 

B \ 6 = jQu>z{l + j) dz 

Ble =Jq ^ z ( 1 + 2 ^)^ + ^ dz 

*-j 

^66 =^Q662(i + 2ipn ^ 

D\ e =J Quz 2 {\ + j)dz 

p; 6 =/«. 4 <i2 

#26 =^< 5262 2 (1 + ^X 1 + ^ 


( 2 . 31 ) 


( 2 . 32 ) 
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The lamina, material law relating transverse shear stresses and strains is 



(2.33) 


where 

C44 = + G , 23 5in*'cv 

C45 =(Gi 3 - G , 23 )C , osa 5 ma 
C55 + G'x35 , in 2 a 

in which a is the ply orientation angle. Substitution of Eq. (2.33) into the last two of 
Eqs. (2.24), in conjunction with Eqs. (2.9) for the transverse shear strains, results in 

the following linear elastic constitutive law for a laminated composite shell wall relating 

transverse shear resultants and strains: 



^44 ^45 ( Ixz 1 

^45 ^55 . \ 7 Sz ) 


(2.34) 


The transverse shear stiffnesses, ^44, A45, and ,455 in Eq. (2.34) are given by 


A44 


7 c “ (1 + i 


^45 = / C\$dz 


*55 


■/« 

-l 


)dz 


-1 


(2.35) 


^55(1 + ^) d z 


Since Eqs. (2.9) represent average values of the transverse shear strains, the constitu- 
tive law relating transverse shear resultants and strains, Eq. (2.34), can be viewed as 
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a Hooke’s law based on the assumption of constant transverse shear strain distribution 
through the thickness. Alternatively one can obtain the constitutive law relating trans- 
verse shear resultants and strains based on the assumption of constant transverse shear 
stress distribution through the thickness. A detailed discussion on the subject is given in 
Chapter 2 of the text by Vasiliev 26 . However, both the methods result in a shear correction 
factor of one (as opposed to 5/6) for isotropic materials. Cohen 27 derived the transverse 
shear stiffnesses of laminated anisotropic shells without making either of the assumptions 
mentioned above. He employed Castigliano’s theorem of least work to minimize the shear 
strain energy, and obtained the desired constitutive law, which for homogeneous isotropic 
materials gives a shear correction factor of 5/6. 


2.2. 1.4 EQUILIBRIUM EQUATIONS 

The equilibrium equations for the shell can be derived using the principle of virtual 
work which is stated as 

6W& 11 = sw£f u , ( 2 - 36 ) 


for every kinematically admissible displacement field. For the purpose of deriving the 
equilibrium equations of the shell, the contribution of the augmented virtual work due 
to interacting loads is neglected in Eq. (2.28) for the external virtual work. Thus, sub- 
stituting Eqs. (2.27) and (2.29) for internal and external virtual work, respectively, into 
Eq. (2.36), using the definitions of the strain-displacement relations given by Eqs. (2.10) 
through (2.15), performing integration by parts, and recognizing the arbitrary nature of 
the first variation of the displacements, results in the following set of equilibrium equa- 
tions, or Euler equations, for the shell: 
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9N XX 1 dN 0x 

9x R 99 


6v : 


9N x e 1 gjVga 
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(2.39) 
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(2.40) 

(2.41) 


2.2. 1.5 BOUNDARY CONDITIONS 


The derivation of Euler equations from the principle of virtual work results in the 
following boundary integrals: 
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(N xx —)6u + N x $6v + Q x 6w -f M xx b<\> x + M x q84>$ 
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(2.42) 


+ J N$ x 6u + Neetv + Qe6w + MsxHx + Mee^4>e 

-i 


+0 


dx = 0 


-© 


Boundary integrals in Eq. (2.42) can be made to vanish individually by specifying the 
boundary conditions in two ways. One way is to prescribe periodic boundary conditions; 
alternatively either an essential or a natural boundary condition can be prescribed. In the 
first case, the periodic boundary conditions at x = ±1 are expressed as 

N xe (l,0) = N xe (-l,0), Sv(l,0) = 6v(-l,8) 

QALO ) = QA-U9 ), MM) = 6w(-L0) 

M xx (l,0) = M xx (-l,0 ), M(/,0) = M(-/,0) 

M x0 (l,O) = M x $(—l,0), 6MW) = He(-l,0) 0 € [-0,0] 

For the closed-end pressure vessel effect, prescribe 


(2.43) 


N xx (l,0) = ^ and N xx (—l,0) = ^ 0 € [-0,0] 


(2.44) 
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Thus, 6u(±l,6) is not prescribed to vanish. Periodic conditions at 0 - ±0 edges are 
N er (x,Q) = Nex(x, -0), M*.0) = 


Noe(x,&) = JV^(x,-0), Sv(x,0) = «t>(x,-0) 

Q g(x ,Q) = <? e (x,-0), M*, 0 ) = Mx,-0) 


(2.45) 


Mee(x,&) = M««(x,-0), &<j> x {x,Q) - 6<f> x (x, 0) 

M flx (x,0) = M e ,(x,-0), ^(x,0) = 3 ' e I"M] 

In the second case, the associated boundary conditions at x = ±/ edges are to 

prescribe 

either N xx - f oru but not both, 
either N x e or v but not both, 


either Q x or w but not both, 


either M xx or 4> x but not both, and 
either M x $ or <Ps but not both. 


The associated boundary conditions at 0 = ±0 edges are to prescribe 

either N 9x or u but not both, 
either N$$ or v but not both, 
either Qe or w but not both, 
either Mg x °r 4*x but not both, and 
either Mee or <j>e but not both. 


2.2.2 STRINGER 

Let u,(x) and w 3 (x) denote the axial and normal displacements, respectively, of a 
material point on the stringer reference axis, and let <f> 9s (x) denote the rotation of normal 
as shown in Fig. 2.2. Thus, the axial and normal displacements of a generic material 

point of the stringer are given by 

t/*(x,£) = ti,(x) + ^*,(x) (2,46) 
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W r .(i,0 = «>,(*) 


(2.47) 


respectively. The coordinate system (x,9,Ois located at the centroid of the stringer as per 
the right-hand rule (see Fig. 2.2), in which £ is the normal coordinate. Using Eqs. (2.46) 
and (2.47) and assuming small displacement gradients, the three-dimensional engineering 

strains are 

e xx = Cxs + £, k «s e\ z - 0 e xz ~ T** (2.48) 

which are independent of the 0-direction coordinate because of the symmetric deformation 
assumption. In Eq. (2.48), the one-dimensional strain-displacement relations are defined 

by 

c xs = u' s nes = 4>'os 7 2s = 4>6s + w ' s (2-49) 

in which c rs is the normal strain of the centroidal line, the product £n es is the portion of 
the axial normal strain due to bending, lzs is the transverse shear strain, and the prime 
denotes an ordinary derivative with respect to x. 

The physical force and moment resultants for the stringer in terms of stress compo- 
nents of the symmetric stress tensor are given, in usual way, by 

{N xs , Mgs) =[[ (1-04 dA s 

J J A, (2.50) 

V zs = Jj ^ Kz dA s 

in which N xs is the axial force in the stringer, Mg s is the bending moment, V zs is the trans- 
verse shear force, and is the cross-sectional area of the stringer. Based on transverse 
shear deformation theory, the internal virtual work expression for the stringer is 


sstringer 


6W”t - = t 1^X5 

-/ 


l 

J [N xs 6t xs + Messes + Uzs^Tzs]^ 3 -’ 


(2.51) 


and the Hooke’s law is 

Ajs = {EA)st X s 


A4$s — {EI),kb, 


V^s — {(jA)s'fz 


(2.52) 
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2,2.3 RING 


The structural model is based on transverse shear deformation theory and includes 
cross-sectional warping due to torsion. Warping is a distinctive feature of thin- walled, 
open section beams. Restraint of cross-sectional warping, as occurs here in the ring 
due to its contact with the shell, leads to additional longitudinal normal strains in the 
ring as a result of torsion. The extension of classical thin-walled, open section, curved 
bar theory to laminated composite materials was developed by Woodson, Johnson, and 
Haftka" 8 . However, Woodson et al. did not consider transverse shear deformations. Most 
of the developments for the ring theory presented here are obtained from Woodson’s 
dissertation 29 . The coordinate system (a,0,£) is located at the centroid of the ring as 
per the right-handed system as shown in Fig, 2.3, in which £ is the normal coordinate 
in the radial direction. Let the displacements of a material point on the ring reference 
axis in the x-, 0 and ^-directions be denoted by u r (9 ), u r (0), and w r (8 ), respectively. 
Let the rotations about the x-, 0-, and ("-axes be denoted by <j> xr {0), 4>9 r {8), and <t> zr (0 ), 
respectively. See Fig. 2.3 for the positive sense of these quantities. The displacements of 
a generic point in the cross section are related to the displacements and rotations of the 
point on the reference surface by the approximations 

Ur{x,6 , 0 = u r {8) + C <t>Br{0 ) (2.53) 

hr(x,0,O = V r (9) + (<j> xr (0) + X(f> zr {0) - u{x,()T r {0) (2.54) 

W r (x,8,C) = W r {8) - X(f> 0 r(d) (2.55) 

The cross section of the ring is normal to the 0-axis, that is, the x - ( plane, so that U r 
and W r are interpreted as in-plane displacement components, and V r is the out-of-plane 
component. It is assumed that the cross section is rigid in its own plane. Hence, in-plane 
displacements in Eqs. (2.53) and (2.55) are composed of a translation of the cross-sectional 
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origin plus a small rigid body rotation <f> $r about the 0-axis. The out-of-plane displacement 
given by Eq. (2.54) is composed of a translation of the origin v r , a component (<p xr due 
to bending about the ar-axis, a component x<p zr due to bending about the (“-axis, and a 
component —u>r r due to warping of the cross section out of the flexural plane. In Eq. 
(2.54), w(z,(“) is the warping function for the ring’s cross section, and r v (6) is the twist 
rate which is given after the next equation. The internal virtual work is 
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T M xr 6t\ xr + M zr bn zr + T sr 6r r + M u;r 6(T r / Rq) -f V xr S*( x 


(2.56) 


T ^ zr^*t zt\Ro d0 

in which N$ r is the circumferential force, M xr is the in-plane bending moment, M zr is the 
out-of-plane bending moment, M ur is the bimoment, T sr is the St. Venant’s torque, V xr 
is transverse shear force in the a-direction, V zr is transverse shear force in the (“-direction, 
eor is the circumferential normal strain of the centroidal arc, K xr is the in-plane bending 
rotation gradient, n zr is the out-of-plane bending rotation gradient, 7 rr is the transverse 
shear strain in x-6 plane, *) ZT is the transverse shear strain in 0-( plane, and Rq is the 
radius of ring reference arc. The rotations and strain-displacement relations are 
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in which the over-dot denotes an ordinary derivative with respect to 6. The material law 
is based on the assumption that the shear forces are decoupled from extension, bending, 
and torsional deformations of the ring. Thus, Hooke’s law for the ring is 
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The stiffnesses in Eq. (2.58) are commonly referred to as modulus-weighted section 
properties. The “EH” terms are unique to laminated thin-walled beams (see Bauld and 
Tzeng 30 ). If the laminate construction for each branch of the ring is specially orthotropic 
with respect to X, 8-. and (-directions, then the “EH” terms are all equal to zero. The 
stiffness elements are evaluated from a computer code developed b, Woodson. The reader 
is encouraged to refer to Chapters 2 and 3 of Ref. [291 for further details on this subject. 

The transverse shear stiffness elements in Eq. (2.59) are given by 


GA T e =[A 5 5 h-w + X/^ 66 ) fc 
1 
K 

G A xz —[A4 5 hw + y^( A^)k 
k = i 
K 

GA z e —[^66 hw "1" ^ ^ (^-55 )k 

I t 


(2.60) 


k = 1 

in which the transverse shear stiffnesses, dq, .'is, and dss are calculated based 01 
assumption of constant transverse shear strain distribution through the thickness, and 
are given by Eq. (2.35). In deriving the transverse shear stiffness elements given by Eqs. 
(2.60) above, it is assumed that cross section of the ring is made up of a vertical web and 
horizontal flanges. That is, the web is assumed to be parallel to the q-axis, and flanges 
are assumed to be parallel to the x-axis. In Eqs. (2.60), the parameters h„ and K denote 
the web height and flange width, respectively, and A’ is the total number of flanges in the 

ring cross section. 

For structural models in which the effect of warping of the ring cross section is 
excluded, the contribution of the bimoment, M„ r , to the virtual work of the ring in Eq. 
(2.56) is neglected, and the fourth row and column of the stiffness matrix, Eq. (2.58), are 

ignored. Also, the warping function w(*,C) is taken as zero ’ 
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2.3 CLASSICAL FORMULATIONS 


2.3.1 SHELL 


The shell is modeled with Sanders’ theory 31 , in which first approximation thin shell 
theory is used; i.e., the effects of transverse shear and normal strains are neglected. The 
displacements at an arbitrary material point in the shell are approximated by Eqs. (2.1) 
to (2.3), in which the rotations <p x and <f>$ are related to the displacements by Eqs. (2.16) 
and (2.17), respectively. Thus, transverse shear strains e xz and e$ z in Eqs. (2.9) vanish. 
For small displacement gradients, the three-dimensional engineering strains in Sanders’ 
theory are given by Eqs. (2.7) and 

li+j7L + + z ( l + Tr) kS x6 


€x6 — 


(2.61) 


U + R ) 

in which the quantity k x9 is the twisting strain measure in the Sanders’ theory, which is 
defined in the following sub-subsection. 

2.3.1.1 STRAIN-DISPLACEMENT RELATIONS 

Define a generalized strain vector in terms of the shell strain measures by 

^ shell = [Crxi ^661 ')x 9 i ^xxi ^661 ( 2 . 62 ) 

The first five strain measures of the shell reference surface in Eq. (2.62) are related to the 
displacements by Eqs. (2.10-2.12), and the sixth strain measure, k% 6 , is given by 

(2.63) 


_ dffo 1 d<t> x 1 
16 dx + R d 8 + f?. 


in which the rotation about the normal, <f>~, is given by 

1 dv 1 du 

~ 2 dx~RdO^ (2 ' 64) 

The Donnell- Mushtari- Vlasov (DMV) approximation, or quasi-shallow shell theory is ob- 
tained by neglecting the term ^ in Eq. (2.17) for the rotation <p e , and the rotation about 
the normal 4> z in Eq. (2.63). 
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2.3.1.2 VIRTUAL WORK 


Define a generalized stress vector in terms of the stress resultants and couples of 
Sanders' theory by 

& shell = M xx , Me$, M x q] ( ^ 

such that the interna! virtual work is still given by Eq. (2.21) except that the stress 
and strain vectors are 6 X 1 vectors in Sanders theory. Quantities N‘, and M’, are the 
modified shear and twisting moment resultants. In terms of physical stress and moment 
resultants of the shell these are given by 


Ni 9 = \{N x e + N»x) + }^ M * e ~ Mex) 

M° e = M x6 + Mtx) 


( 2 . 66 ) 

(2.67) 


In the Sanders' original paper 3 ' the term - M„) in Eq. (2.66) was considered 

to be small as compared to J(j V„ + IV.,), and was, therefore, neglected. However, this 
approximation is no. made here. For infinitesimal virtual displacements, the internal 
virtual work for the shell can be obtained by substituting Eqs. (2.62) and (2.65) into Eq. 
(2.21), which results in 

6w mi = IJ [f/ I ,6(„ + N, l Sci + N; l l-l,i + M,rS^, + M»6^i + M^<»} ds < 268 > 

where S denotes the area of the reference surface. The external virtual work expression 
for the classical shell theory is still given by Eq. (2.28). 


2.3.1.3 CONSTITUTIVE RELATIONS 

Consider the material law for an orthotropic lamina given by Eq. (2.30). To get the 
material law for the shell, substitute Eq. (2.30) into the definitions of the resultants in 
terms of stresses, Eqs. (2.24); substitute Eqs. (2.7) and (2.61) for the three-dimensional 
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strains; and then perform the integration with respect to the thickness coordinate. Using 
the definitions of the modified resultants in Eqs. (2.66) and (2.67) gives the final form of 


the material law as 
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(2.69) 


where the stiffnesses An, /li 2 > A 22 , Bn , i?i 2 , B 22 , D\\ , D \2 and D 22 in Eq. (2.69) are given 
by the first three of Eqs. (2.32), with the remaining stiffnesses defined by 

(Ai6,-B 61 ) = J(l,z)Qi e [l + — + 

(A 26 ,B 62 ) =J{1,z)Q 26 [1 + 7^ + 7 WK 1 + 77 ) dz 
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(2.70) 


2.3.2 STRINGER 

The stringer is modeled with Euler-Bernoulli beam theory thereby neglecting the 
transverse shear strain. Hence, equating j zs in the last of Eqs. (2.49) to zero results in 
the following expression for <f>$s : 

4>es = -w' s (2.71) 
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It may be noted that neglecting the transverse shear strain would also modify the virtual 
work statement given by Eq. (2.51), and the third equation in the Hooke’s law, Eq. (2.52), 
is neglected. 

2.3.3 RING 

The ring is modeled with thin-walled, open section, curved bar theory developed by 
Woodson, Johnson, and Haftka 28 . For classical formulations, the transverse shear strains 
are neglected. Hence, equating 7 XT and 7 zr in the last two of Eqs. (2.57) to zero results 
in the following expressions for the rotations <j> xr and (f) zr . 

<t>xr = 2-( Ur - W r ) 4*zr = (2.72) 

It may be noted that neglecting the transverse shear strains would also modify the virtual 
work statement given by Eq. (2.56), and Hooke’s law for the shear resultants, Eq. (2.59), 
is neglected. 

2.4 DISPLACEMENT CONTINUITY 

In order to maintain continuous deformation between the inside surface of the shell 
and stiffeners along their lines of contact, the displacements and rotations should be 
continuous at the shell-stiffener interface. For a symmetrical section stringer, the unit cell 
model is symmetric about x-axis, and the only non-zero displacements for the stringer 
are the axial and normal displacements. The axial and normal displacements at the top 
flange of the stringer in contact with the shell are obtained from Eqs. (2.46) and (2.47) 
for £ = e s , where e s is the radial distance from the stringer centroid to the contact line 
along the inside surface of the shell. Similarly, the corresponding shell displacements 
at the inside surface of the shell (i.e., at z = -tj 2) are obtained from Eqs. (2.1) and 
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(2.3). Hence, the following displacement continuity constraints are imposed along the 
shell-stringer interface (i.e., -I < x < l, 0 = 0). 

g X s = u(x, 0 ) - £<Mz, 0) - [«,(*) + e s <j) es {x)\ = 0 (2.73) 

g za = w{x, 0) - w s (x) = 0 (2.74) 

The asymmetrical section ring bends out-of-plane and twists, in addition to in- plane 
bending and stretching along its circumference. Hence, the displacement field for the ring 
consists of axial, circumferential, and normal components given by Eqs. (2.53) to (2.55). 
From Eqs. (2.54) and (2.55), it can be observed that the circumferential and normal 
displacements of the ring vary along the width of the attachment flange. 

Point-wise continuity of the circumferential displacement between the inside surface 
of the shell and the attachment flange of the ring implies 

V(x,0, — t/2) — V r (x,0, e r ) x € (—6yi, 6^2)1 0 E ( — 0,0) (2.75) 

in which bj\ + bjo = bj > 0 where bj is the width of the attachment flange, and e r is 
the distance from the ring reference arc to the contact line along the inside surface of the 
shell. Since the kinematic assumptions in the ring theory give V r as an explicit linear 
function of x , Eq. (2.54), and the ^-distribution of the shell displacement V, Eq. (2.2), is 
not known apriori, pointwise satisfaction of Eq. (2.75) across the width of the attachment 
flange cannot be achieved. To proceed, the shell displacement is approximated in a Taylor 
series in x about x = 0. That is, 

8V 

V (x,0, -t/2) = V (0,0, —t/2) + x ~q~^ U=o + 0(x~) (2.76) 

Substituting Eq. (2.76) into Eq. (2.75), the continuity of the circumferential displacement 
across the width of the attachment flange can be achieved through terms of order x. Thus, 
V(O,0,—t/2) = V r (O,0,e r ) leads to 

96r = v(o,0)- -4>g( 0,0)- [v r {6) + e r <j> xr (6) - u>oT r (6)] = 0, (2.77) 
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The constraint G 2r = 0 imposed through Eq. (2.78) also implies that the rotation about 
-axis of the shell’s line element tangent to x-curve, f£, equals the rotation of the ring 
around C-axis along their contact line. In Eqs. (2.77) and (2.78), parameters u* and 
are the constant coefficients in the contour warping function, w(*,0 = + *Wi, for the 

attachment flange of the ring. (Thickness warping is neglected and C ^constant along the 
flange contour.) For structural models in which the effect of warping of the ring cross 
section is excluded, the contour warping function w(*,C) is taken 35 zero ‘ 

Similarly, point-wise continuity of the normal displacement between the inside sur- 
face of the shell and the attachment flange of the ring implies 


W(x,0,-t/2) = W r (x,6,e r ) x e{-b fl ,b f2 ), 9 (2.79) 

Since the kinematic assumptions in the ring theory give W r as an explicit linear function 
of x, Eq. (2.55), and the x-distribution of the shell displacement W , Eq. (2.3), is not 
known apriori, pointwise satisfaction of Eq. (2.79) across the width of the attachment 
flange cannot be achieved. To proceed, the shell displacement is approximated in a Taylor 

series in x about x = 0. That is, 

W(x, 0 ,-t/ 2) = W(0,6,-t/2) + x^\ x =o + 0 (x 2 ) (2.80) 

Substituting Eq. (2.80) into Eq. (2.79) the continuity of the normal displacement across 
the width of the attachment flange can be achieved through terms of order x. Thus, 
W(0,0,-t/2 ) = W r ( 0,0, e r ) leads to 

g zr = w(0,0) - w r (0) = 0 ( 2 - 81 ) 
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and x^~\ xr=0 = W r (x, 0, e r ) — W r ( 0,0, e r ) leads to 


Ger= ~fa\x=0 + 4>6r{0) = 0 (2.82) 

The constraint G$ r = 0 imposed through Eq. (2.82) also implies that the rotation about 
0-axis of the shell’s line element tangent to x-curve, is equal to the twist of the ring, 

4>& r , along the contact line. 

Point-wise continuity of the axial displacement between the inside surface of the 
shell and the attachment flange of the ring implies 

(7(x,0,-i/2) = U r (x,0,e r ) x e (- 6 / i , 6 / 2 ), 0 €(-©,©) (2.83) 

Since the kinematic assumptions in the ring theory give U r independent of x, Eq. (2.53), 
and the shell displacement {/, Eq. (2.1), is an arbitrary function of x, pointwise satisfaction 
of Eq. (2.83 ) across the width of the attachment flange can be achieved only through terms 
of order x°. Thus, [7(0,0, -i/2) = U r { 0,0, e r ) leads to 

9xr = u{ 0,0) - ^*(0,0) - [u r {6) + e r <f>$ r (0)\ = 0 (2.84) 

In order to include the axial load sharing between the shell and stringer due to closed- 
end pressure vessel effects directly into the analysis, a seperate constraint is imposed. This 
constraint is that the elongation of the shell at 0 = 0, z = 0 and the elongation of the 
stringer at £ = 0 are the same; i.e. 

[U(l,0,0)-U(-l,0,0)]-[U s (l,0)-U s (-l,0)] =0 (2.85) 

Substituting Eqs. (2.1) and (2.46) into Eq. (2.85) leads to 

[u(/,0)-m(-/,0)] - [u s (l)-u s (-l)} =0 (2.86) 
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2.5 AUGMENTED VIRTUAL WORK FOR THE ASSEMBLY 

Inter-element continuity is enforced by augmenting the virtual work functional with 
the integrals of Lagrange multipliers functions times the variations in the displacement 
constraints. Refering to Fig. 1.3 the Lagrange multipliers are interpreted as the compo- 
nents of the interacting line loads between the stiffeners and the shell, and are defined 
positive if acting on the inside surface of the shell in positive coordinate directions. Thus, 
for the shell, the augmented (or external) virtual work due to the interacting loads is 
l . 

SW^ kel1 = J |a xs (x)[6w(x,0) - -6<f> : r(z,0)] + A zs (x)<5w(x, 0) j dx 
-/ 

© t 
+ J {A xr (0)[MO,0) - ^r(0,*)] + A^)N0,9) - -SM0.9)] (2.87) 

-e 

O w \c 3^ i t 3(f) $ j \ 1 

+ A 2r (0)6to(O,0) - A^ r (0)6( — I^q) + ^zr{9)H^\ x= o “ 2~£h ' 

(R-L)d0-Q[6u(l,O)-6u(-l,O)} 

The axial force Q in Eq. (2.87) is an additional Lagrange multiplier that accounts for axial 
load sharing between the stringer and shell. Similarly, for the stringer, the augmented (or 
external) virtual work due to the interacting loads is 


stringer _ 


A**(®) [£«#(*) + e s 6<f>6s(xj\ + A 2S (x)6tn s (a:)| dx 


+ Q[6u s (l) - 6u s {-l)] 


and. for the ring is given by 


(5>V\ ins = - J (A I . r (0)[<5u r (^) + e r 6(f>er{0)] + A« r (0)[^ w r(0) + e r 6<j) xr (0) u;o6r r (0)] 

+ A, r (0)*tl»r(*) + Aflr {8)64>0t(O) + A, r (0) [W,r(*) - Wl«*v(0)] }(1 + ^)* 


R 0 dd 


(2.89) 
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The displacement constraints (Eqs. (2.73), (2.74), (2.77), (2.78), (2.81), (2.82), 
(2.84) and (2.86)) are enforced by vanishing of the inner product of these equations with 
the variations in the Lagrange multiplier functions. The variational form of these con- 
straints are 

/ 

/ [« xs Qxs T ^^^ 5 ^ 25 ] dx — 0 (2.90) 

-/ 

0 

J \b^xr9xr + ^ffrSBr + ^irffzr + ^A# r GW + 6A zr G zr ] (Rq + e r ) dd = 0 (2.91) 

-© 

«(-/, 0 )] - [«,(/) - «,(-/)]} = 0 (2.92) 
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CHAPTER 3 

GOVERNING EQUATIONS FOR NONLINEAR ANALYSIS 


3.1 ANALYTICAL MODEL AND ASSUMPTIONS 

A geometrically nonlinear elastic analysis is carried out for the orthogonally-stiffened 
cylindrical shell subjected to internal pressure. The structural repeating unit (or unit cell 
model) of Fig. 1.2 is analyzed to obtain response of the entire structure. The shell is 
modeled with Sander’s nonlinear theory of thin shells, and the stiffeners are modeled with 
a -nonlinear Euler- Bernoulli beam theory. The purpose of nonlinear elastic analysis is 
twofold: First, the distributions of interacting loads between the shell wall and the stiff- 
eners are obtained and compared with those obtained from a geometrically linear elastic 
analysis. Second, the influence of geometric nonlinearity on the the stress concentration 
in the shell adjacent to the stiffeners due to “pillowing” is studied. 

Only stiffeners with symmetrical cross sections are considered for the nonlinear re- 
sponse. Hence there is no out-of-plane bending and torsion of the ring, and no rotation of 
the joint at the stiffener intersection about the circumferential axis. Also, on the basis of 
the symmetry about the x- and 6 - axes for the repeating unit, only the interacting line load 
components tangent and normal to the stiffeners are included in the analysis. The shell- 
stringer interacting force components per unit length along the contact lines are denoted 
by A xs (x) for the component tangent to the stringer, and A 2S (x) for the component normal 
to the stringer. The two shell-ring interacting load components per unit length along the 
contact lines are denoted by A flr (0) for the component tangent to the ring and \ zr {B) for 
the component normal to the ring. These interacting loads acting in a positive sense on 
the inside surface of the shell are shown in Fig. 1.3. In this chapter nonlinear formulations 
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are presented for the shell and stiffeners based on classical theory. The following general 


assumptions are made for nonlinear elastic analysis of the repeating unit: 

1. Normals to the undeformed reference surface remain straight and normal 
to the deformed reference surface, and are inextensional. 

2. Material behavior is linearly elastic. 

3. The thickness normal stress is assumed to be small with respect to the normal 
stresses in the axial and circumferential directions, and hence it is neglected 
in the material law. 

3.2 SHELL 


3.2.1 STRAIN-DISPLACEMENT RELATIONS 


Sanders 32 nonlinear theory of thin shells is employed to model the shell. The gen- 
eralized strain vector for the shell is given by Eq. (2.62). Assuming the strains are small 
and rotations are moderately small, the membrane strain-displacement relations are 



in which the rotations <p x , 4>e and (p z are given by Eqs. (2.16), (2.17) and (2.64), respec- 
tively. The change in the normal curvature components, k xx , k^, and k s x6 in terms of the 
shell rotations are linear and are given by Eqs. (2.10), (2.11) and (2.63), respectively. The 
Donnell- Mushtari- Vlasov (DMV) approximation, or quasi-shallow shell theory, is obtained 
by neglecting the rotation about the normal in the strains of Eqs. (3.1) and (2.63), and 
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the term ^ in rotation <f)$ of Eq. (2.17). In this work, it is assumed that the contribution 
of the rotation about the normal to the strains is negligible for the shell. However, the 
DMV approximation is not assumed initially. 


3.2.2 INTERNAL VIRTUAL WORK 


The generalized stress vector in terms of the stress resultants and couples of Sanders' 
theory is defined by Eq. (2.65), and the internal virtual work for the shell in deformed 
state is given by Eq. (2.68). However, the stress resultants and stress couples are now 
defined in terms of second Piola-Kirchhoff stress tensor, which is based on the undeformed 
configuration of the body. Substituting for strains from Eqs. (3.1), Hxx from (2.10), n$e 
from (2.11), and k s x6 from (2.63) into Eq. (2.68), in conjuction with Eqs. (2.16) and (2.17) 
for the definitions of rotations, results in the following expression for the internal virtual 


work for the shell: 
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(3.2) 


3.2.3 EXTERNAL VIRTUAL WORK 

Mathematically the stiffened shell is considered closed at x — ±oo. The work done 
by hydrostatic pressure on an enclosed volume can be derived from a potential energy 
functional, since hydrostatic pressure is a conservative load. However, the infinite volume 
is inconvenient to deal with. Instead, a repeating volume element of the structure that is 
of finite size can be considered, and the potential energy functional for it can be derived. 
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An enclosed volume to contain the pressurized medium can be modeled by joining the 
edges of the repeating unit to the axis of revolution of the cylinder as shown in Fig. 
3.1. The shaded portions are assumed to represent diaphragms (or end caps) enclosing 
the so-formed pie-shaped volume. These diaphragms do not resist deformation of the 
repeating unit, but act to transmit loads normal to the edges of the repeating unit due to 
the internal pressure. Under the action of uniform internal pressure, the right cylindrical 
sector expands to acquire larger volume. Thus, for the deformed volume of Fig. 3.1, the 
external virtual work for the shell due to internal pressure is written as 

m; heil = m; h a eli + m s p h d el1 ( 3 - 3 ) 

where SW shel/ is the virtual work done by the internal pressure acting on the deformed 

ps 

panel area A' B'C'D ' , and iW pd el1 is the virtual work done by the internal pressure acting 
on the deformed diaphragm areas O’A'B' and 0"C'D'. The virtual work done by the 
internal pressure acting on the deformed panels O'A'C'O" and O'B'D'O" is zero due to 
the periodicity of displacement boundary conditions at 8 = i0. 

Expression for 6YV pa elt 

Consider a material point at P on the undeformed reference surface of the shell that 
goes to point P' on the deformed surface. Let the position vector of point P be R and 
of point P' be R* , and further let the displacement vector from point P to P' be U, as 
shown in Fig. 3.2(a). In terms of components along the directions of unit vectors i, t(0), 
and n(0) of a generic point on the undeformed surface, these vectors are represented as 


R = a; \ + R h(0) 

(3.4) 

U = u(x,0 ) i + v(x,0) t{9) + w(x,6) n(0) 

(3.5) 

R* = R -1- U 

(3.6) 
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Using the geometric property of the cross products of position vectors (Fig. 3.2(b)), an 
elemental area dS* of the deformed surface can be written as 

(3.7) 


. c>R* , 5 R* 

dS n* = — - — dx X — — d6 


dx " 86 

where n* is the unit vector normal to the deformed surface element. Substituting Eqs. 
(3.4) and (3.5) into Eq. (3.6), using the resulting expression for R* in Eq. (3.7), and 
noting that || = t and || = -n, the following expression for the elemental area of the 


deformed surface is obtained. 

dS* n* = 
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1 J. i ±( 2 . _ 77 ) 

R 36 1 + R'- 98 + W > R' de v > 

where dS = RdO dx is the elemental area of the undeformed surface. For the hydrostatic 

pressure p acting on the deformed surface A'B'C'D', the external virtual work is given by 


6 W. 


i hell 


= JJ (. pdS * n*) • $U 


K ps —11 \r “ ; * ' “ (3-9) 

where <!>U = 6u(x,6) i + 6v(x,6) i(6) + 6w(x,6) h(6) is the virtual (infinitesimal) displace- 
ment vector, and S * is the area of the deformed reference surface. Upon substituting Eq. 
(3.8) into Eq. (3.9) and carrying out the algebra, the following expression for the first 
component of external virtual work is obtained: 
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(3.10) 

in which S denotes the area of the undeformed reference surface. Equation (3.10) is 
expanded and the terms are rearranged to obtain the following intermediate expression 
for the first component of external virtual work: 
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The terms containing derivatives of the virtual displacements in Eq. (3.11) are integrated 
by parts, and then the terms are further rearranged to obtain the following expression for 
the first component of external virtual work functional: 
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(3.12) 


Expression for SW ^ el1 

The stiffened shell has a periodic symmetry about its axis of revolution in the sense 
that a rotation of the structure through integer multiplies of 20 brings the structure into 
self-coincidence. That is, after a rotation through 20 about the axis of revolution the 
stiffened shell is indistinguishable from its original position. Since the pressure load is 
spatially uniform, the deformed structure exhibits the same periodic symmetry about the 
axis of revolution. That is, a rotation of the deformed stiffened shell through an integer 
multiple of 2© about the axis of revolution brings the deformed shell into self-coincidence. 

The 0-curves of the shell’s reference surface at x = ±/, both in undeformed and 
deformed states, serve to define the end diaphragms. In the undeformed state the 0-curves 
are circles and the end diaphragms are circular areas of radius R. The diaphragms in the 
deformed state are defined by the displacement of diametrical lines. Deformed images 
of diametrical lines remain straight and pass through the axis of revolution. Consider a 
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Axis of revolution 



Q' 


Fig. 3.3 Image of diameter POQ in the deformed state is 
P'O'Q' with O' on the axis of revolution. 
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typical diameter PQ in the undeformed diaphragm as shown in Fig. 3.3. The material 
point at P displaces to P' and the material point Q displaces to Q' . The displacements 
of particles P and Q are such that the vector P r Q' passes through the axis of revolution 
at point O'. The deformed image of every diametrical is a straight line through the same 
point on the axis of revolution. Thus, the point O' on the axis of revolution is common to 
every displaced diametrical line. The displacement of the center of the diaphragm from 0 
to O' is denoted by uoi. The position vectors of points P' and Q' relative to point 0 are 


rp = up i + vp t(0) + (R + wp) n(0) (3.13) 

tq = uq i — vq t(0) — ( R + wq) A(0) (3-14) 

in which unit vectors t(0) and A(0) are tangent and normal to the 0-curve at P. (Note that 
point Q is on the opposite end of the diameter so its tangent and normal are — 1(0) and 
— n(0), respectively.) The displacement of the center of the diaphragm is determined by 
the fact that points P\ O', and Q' lie on a straight line. In vector notation this straight 
line condition is 

(rp - u 0 i) x (rQ - u 0 i) = 0 (3.15) 

Substituting Eqs. (3.13) and (3.14) into (3.15), and carrying out the cross product results 

in two independent conditions, which are 

UQ — U 0 _ R+JWQ 
Up — Uq R + Wp 

VQ _ R + WQ 
vp R + wp 

Solve Eq. (3.16) to get 

uo = fl Up + fo UQ 


(3.16) 

(3.17) 

(3.18) 


where 


fi = 


R + WQ 
2R + wp + wq 


(3.19) 
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(3.20) 


_ R + wp 

2 2 R + wp + wq 

From Eqs. (3.19) and (3.20), it can be seen that the sum /i + h = 1- Equation (3.18) 
shows that the displacement of the center of the diaphragm can be written in terms of the 
displacements of material points on the opposite ends of the diameter of the undeformed 
diaphragm. Since all deformed diametrical lines pass through O', displacement u 0 is 
independent of 0. 

Since displacement uo is independent of 0 , the periodicity of the deformation pattern 
of the 0-curve can be used to write u 0 in terms of the displacements of points on the 0- 
curve in the first unit cell. Consider an odd number of unit cells, such as N r = 3 as shown 
in Fig. 3.4(a), and an even number of unit cells, such as N r = 4 as shown in Fig. 3.4(b). 
For the case of N r = 3. identify point P as C and point Q as F. A clockwise rotation 
of 20 = 2 tt/ 3 brings diameter CF to position EB, and the conditions of self-coincidence 

implies 

u F = u B ,V F = V B ,w F = w B (3.21) 

where point B is at the stiffener (0 = 0) in the first unit. Using Eq. (3.18) the displacement 
of the center of the diaphragm can be written as 


Uq = fl u(0) + fi u (0)» N r odd 


(3.22) 


^ '2R + w{Q) -(- w{ 0) 

R 4- w(0) 

2 ~ 2 R + w(Q) -I- «i(0) 

in which displacements at C correspond to 0 = 0 and displacements at B correspond to 
0 _ o ^ the first unit. For the case of N r = 4, diametrically oppposite points must have 
the same displacement components due to periodicity. That is, a rotation of diameter CG 
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in Fig. 3.4(b) through N r / 2 multiples of 20, i.e., 7r radians, implies uc - v c v g, 
and w c = WG, with 8 = 0 for point C in the first unit. Identify points C and G with 
P and Q, respectively, in Eq. (3.18) such that the periodicity condition gives u 0 = u(Q). 
Similar considerations for diameter BF connecting the two stringers leads to u 0 = u( 0). 
Therefore, for even number of unit cells one gets 

u 0 = u(0) = u(0), N r even (3.23) 

The condition for displacement no for N r odd is more complex than for N r even. 

Consider a point S in the circular sector OAB of the diaphragm at x = l as shown 
in Fig. 3.1. The polar coordinates of S are r and 9 with 0 < r < R and -0 < 9 < 0. 
The particle at S displaces to point 5' in the deformation. The position vector of point 

S' relative to point 0 is 

r *(l,r,d) = r n (0) + u(/,r,0) (3.24) 

in which r n(6>) is the position vector of point S, and u (l,r,0) is the displacement vector 
of S. Since the deformed image of a diameter is a straight line, the displacement vector 
for point S can be interpolated by 

u(/, r, 0) = \ + v(l, 8) i{9) + w(l,9) n(tf)] + (1 ~ ^) u o(0 i (3-25) 

in which u 0 (l) is the axial displacement of the center of the diaphragm at x = /, and 
it is given by Eq. (3.22) if N r is odd and Eq. (3.23) if N r is even. At r = R the 
displacement vector for the diaphragm coincides with displacement vector of the shell s 
reference surface at x = l. The virtual work of the pressure acting on the deformed 
diaphragm O'A'B' referenced to the undeformed diaphragm OAB is 

Wp r'L_ t = / Jr(% dd x lF dr) *= l {3 - 26) 

-0 0 
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Substituting Eqs. (3.24) and (3.25) into Eq. (3.26) results in the following expression for 
the internal virtual work. 


© R 




= If, 

at x~ l J J 


~0 0 


Suo + ji[6u(L0) - Suo] 

x Ou 

, r m 

- u 0 ] 


1 / ck' 


T R bw 


1 + *(§? + ») Titw-”) 


1 4 - - 
1 ^ /? 


rdr cL6 


x=l 


(3.21 


A similar procedure is used on the circular sector of the undeformed diaphragm at 
x = - / . Note that the displacement of the center of the diaphragm at x = -/ is 
and it is given by either Eq. (3.22) or Eq. (3.23) as long as the displacements of the 
0-curve at x = — / are used in these formulas. The external virtual work of the hydrostatic 
pressure on the deformed sector O n C f D f is 

© R 

= -J J p (% dd x % dr ^=-‘ ( 3 - 28 ) 


swtf e “ 


at t~ — l 


-0 0 


where l is replaced by -l in Eqs. (3.24) and (3.25). The total external virtual work of the 
pressure on the diaphragms is the sum of the virtual work of the pressure on each of them. 
Expanding the determinate in Eq. (3.27) and performing integration with respect to polar 
radius r (since r appears explicitly in the integrand), and following a similar procedure for 
the virtual work done by the pressure acting on the end diaphragm 0"C'D' (Eq. (3.28)), 
the second component of the external virtual work is written as 
© 

, .. f { r 1 , - 1 , dv dw . dtn ,1 c 1 c , 

=P J { L* + 2» + ^ (»■ + » 2 ) + ("gj - ( 3 4 ” + 

( dw v , x du 


-© 

1 

+ Jr 
1 

+ Jr 


~ v )( u ~ u °) ~ ~&j( R + w }] 6v 
v % ~ ( u ~ U °^ R + w + , R 


d9 


(3.29) 


The total external virtual work for the shell due to internal pressure can now be 
obtained by substituting Eqs. (3.12) and (3.29) into Eq. (3.3). In this substitution process. 
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the line integrals at edges x = ±1 from Eqs. (3.12) and (3.29) combine and simplify. 
(Several terms combine to add to zero, and integration by parts in 9 are performed to 
eliminate derivatives of the virtual displacements. Periodicity of the displacements at 
0 = ±0 result in vanishing of the boundary terms in the integration by parts.) After 
these manipulations, the following expression for the total external virtual work for shell 
results: 
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(3.30) 


On the basis of periodicity of the shell’s displacements at 9 = ±0 edges, it can be shown 
that the first of the line integrals in Eq. (3.30) vanishes, and the final result for the 
external work functional for the shell under hydrostatic pressure is 
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(3.31) 
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3.3 STRINGER 


For a symmetrical cross section, the internal virtual work statement for the stringer 


nV stringer = j + M$s 6K 6s \dx (3.32) 

in which N xs is the axial force in the stringer, Me s is the bending moment, c xs is the 
extensional normal strain of the centroidal line, and k$ s is the change in curvature of the 
centroidal line. Based on Euler-Bernoulli nonlinear beam theory, the strain-displacement 
relations for stringer are 


€xs = u's + f («;' ) 2 ties = - w " (3.33) 

in which the prime denotes an ordinary derivative with respect to x. Substituting Eq. 
(3.33) into Eq. (3.32) results in the following expression for internal virtual work for the 
stringer. 


X ^stringer 
v r v int 
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OX 
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dx 


Hooke’s law for the stringer is 


(3.34) 


A'xs — (E A) s^XS — (^E) s K^5 


(3.35) 


3.4 RING 

For a ring with symmetrical cross section, there is no out-of-plane bending and 
torsion. Hence, the statement of internal virtual work for the ring is 

e 

6W™ 9 = I [Norteer + M xr Sn xr ]R 0 d8 (3.36) 

-e 
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in which AV is the circumferential force, M xr is the bending moment, e* r is the circum- 
ferential normal strain of the reference arc, k xt is the change in curvature of the reference 
arc, and R 0 is the radius of ring reference arc. Based on small strains and moderate 
rotations, and the kinematic relations for the ring are 


1 1 2 

€&r = + Wr ^ + 2^ x1 

1__ 

Ro 

<Pxr — 77 “ ( W r ) 
ito 


* xr — 


(3.37) 


in which the over-dot denotes an ordinary derivative with respect to 6. Substituting Eqs. 
(3.37) into Eq. (3.36) results in the following expression for internal virtual work for the 

ring: 
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(3.38) 


Hooke’s law for the ring is 


Ner — (EA.)rC$r 


M x r — (El)r^z 


(3.39) 


3.5 INCREMENTAL VIRTUAL WORK 

The internal and external virtual work functionals obtained above contain terms 
that are nonlinear in the displacements times virtual displacements of the first degree. 
Newton’s method is employed to solve the nonlinear equilibrium equations that result 
from the principle of virtual work. It is convenient to derive the incremental virtual 
work functionals for each structural element for implementation of the update procedure 
in Newton’s method. At fixed values of external loads, consider replacing the (actual) 
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displacement field in the virtual work functional of an element bv an approximation. Let 
this approximation be a known displacement field plus a small incremental displacement 
field. The approximation is substituted into the virtual work functional of the element and 
the functional is linearized in the incremental displacement field. As an example of the 
replacement procedure, the case of shell is explained in detail. Let the shell displacement 
field 


{u} — {u} + {Aa} 


(3.40) 


in which {u} on the right hand side of Eq. (3.40) is assumed to be known and (Au) is 
an unknown small increment. As a result of replacement illustrated bv Eq. (3.40), the 
strains and curvatures become 


{c} — - {cl + {Ael 

(3.41 

{k} — - {k} + {An} 

(3.42 


in which {Ac} and {A/c} are linear in {Au}. The stress and moment resultants are 
replaced as 

{N} — {N} + {AN} (3.43) 

{M} — * {M} + {AM} (3.44) 

in which {AN} and {AM} are determined by the constitutive law (Eq. (2.69)) by using 
incremental strains and curvatures; viz., {AT} and {Acc}. 

The virtual displacements are not incremented since they are any of a set of kine- 
matically admissible (test) functions and it is only the (actual) displacements (in the space 
of trial functions) that are being determined by iteration. Consequently, 


b (Au) = 0 


(3.45) 
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Substituting the replacements (Eqs. (3.40) through (3.44)) in the virtual work functional 
(Eq. (3.2) or (3.30)), and linearizing in the increments results in the incremental virtual 
work functionals. The above procedure is repeated for the stiffeners sis well to obtain 
respective incremental virtual work functionals. 


3.5.1 SHELL 


The linearized incremental strains and curvatures for shell are obtained by following 
the incremental procedure described above. Substituting Eq. (3.40) to (3.42) into Eqs. 
(3.1), (2.10), (2.11), and (2.63) in conjunction with Eqs. (2.16) and (2.17), results in the 
following linearized incremental strains and curvatures for shell: 
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Similarly following the replacement procedure explained above, from Eq. (3.2) the lin- 


earized incremental internal virtual work for the shell is obtained as 
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and from Eq. (3.30), the linearized incremental external virtual work for the shell due to 
internal pressure load is obtained as 
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(3.53) 


It may be noted that Eq. (3.53) leads to a symmetric load stiffness since constant pressure 
acting on a fully enclosed volume is a conservative force system. 


3.5.2 STRINGER 


Following the incremental replacement procedure Eq. (3.33) results in the following 
linearized incremental strains and curvatures for the stringer: 
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Ak X s — 


w s 


8 2 A 
dx 2 


(3.55) 


Simularly, from Eq. (3.34) the linearized incremental internal virtual work for the stringer 
is obtained as 
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3.53 RING 

The application of incremental replacement procedure to Eq. (3.37) results in the 
following linearized incremental strains and curvatures for the ring. 
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(3.58) 


Similarly, from Eq. (3.38) the linearized incremental internal virtual work for the ring is 
obtained as 
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3.6 DISPLACEMENT CONTINUITY 


For the symmetrical cross section stiffeners, the repeating unit is symmetric about 
x- and 6- axes, and the only non-zero displacements for the stiffeners are the axial and 
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normal displacements. For a symmetrical cross section, there is no out-of-plane bending 
and torsion of the ring. Thus, axial displacement of the points on the reference arc of the 
ring, U Ty in Eq. (2.53) is zero. In addition the rotations about 6 - and C-axes, 0$ r and 
<f) zry in Eqs. (2.54) and (2.55) for the circumferential and normal displacements of the 
ring, are set to zero. Hence, for the symmetrical cross section stiffeners, the displacement 
continuity constraints imposed between the shell and stringer at their interface are given 
by Eqs. (2.73) and (2.74), and the constraints imposed between the shell and ring at their 
interface are given by Eqs. (2.77) and (2.81). Note that in Eq. (2.77) the twist rate r r is 
equal to zero for a symmetrical cross section ring. 

The advantage of enforcing the displacement continuty constraints, as opposed to the 
strain compatibility constraints (as was done by Wang and Hsu 19 ), between the discrete 
elements is more evident in the nonlinear analysis. The displacement based continuity 
constraints result in equations which are linear in displacement and rotations of the el- 
ements (Eqs. (2.73), (2.77), etc.). On the other hand the strain based compatibility 
constraints would yield nonlinear equations in terms of the displacements and rotations 
of the structural elements (see e.g.. Eqs. (3.1), (3.33) and (3.37)). 

3.7 AUGMENTED VIRTUAL WORK TERMS 

Based on the discussion for the displacement continuity constraints in the foregoing 
section, the augmented (or external) virtual work terms for the shell due to the interacting 

loads are given by 
/ 

wf e,i = J |a ts (x)[^w(x, 0) - l<50 x (x,O)] + A 2S (x)£t«(x,0)) dx 

-i 

0 (3.60) 
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Similarly, the augmented (or external) virtual work terms for the ring due to the inter- 
acting loads are 


© 

s w[ in9 = - J {A fir (0)[aM0) + ^*r(0)] + A, r (0)tftDr(0)} ^ + Y 0 )R ° de (3 ' 61) 
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The augmented (or external) 
still given by Eq. (2.88). The 
are 


virtual work for the stringer due to the interacting loads is 
variational form of the constraints in the nonlinear analysis 
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CHAPTER 4 

FOURIER APPROXIMATIONS AND SOLUTION PROCEDURE 

4.1 INTRODUCTION 

The structural model is periodic in the circumferential and longitudinal directions 
both in geometry and in material properties. With respect to the applied internal pressure 
load, which is assumed spatially uniform, the model is periodic in the circumferential 
direction. For a closed-end pressurized shell, in which closure is mathematically presumed 
to occur at x = ±oo, there is an axial stretching that corresponds to a non-periodic 
axial displacement field. However, for each equal length segment of the shell obtained 
by sectioning it perpendicular to the x-axis, the total axial force is the same. Thus, the 
elongation of each segment is the same. These periodic conditions and uniform axial 
extension due to a closed-end pressure load permit definition of a structural repeating 
unit or unit cell model as shown in Fig. 1.2. 

The periodic nature of the repeating unit or unit cell model requires that the stress 
and moment resultants, and the conjugate displacements and rotations for the discrete 
elements are also periodic in nature so that the repeating units, when placed together, 
form the complete stiffened cylindrical shell model. An analogy to this approach can be 
a jigsaw puzzle with “fully interlocking’' identical pieces. This implies that the boundary 
conditions for the repeating unit are periodic in nature. The periodicity of the forces and 
displacements at the edges of the repeating unit (or in other words, the periodicity of 
boundary conditions) is ensured by assuming a Fourier Series solution to the linear elastic 
and nonlinear elastic response of the repeating unit. The Ritz method is utilized and the 
principle of virtual work is applied separately to each structural element. Displacements 


67 



are seperately assumed for the shell, stringer and ring. The periodic portions of the dis- 
placements and rotations are represented by truncated Fourier Series having fundamental 
periods in the stringer and ring spacing. The non-periodic portions of the displacements 
due to axial stretching are represented by simple terms in x . The Fourier Series reflect 
symmetry about the r-axis only for the repeating unit since the ring can be either sym- 
metrical or asymmetrical in cross section. In this chapter the assumed Fourier Series 
approximations for the actual, virtual and incremental (for nonlinear analysis only) dis- 
placement fields, and the actual and virtual interacting line load intensities are presented. 
The system of discrete equations is obtained and their solution procedure is outlined. 

4.2 DISPLACEMENTS AND ROTATIONS APPROXIMATIONS 


4.2.1 SHELL 

For the shell, actual displacements of the middle surface (see Fig. 2.1) are repre- 
sented as 


M N M N 

u(x,0) = «i mn Sin(a m x)Cos(/3 n 0) + EE u 2mn Cos(a m x)Cos(f3 n 0 ) 

m = l n=0 m=l n=l 


M N 

M N 

(4.1) 

v{x,0) = vi mn Cos(a m 

l x)Sin(/3 n 0) + J2 v ^Sin(a m x)Sin(M) 

(4.2) 

777 = 0 71 = 1 

m= 1 n = 1 


M N 

M N 


w(x,0) = w imn Cos(a m 

x)Cos((3 n 0) + ^ '^2,w 2 mnSin{a m x)Cos(i3 ri 0), 

(4.3) 


m=0 n=0 m = 1 n= 1 


and rotations of the normal are 

M N M N 

*.(*.*)= E E 4>x i mn Sin(a m x)Cos(l3 n 0) + EE 4 > x2mnC os(d Tn x)C o${{3 n 9) (4.4) 

m = 1 n=0 m = l n= 1 

M N M N 

^ ^ )Si7l(fi n 0 ) -j- ^ ^ ^ ^ 4*02mn$ iTl((XmX) S I7l{j3 n 6} (4.5) 

m=0 n=l m = l n=l 
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in which a m = - 1 and (3 n = ~- where m and n are non-negative integers. Coefficient 
qo in the axial displacement field of the shell represents elongation caused by either an 
axial mechanical load or due to closed-end pressure vessel effects. Note that some terms 
in the truncated Fourier Series of Eqs. (4. 1-4.5) have been omitted. The coefficients of 
the omitted terms are u 2 0Oi ^2m0> u 2 on* w 2 mO, <£r 200 , <f>x2mQ , and <f) x2 o n, in which m € Sm 

and n 6 Sn where Sm = {1,2 M} and Sk = {1,2, ...,7V}. The rationale for their 

omission is discussed in Section 4.4. 

The test space of the virtual displacements and rotations is the same space used for 
actual displacements and rotations. The virtual displacements are represented as 

Sq x 

6u(x,0) £ { ~ 2 i~ » { ,u ipoSin(a p x),{ni lpq Sin(a p x)Cos(fiqd),6u2 Pq Cos(a p x)Cos(i3 q B)} 

(4.6) 

6v(x,0) £ {Sv-iOqSin(p q 8),6vi pg Cos(a p x)Sin(P q 0),Sv2p q Sin(Q p x)Sin(/3 q 6)} (4.7) 
6w(x, 8) £ pwioo.^ipoCo^cipa:), 6wio q Cos(P q 8), 6w\ pq C os{a p x)C os(0 q 8), 

(4.8) 

Sw 2 pq Sin{a p x)Cos(p q 0 )} , 
and the virtual rotations of the normal are 

6<j> x (x,8) £ { b4>x\poSin{a p x),64> x \ pq Sin(a p x)Cos{fi q d),b4) x 2 pq Cos(oL p x)Cos(p q 8 )} 

(4.9) 

S(f>e (x,8) £ { 84>8WqSin(j3 q 0),6(p8ipqCos{G p x)Sin(l3 q 0),6<t>82pqSin(a p x)Sin(0 q 8 )} (4.10) 
in which a p = %j- and 0 q = ^ where p £ Sm and q £ Sn. 

4.2.2 STRINGER 

The actual displacements of the centroidal line of stringer (see Fig. 2.2) are 

M M 

^fl X — "V r — r 

= -fi- + u simSin(a m x) + ^ u s2m Cos(a m x) (4.11) 

m — 1 771 = 1 
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(4.12) 


M M 

w s (x) = ^ w s i m Sin(a m x) + w S 2mCos(a m x), 
m=l m= 1 

and the rotation of the normal of the stringer about the 0-axis is 

M M 

4>6s{x) = Z 4 >$s\mSin{a m x) + <f>es 2 mCos(a m x) (4.13) 

m = 1 m=l 

where the coeffcients u s2 o , w s20 and <j> es2 o are omitted. Coefficient q 1 in the axial dis- 
placement field of the stringer represents elongation caused by either an axial mechanical 
load or due to closed-end pressure vessel effects. 

The virtual displacements of the centroidal line of stringer are 

6 u s (x) e {^p-, 6 u s i p Sin{a p x), 6 u S 2 pCos(a p x)} (4.14) 

6w s (x)£ w s ipSin(a p x), 6 w S 2 pCos(a p x )}, (4.15) 

and the virtual rotation of the normal of the stringer about the 0-axis is 

6 d> 6 S (x) G iTi(cy p x ), S(j>$s 2 pCos((XpX )} (4.16) 

where p € S'm . 

4.2.3 RING 

The actual displacements of the reference circle of the ring (see Fig. 2.3) are 

N 

U T (9) = ^UmCoslPnO) (4.17) 

n = l 
A T 

V T (0) = E- Sin((3 n 6) (4.18) 

71 — 1 

N 

Wr (0) = Y, W ™ C °S(M), ( 4 . 19 ) 

n=0 
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and rotations are 


N 

(p rr (0) — ^ " (pxrnS ) 

n=l 

N 

<PeAQ) = T, <P»rnCos(/3 n 0) 

n= 1 

N 

<M*) = £ *Pzrn Sin(f3 n 8) 

?! — 1 

where the coeffcients u r o and <j)« r o are omitted. 

The virtual displacements of the reference circle of the ring are 

6u r (8) £ {6u rg C'os(f3 q d)} 

Sv r {8) £ {6v rq Sin(d q 0)} 
t>iv r (8) £ {6w r0 ,6w rq Cos(f3 q 8)}. 

and the virtual rotations are 

bOj-r ( 9 ) £ {6<j> xrq Sin(f3 q 0)} 

&<P$A8) £ {fi<t>0 r qCos((3 q 6)} 
t0zA8) € {(><i>zrqSin(l3 q 0)} 

where q £ S.\ . 

4.3 INTERACTING LOAD APPROXIMATIONS 

4.3.1 SHELL-STRINGER 

The distributions of the actual interacting line loads bet ween the shell 

(see Fig. 1.3), or Lagrange multipliers, are taken as 

hi M 

Aj* 5 (x) — ^ ' ^xsl??i *5 t'?l( (XjnX ) T ^ ' ^xs2m L- 0.s( Cl m X ) 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 


and stringer 
(4.29) 
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(4.30) 


M M 

A XJ (x) = ^2 Xz S i m Sin(a m x) + ^ A Z s 2 mCos(a m x) 
m = 1 m = l 

where the coefficients A X 2 s 0 and A^so are omitted. 

The test space of the virtual interacting loads is the same space used for actual in- 
teracting load approximations. The shell-stringer virtual interacting loads are represented 
as 

^A x $(x) £ { ^ A x $ j pSi ti( otpX ) * S A x $2pCov»(OpX ) } (4*31) 

<5A. s (x) € {(?A, J! ip5m(Qpx),6A ; : 5 2pC'0'S(apx)}, p £ Sm (4.32) 


4.3.2 SHELL-RING 


The distributions of actual interacting line loads between the shell and ring (see Fig. 
1.3) are taken as 

N 

\ Tr ( 0 ) = ]T A X mCos(l3 n d) (4.33) 

n = l 

A r 

Ag r (6>) = ^2 A 0rnSin{f3 n 6) 

n = 1 

(4.34) 

N 

\. r (9) = ^2 >>zrnCos((3 n 0) 
n = 0 

(4.35) 

N 

A#r(0) = Agrn Cos(j3 n 6) 

n= 1 

(4.36) 

N 

A. r (0) = Azrn Sin((3 n 0) 

n — l 

(4.37) 

where the coefficients A xr o and A^o are omitted. 

The test space of the virtual interacting loads is the same space 
interacting load approximations. The distribution of the shell-ring 

used for the actual 
virtual interacting 


loads are taken as 


6\ xr (8) 6 {<5A xr(J Co.s(/3,j#)} 


(4.38) 
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SXeAO) £ {SX 6rq Sin(fi g 0)} 


(4.39) 


6\~ r (6) £ {SX. r o,6\ zrq Cos((3 q d)} 

(4.40) 

SAer(O) 6 {6A erq Cos(f3 q 0)} 

(4.41) 

6A zr (0) £ {6A ZTq Sin(/3 q 0)}, q £ Sk 

(4.42) 


4.4 TERMS OMITTED IN THE FOURIER SERIES 

Terms omitted in the truncated Fourier Series for the displacements, rotations, and 
interacting loads are determined from the rigid body equilibrium conditions for the ring 
and stringer, and from the displacement continuity conditions between the shell and stiff- 
eners. The external virtual work for the stringer and ring must vanish for any possible 
rigid body motions of these elements. 

4.4.1 RIGID BODY EQUILIBRIUM FOR RING 

Consider the ring in its entirety, that is, as made up of an integer number of repeating 
units around its circumference. Let 0 be the global coordinate such that 0 < 0 < 2ir (see 
Fig. 3.4). The global coordinate 0 is related to the local coordinate 6 in the i th repeating 
unit by 0 = (p { + 0 , i = 1 , . . . , AT where AT is the total number of repeating units and 
4>i = (i — 1 )^r-. Rigid body motion of the ring is considered as a translation of its center 
point 0 to O' plus an infinitesimal rotation about point O'. The displacement of a generic 
point P in the ring is then given by 

U Rigid — U Translation 4” U Rotation (4.43) 

The ring is referred to a fixed cartesian system (x,y,z) with unit basis vectors i, j, and 
k as shown in Fig. 4.1. Also it is convenient to use cylindrical coordinates (.r,0, R 0 + Q 
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with unit vectors i, t (<^>), and ri(d>) as shown in Fig. 4.1. The translational displacement 

of point 0 is given by U Translation = U T i + U y j + f r 2 k. For infinitesimal rotations, the 

rotation can be represented as a vector $1 = Sl x i + Q y j + ft z k. The displacement of P 
relative to 0 due to rotation is 

U Rotation = SI X P (4.44) 

where P = x i + (Ro + C) Sin<f> j + ( Ro + C) Coso k. Thus, the total displacement vector 
in cartesian components becomes 

U Rigid = [U x + (Ro + OttyCoscp - (Ro + £)Sl s Sin<f>] i 

+ \U y - (Ro + O^xCoso + xil z ] j (4.45) 

*F [F z -|- (Ro 4~ £ )Q x Sin(f) — k 

In Chapter 2, the displacement vector of a generic point P is represented in cylindrical 
coordinates as 

U Ring = U r (x,<f>,Qi + V r (Xi<f>iQi(<f>)+ W r (x, <f>, C)n(0) (4.46) 

where the cylindrical coordinate components, f/ r , V r and W r , are given by Eqs. (2.53) to 
(2.55). Note that the local polar angle 6 has been replaced by the global polar angle <p. 
The direction cosines between the unit vectors in cartesian and cylindrical coordinates are 



i 

j 

k 

i 

i 

0 

0 

t 

0 

Cos<t> 

-Simp 

ri 

0 

Sin cp 

Cos0 


Thus, U Rigid in cartesian components, Eq. (4.45), can be transformed to cylindrical 
components, and then like components can be equated between it and Eq. (4.46). As a 
result of this process the displacements of the reference arc of the ring become 

u r —U x + Ro(£lyCos(f) - £t z Sin4>) 

ty =U y Cos<p — U z Sin<f> — Ro&x (4.47) 

uy —UySin(p + U z Cos<J), 


75 


and the (infintesimal) rotations become 

< Pxr = — ^x 

(p 0T =£l y Cos4> — Q z Sin<t> (4.48) 

<p ZT =£l y Sina + Q z Cos<t> 

From the last two of Eqs. (4.48). the rate of twist t t as given by Eq. (2.57) vanishes. It 
may be noted that the strains and changes in the curvatures in Eqs. (2.57) also vanish 
identically for all values of (7xi F y . Fz, fin fly, and Hz* 

Substitution of displacements and rotations from Eqs. (4.47) and (4.48) above into 
the external virtual work for the ring, Eq. (2.89), results in 




[F xr 6 U x + F yr 6U y + F sr 6U t + C x 


m, 


4 “ Cyr&Qy 4 " 


(4.49) 


where 

2 7T 

F xr — J A rr ( Rq + e r ) d<f> (4.50) 

o 

2?r 

F yr = J [A# r Cos(fc + A zr Sin<j>\ (Rq + e r ) d(f> (4.51) 

o 

2 TT 

F zr — J [A rr Cosd) - \$ r Sin<f>] (Ro + c r ) dcp (4.52) 

o 

2 7T 

C xr = J A $r {Ro+e r ) 2 d<f> (4.53) 

0 


2 it 

C yr - I {[*oA„ + \er}Co S <p + \ ZT Sm4>}(Ro + er)d* (4.54) 

o 

2r 

C- r = J {-[Ro\xT + h»r]Sin<j> + \ zr Cos<t>}(Ro + e r )d<f> (4.55) 

o 

Equation (4.49) represents the external virtual work for the complete ring, and is obtained 
by replacing the local polar angle 6 in Eq. (2.89) by the global polar angle <$>. Since the 
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external virtual work for the ring in Eq. (4.49) must vanish for any possible rigid body 

motions, the components of the force and moment resultants defined in Eqs. (4.50-4.55) 

should vanish individually; i.e., F xr = F yr = F zr — C xr = C yr = C ZT = 0. Consider the 

force resultant in redirection for the complete ring as given by Eq. (4.51). Substituting 

(f> = di + 9 and writing the integral as a sum over all units, one gets 

N r 

F yr = Y,[FtCos<Pi + F n Sin<j>i] (4.56) 

i = 1 


in which 



Cosd + \- r Sind (Ro + e r ) dd 



-© 


In the i ih unit the force component Ft is tangent, and component F n is normal, to the 
ring at 0 = <p t . Since the magnitudes of the components F t and F n are the same for each 
repeating unit, Eq. (4.56) can be rewritten as 

;V r A> 

Fyr = F t y CosQi + F n ^2 Sin<t>i (4.57) 

i — 1 t=i 

\ f \ r 

For the complete ring Sincpi - 0. Hence, the total force resultant F yr 

is equal to zero. Similarly, it can be shown that the force component F zr in Eq. (4.52), 
moment component C yr (Eq. (4.54)), moment component C ZT (Eq. (4.55)), all vanish 
for a complete ring. Thus, the rigid body motions that lead to lion-trivial equilibrium 
conditions are the force and moment resultants in the ^-direction. These resultants must 
vanish for each repeating unit. Hence, the non- trivial 2 -direction equilibrium equations 

are e 

J ^rr( d ) (Ro + e r ) dd = 0 (4. 08) 

-0 


and 


0 

J \„ r (9) (Ro + e r ) 2 (16 = 0 


(4.59) 


i i 


-0 



4.4.2 RIGID BODY EQUILIBRIUM FOR STRINGER 


For the stringer the rigid body motions are spatially uniform z-direction and c- 
direction displacements. (A rigid body rotation of the stringer in the x-z plane is not 
considered since this motion would violate longitudinal periodicity of the repeating units.) 
Vanishing of the external virtual work for an arbitrary rigid body displacement of the 

stringer in the axial direction leads to the ^-direction equilibrium equation 

/ 

J A 7S (r) dx = 0 (4.60) 

-/ 

Similarly, the equilibrium equation for a rigid body displacement of the stringer in the 
r-direction is 

J \ zs (x) dx = 0 (4.61) 

-i 

4.4.3 TERMS OMITTED 

Equilibrium Eqs. (4.58), (4.60) and (4.61) imply that coefficients 

A X r0 = 0 A xs oo = 0 A 35 20 — 0 (4.62) 


in the Fourier Series for the interacting loads, and these conditions have been taken into 
account in Eqs. (4.33), (4.29), and (4.30). The sine series for A$ r given in Eq. (4.34) 
satisfies the equilibrium condition given in Eq. (4.59). 

Consider the variational form of the constraints, Eqs. (2.90) and (2.91), for the 
spatially uniform components of the virtual interacting line loads. These equations are 


<*200 


/V 


“<£ x 200 + £(«20n “ “0x2On) — (<*520 + e S <P»s2o) 


sx 


is 20 


= o 


(4.63) 


n~ 1 


N 

[ E wi ° n _ Ws2 ° 

n — 0 




0 


(4.64) 
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6\xrO — ^ 


(4.65) 


^200 - “0x200 + “ -0x2mo) “ i u rG + €r<t>9rO) 

m— 1 

Since these equations are satisfied on the basis that S\ xs 20 = 0, S\ Z s20 = 0 and 

SX xr0 = 0, consistent with Eq. (4.62), the bracketed terms in Eqs. (4.63) to (4.65) 
do not necessarily vanish. The implication that these bracketed terms in Eqs. (4.63) to 
(4.65) do not vanish is that displacement continuity conditions are not satisfied pointwise. 
Pointwise continuity can be achieved by taking each Fourier coefficient appearing in the 
bracketed terms of Eqs. (4.63) to (4.65) to be individually zero. Fourier Series given in 
Eqs. (4.1), (4.4), (4.11), (4.13), (4.17), and (4.21) reflect this choice. Moreover, Fourier 
coefficients t/ooo, u s20> and u r 0 represent rigid body displacement in the axial direction 
for the shell, stringer, and ring, respectively, and setting them to zero can be justified 
on the basis that rigid body displacement does not contribute to the deformation of the 
structural elements. Since Fourier coefficient w s 20 represents rigid body displacement of 
the stringer in the --direction, it would seem that it should be set to zero as well. However, 
to maintain continuity between the stringer and the shell in the --direction, the condition 

N 

ttqon — w 5 20 = 0 (4.66) 

ri=0 


is imposed to determine w s 20 after obtaining the solution for the displacement components 
that deform the shell; i.e., Fourier coefficients Wiq h , n — 1 A 7 , are taken to be non-zero 
independent degrees of freedom since the stringer coefficient w s 20 is not a part of the 
solution vector. 

Finally, consider the constraint equation associated with the spatially uniform 

component of the interacting moment intensity, which is omitted in the series given by 
Eq. (4.36). Derived from Eq. (2.91), this constraint equation is 


[£ 


CXrn^2mO 4 " fydrO 




rO 


= 0 


(4.67) 


79 



The constant component of the twist, o^ r0 , is equated to zero from the considerations 
associated with Eq. (4.65). Consequently, a non-zero value of the uniform component of 
the interacting moment intensity, A« r0 # 0. would not contribute to the equilibrium of the 
ring, since Ag r o and 4>e r o are conjugate variables in the external work for the ring (refer to 
Eq. (2.89)). Since <f>» r o' = 0. it is necessary that A« r0 = 0 to achieve consistent conditions 
for the torsional and out-of-plane bending equilibrium of the ring. With M#ro = 0 in Eq. 
(4.67), the bracketed term does not necessarily vanish, and as a result pointwise rotational 
continuity betwen the shell and the ring is not assured. Pointwise rotational continuity is 
achieved if the coefficients te 2m o = 0, m = 1 M, as is done in the Fourier Series for the 
normal displacement of the shell given by Eq. (4.3). 

4.5 DISCRETE EQUATIONS AND THEIR SOLUTION 

The Ritz method is used to obtain the system of discrete equations. The principle 
of virtual work is applied separately to the shell, stringer, and ring. The virtual work 
functionals are augmented by Lagrange multipliers to enforce kinematic constraints be- 
tween the structural components of the repeating unit. Discrete equations for the linear 
analyses are solved directly. For the nonlinear analysis an iterative solution procedure is 
employed. 

4.5.1 LINEAR ANALYSES 

The approximations in Eqs. (4.1) through (4.28) for the actual and virtual displace- 
ments, and Eqs. (4.29) through (4.42) for the actual and virtual interacting loads are 
substituted into the virtual work functionals for each structural element, and also substi- 
tuted into the variational form of displacement continuity constraints. Then integration 
over the space is performed. This process results in a. 10M N + 13A/ + 14A r + 6 system of 
equations for the transverse shear deformation model and GMN + 10A/ + 11 A' +6 system 
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of equations for the classical model, governing the displacements and the interacting loads. 

These governing equations are of the form 

"An 0 0 B\\ B\2 Bn~ r & shell r An 

0 A 22 0 B 2 \ 0 ^23 sir 0 

0 0 A 33 0 B 22 0 &ring 0 m 

B? B[j 0 0 0 0 | \ 3tr ' 0 

Bn 0 Bf 2 0 0 0 A r , np 0 

. Bjl B? 3 0 0 0 0 J Q J 0 , 

in which submatrices A'n, K 22 and A ’33 are the stiffness matrices for the shell, stringer, 
and ring, respectively, obtained from their respective internal virtual work statements. 
The submatrices B{jJ,j = 1,2,3, in Eq. (4.68) are determined from the external (or 
augmented) virtual work terms involving the interacting loads, and the constraint Eqs. 
(2.73), (2.74), (2.77), (2.78), (2.81), (2.82), (2.84) and (2.86). The vector 011 the right- 
hand-side of Eq. (4.68) is the force vector, determined from the external virtual work 
terms involving pressure. The constraint equations correspond to the last three rows of 
the partitioned matrix in Eq. (4.68). The non-zero elements of the matrices Ahj, B t j and 
An for the transverse shear deformation and classical models are given in Appendices A 
and B, respectively. The symbolic manipulation software Mathematica is used to derive 
these submatrices. The discrete vectors of unknown variables (i.e. shell and stiffeners' 
displacements, and shell-stiffener interacting loads) in Eq. (4.68) for the transverse shear 
deformation model and classical model are given in the following two sub-subsections. 

4.5.1. 1 TRANSVERSE SHEAR DEFORMATION MODEL 

The discrete displacement vector for the shell is the (10MA" + 3A/ + 3N + 2) x 1 
vector 

Ushell = [uo,uJ ,...,Um] T ( 4 . 69 ) 

in which subvectors are 

T 

Uo - [</o, it’ioo, ^ 101 ,^ 101 , $ 0101 , •••,tfioN,wiojv,$ 0 ioiv] (4.70) 
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U m — [ w lmO' ^lmO^xlmOi ^1 ml i w 2ml i V 1 m 1 i v 2ml t ^’lml i u ’2mli 0*1 ml ’ 0x2 ml * 

0^1 ml » 4 > 82m\ ? ••M^lmA r !^2mA T ^lmA' i ^2mN j ^lmA r i W) 2mA'i 0xlmA ; ^ (^* * 

T* 

t Px2mN i &0\ mN i ^()2mA’] 

where m = 1 M 

The (6A/ + 1) X 1 discrete displacement vector for the stringer and (6 N + 1) x 1 vector 
for the ring are 

Ustr = [<7l • Usll.U s 2l,Wsll. U's21 • </>0sll • </>0s21 , u slM • u s2M i 

(4.72) 

W'slA/i tt, s2A/i < ^«4lA/» < / > es2A/J 

T 

U r i n g = [in r o, «rl , »rl , UVl , ^*rl , 4>»rl , 4>zrl , — U r fi,V r N* WrNi <f>xrN, 00 tN,<I>ztn) (4.73) 

in which the term to s 0 for the stringer has been omitted as discussed in reference to Eq. 
(4.66). The AM x 1 discrete interacting loads vector for the shell-stringer interface and 
(5A' -|- 1) x 1 vector for the shell-ring interface are 


A sir — |/^xsl 1 *> A x $2 1 t A^sl 1 i 

A- s 21 i • 

T 

• * i A x 5 1 M ? A ts o A/ i A ^ 5 1 Af ? A - 52 Af ] 

(4.74) 

A ring “ [ A zri) ** A X rl , A$ r ] , A 2 rl ? A#rl 

, A zr i , 

T 

• ♦ • A xrA/ ' A^ r yv ? ^zrN •> ^-OrN * ^zrN ] 

(4.75) 
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where m = 1 M 


The (AM + 1) x 1 discrete displacement vector for the stringer and (4A r + 1) x 1 vector 
for the ring are 


T 

U sir — [</l ? ^sl 1 ' ' W s i i . W s o 1 * • *•- ^slAf * ^s2\l * ^slAI * ^V2 A/] (4./ 9) 


ing — ^rO * ^rl * ^rl ^ ^rl •> 0Br\ * ^rN i VrN i ^-VA’- 0#rA T ] 


( 4 . 80 ) 


The 4 A/ x 1 discrete interacting loads vector for the shell-stringer interface and (5jY + 1 ) x 1 
vector for the shell-ring interface are the same as for the shear deformation model and are 
given by Eqs. (4.74) and (4.75). 

The stiffness submatrices A’n, A’oo and A'33 in Eq. (4.68) are symmetric. From 
Appendices A and B, it may be noted that the elements of the shell stiffness matrix A’n 
are sparsely populated because of the orthogonality of the Fourier Series. The structure 
of A’n is block diagonal; i.e.. nonzero elements occur for index p of the test function equal 
to index m of the trial function, and index q of the test function equal to index n of the 
trial function. See Fig. 4.2. Furthermore, there is a decoupling of the elements obtained 
from the displacement and rotation approximations using single and double Fourier series. 
This special nature of the shell stiffness matrix A'n makes it possible to invert the matrix 
in blocks (i.e., as blocks of 2 x 2 or 6x6 submatrices, etc.) instead of inverting the 
matrix as a whole (i.e., as (6 M A + 2 M + 2 A + 2) x (6M A + 2 M + 2A + 2) matrix). This 
results in considerable saving of the CPU time. As an example consider the linear analysis 
with Fourier Series approximations truncated at 8 harmonics in and ^-directions, i.e.. 
M = A — 8. This results in 418 x 418 A'n matrix for the classical structural model. It 
takes 20 seconds of CPU time on IBM 3090 machine when the shell stiffness matrix A’n 
is inverted in blocks compared to 71 seconds of CPU time on the same computer when 
the matrix is inverted as a whole. 
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or generic values of indices m an 


Equation (4.68) is first solved for the displacements in terms of interacting loads, 
then this solution is substituted into the constraint equations to determine the interacting 
loads. Thus, the total solution is obtained. An LU decomposition procedure is used to 
invert the blocks in the stiffness submatrices, A'n, K 22 and A’ 33 . 

4.5.2 NONLINEAR ANALYSIS 

As stated earlier, the nonlinear analysis is performed for the repeating unit with 
symmetrical section stiffeners only. On the basis of symmetry of the repeating unit about 
x- and 0 -axes, the displacements and interacting loads approximations are modified for the 
nonlinear analysis. In Eqs. (4.1) through (4.3) for shell displacements, in Eqs. (4.11) and 
(4.12) for stringer displacements, in Eqs. (4.17) and (4.20) to (4.22) for ring displacements, 
in Eqs. (4.29) and (4.30) for the shell-stringer interacting loads, and in Eqs. (4.33), (4.36) 
and (4.37) for the shell-ring interacting loads, the coefficients U 9 mn . u> 2 mm u s 2 m* 

Wsi mi Urn* 00m* tyzm* ^xsim * ^xm * A# rn , and A~ rn , all are set to zero where 

777 = 1,2,.... M and n = 1,2, ...,A\ Note that the corresponding coefficients in the virtual 
displacements and interacting load approximations (Eqs. (4. 6-4. 8), (4.14), (4.15), (4.23). 
(4.26-4.28), (4.31). (4.32), (4.38), (4.41), and (4.42)) for the shell and stiffeners are also 
set to zero. These modified approximations for the displacements are substituted into the 
actual and the incremental virtual work functionals derived for each structural element 
in Section 3.5. The incremental displacements have the same functional forms as the 
actual displacements, with the amplitudes denoted by the prefix A. The integration over 
space is performed after substitution for the Lagrange multipliers, or interacting loads, 
and after substitution for each virtual displacement. This process results in a 371 / N + 

6 M + GN + 6 system of equations governing the increment in displacements (indicated 
by a A preceding the displacement subvector symbol) and the Lagrange multipliers. The 
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governing equations are of the form 


' A U shell ' 


Rshell(&shelh P ) 

Aii s ( r 


Rstr{ V-str ) 

Aii r 

► = 

-Rr(u r ) 

x 

* str 

0 

K 


0 

l Q J 


0 


(4.81) 


where the system matrix [C] is given by 



kill Ushell ) pA( ilshell ) 

0 

0 

An 

A 12 

A |3 ‘ 


0 

A oo ( ^str ) 

0 

A 21 

0 

a 23 

[C] = 

0 

0 

A 3 3 ( u r ) 

0 

A 3 2 

0 

Bu 

Bi ; 

0 

0 

0 

0 


Bl o 

o’ 

Bl > 

0 

0 

0 


13 

Bl 3 

0 

0 

0 

0 . 


Submatrices A' n , A'oo and A" 33 are the tangent stiffness matrices for the shell, stringer, 
and ring, respectively, that are functions of the displacements. These matrices are ob- 
tained from their respective incremental internal virtual work statements. The matrix L 
results from the incremental external virtual work functional for the hydrostatic pressure. 
Eq. (3.53), and is a function of displacements as well. The submatrices = 1.2,3. 

in Eq. ( 4.82) are determined from the external virtual work terms involving the Lagrange 
multipliers, and the constraint Eqs. (2.73), (2.74), (2.77) and (2.81). These B,j submatri- 
ces are not functions of the displacements. The vector on the right-hand-side of Eq. (4.81) 
is the residual force vector. The subvectors of R s heii, Rstr , and A, of the residual force 
vector are obtained from the internal and external virtual work (due to internal pressure. 
Eq. (3.30)) statements of the respective structural element. The constraint equations 
correspond to the last three rows of the partitioned matrix in Eq. (4.82). The elements of 
submatrices A\j and L are given in Appendix C, and the elements of subvectors R s helh 
R str and R r are given in Appendix D. The elements of submatrices B,j are obtained from 
those given in Appendix B for the classical model by neglecting the elements corresponding 
to the coefficients set to zero based on the symmetry of the stiffeners. 
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The discrete displacement vector for the shell in Eq. (4.81) is the (3A/A : + 2 M + 


2 N + 2) x 1 vector 


^ shell 


~T - T - T 1 - 

= [ u 0 ' u l “mJ 


(4.83) 


in which the subvectors are 


*0 — [#0* ^00* Voi . U;oi , t’OAN W()A r ] 

(4.84) 

^ m [ ^m0 * ^bnO ? ^ml' V m 1 < ^ m 1 ? • • • » ^ 777 A' ■ Vm N * N ] 

where m = L...,A/. In Eq. (4.81) the ( 2 A/ + 1) x 1 discrete displacement vector for the 
stringer and (2N + 1) x 1 vector for the ring are 

Ustr = [</l , U s j . IV s i , .... U s M,W sM ] T 

T (4 - 85) 

fi r — [uvo i t^ri i KM , .... v r ]\f , uva : ] 

The discrete vectors of the Lagrange multipliers are 

a rp 

A str l • ••« 1 

(4.86) 

Ar = [A rr 0 ? Atfrl ? ^zr\ y •*** ^9rN i A, r ;v] 


The iterative method of solution is based on Eq. (4.81). At a fixed value of the 
pressure p, a sequence of displacements is defined by adding an increment to the previous 
member of the sequence to determine the next member in the sequence. For a. good 
initial displacement estimate, the sequence converges to the displacement solution of the 
nonlinear problem. The initial estimate used here is the converged solution at the last 
pressure load step. The update procedure to determine the increment is based on the 
modified Newton method. In the modified Newton method the matrices /\ u , A ’22 and 
A 33, and L in Eq. (4.82) are only computed for the initial displacement in the sequence, 
and are not updated for each new member in the sequence. As shown in Eqs. (4.81) and 
(4.82) the constraints are applied to the increments in the displacements. If the initial 
estimate of the displacement satisfies these constraint equations and the increments satisfy 
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these same equations, then the final displacement in the sequence will satisfy these same 
equations. 

The convergence of the solution at a fixed load step implies that the equilibiiuin 
of the repeating unit is attained at that load step. Thus, at the converged solution the 
residual defined by S>V ext - S>V int = 6u T R, and the change in displacements A u should be 
sufficiently small. A reasonable criterion for convergence test can be to minimize A u T R. 


where the virtual displacements are replaced by the (admissible) incremental displacement 
vector. An error function defined as A u T R is used to check the convergence against a 
preset tolerance, i.e., A u T R < TOL at the converged solution. In the present analysis, 
this error function is formulated as 


A u 7 R = ABS 
+ ABS 


A U^h e ii[R s hell - BnXgtr — BioX r - BuQ]] 

A U^ tr [—Rstr - Bo^Xstr ~ 


+ ABS 


A U r [ — R r — B 32 Xr 


(4.87) 


The stiffness submatrices (An — ?L). A 22 and A33 in Eq. (4.82) are symmetric. 
Prom Appendix C ' . it may be noted that the elements of the stiffness matiix foi the shell. 
A’n. are densely populated inspite of orthogonality of the series. Furthermore, in the 
nonlinear analysis there is a coupling of the elements obtained from the displacement 
approximations using single and double Fourier Series. Thus, the shell stiffness matrix 
{K u — pL) cannot be inverted in the blocks as is done in the case of linear analysis. 
Instead, the matrix is now required to be inverted as a whole (i.e., as (3MN + 2 Al + 
0 V + 2) x ( 3 M N + ‘2AI + ‘2N + 2) matrix), and hence, is computationally more expensive. 
Equation (4.81 ) is first solved for the displacements in terms of interacting loads, then this 
solution is substituted into the constraint equations to determine the interacting loads. 
Thus, the total solution is obtained. An LU decomposition procedure is used to invert 
the stiffness submatrices. 
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4.5.3 VERIFICATION OF NUMERICAL SOLUTION 


Separate FORTRAN computer programs are written for the linear and nonlinear 
analyses. In each program several checks on the numerical results are coded to establish 
necessary conditions for the solution accuracy. These are 

1. The solution for the Fourier coefficients for the displacements and interacting 
loads are substituted into the left-hand side of the Eq. (4.68) to compute the force vector. 
This computed force vector is compared to the input force vector to check the accuracy 
of the numerical solution to Eq. (4.68). 

2. From free body diagrams of the shell, stringer, and ring, overall equilibrium 
conditions are established for each. These overall equilibrium conditions are evaluated 
using the numerical solution to check if they are satisfied. 

3. Point wise equilibrium equations, or the Euler equations for the functionals, are 
not necessarily satisfied by the Ritz method. However, for the stringer and ring in the 
present analysis, the Ritz solution is an exact solution of the Euler equations as well. 
Consequently, the Euler equations for the stiffeners are programmed and evaluated at 
many points using the numerical solution to assess accuracy. 

4. For the nonlinear analysis the accuracy of the numerical results at each load step 
is ensured by checking the error function defined by Eq. (4.87) against a preset tolerance. 
At the converged solution the left-hand and right-hand sides of Eq. (4.81) are individually 
equal to zero, arid this is verified by the error function given by Eq. (4.87). Furthermore, 
for the converged solution at final load step, the final displacements in the sequence are 
substituted into the displacement constraint equations at the shell-stiffener interface to 
verify them. 

All numerical results presented in the following chapters satisfy these checks on 
numerical accuracy. 
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CHAPTER 5 

A UNIT CELL MODEL WITH SYMMETRIC STIFFENERS 

5.1 INTRODUCTION 

Numerical results are obtained for a linear elastic and a geometrically nonlinear 
elastic response of the unit cell model subjected to internal pressure. The stiffeners are 
assumed to have symmetrical cross sections. The purpose of the linear elastic analysis is 
two fold: First, the results obtained are used to validate the structural model employed in 
the analysis by comparing to the results presented by Wang and Hsu . Second, these ie- 
sults for the linear elastic response are subsequently used to compare with a geometrically 
nonlinear elastic response. Data used in the example are representative of the dimensions 
and cabin pressure of a large transport fuselage structure. 

5.2 NUMERICAL DATA 

Numerical data for the example are R = 117.5 in., 21 — 20 in., 2 RQ = 5.8 in., 
t = 0.075 in., Ro = 113.72 in., e s = 1.10 in., e T = 3.78 in., (EA) r = 0.592 X 10' lb., 
( El ) r = 0.269 x 10 8 lb-in 2 ., (EA), = 0.404 x 10 7 lb., (EI) S = 0.142 x 10 8 lb-in 2 ., and 
with the shell wall stiffness matrices given by 

'0.664 0.221 0 

A = 0.221 0.577 0 x 10 6 Ib/in. B = 0 

0 0 0.221 


' 262 

159 

4.33 

D = 159 

210 

4.33 

4.33 

4.33 

159 


This data was originally used in an example by Wang and Hsu 19 . In their analysis, 
Wang and Hsu neglected the contribution of bend-twist coupling terms Die and D 2 6 in 
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the bending stiffness matrix of shell. This implies that the shell wall is assumed to be 
specially orthotropic. The present analysis is also based on the same assumption. All the 
results presented for this example are for an internal pressure p — 10 psi. The Fourier 
Series are truncated at M = N = 16 (unless otherwise indicated), and this results in 966 
degrees of freedom (or equations to be solved) in the structural model. Since 20 = 2.83°. 
the shell in this example is shallow and the DMY shell theory should be adequate. It is 
found that the numerical results using Sanders theory with the rotation about the normal 
neglected and the numerical results using DMV theory were essentially the same. 

5.3 VALIDATION OF STRUCTURAL MODEL 

The results obtained from the linear elastic analysis are compared with those ob- 
tained by Wang and Hsu 19 to validate the structural model. They presented the results for 
the normal displacement and strains of shell, which are limited to linear analysis. Wang 
and Hsu included the interacting loads in their analysis but did not present results for 
them. The distributions of the shell’s normal displacement w in the circumferential and 
axial directions are shown in Figs. 5.1 and 5.2, respectively, for the linear analysis. The 
^-distributions shown for the linear analysis compare very well with those presented by 
Wang and Hsu (see Fig. 5 in Ref. [19]). The distributions of the circumferential and axial 
normal strains on the inner and outer surfaces of the shell from the linear analysis are 
shown in Figs. 5.3 and 5.4, respectively. The values of the e xx and e## strains compare 
very well to those presented by Wang and Hsu (see Figs. 6 and 7 in Ref. [19]), except in 
one respect. The exception is that the circumferential distribution of the axial strain c XT 
at x = — / (Fig. 5.4) does not exhibit a decrease in value as the stringer is approached. 
Wang and Hsu's results, however, show € xx (Fig. 6 in Ref. [19]) decreasing to nearly 
zero as the stringer is approached along the circumference. Several changes to the shell 
displacement approximations were attempted, but any of these attempts could not give 
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x/l or 0/0 


Fig. 5.3 Circumferential normal strain on the inner and 
outer shell surfaces from the linear analysis at 10 psi. 





Programmed from Wang & Hsu (1985) 
solution 


Present analysis at x = - I 





a solution showing a decrease in the axial normal strain at the stringer. As a last resort 
Wang and Hsu's solution was programmed. However, as shown in Figs. 5.1 to 0 . 4 , the 
solution thus obtained only reconfirmed the results obtained from the current analysis, 
and did not show' a decrease in the axial normal strain at the stringer. 

5.4 LINEAR RESPONSE VERSUS NONLINEAR RESPONSE 

5.4.1 PILLOWING 

Circumferential distributions of the normal displacement w for the shell are shown 
in Fig. 5.1 for the linear analysis and in Fig. 5.5 for the nonlinear analysis. Axial 
distributions of the normal displacement for the shell are shown in Fig. o.2 for the linear 
analysis and in Fig. 5.6 for the nonlinear analysis. For reference, the normal displacement 
for the unstiffened shell, or membrane reponse, is w = 0.2287 inches for the linear analysis, 
and w = 0.2290 inches for the nonlinear analysis. The presence of the stiffeners reduces 
the normal displacements from these membrane values as is shown in these figures. The 
pillowing effect is much more pronounced for the linear analysis (Figs. 5.1 and 5.2) than 
for the nonlinear analysis (Figs. 5.5 and 5.6). The largest normal displacement occurs 
midway between the stiffeners, and this value for the linear analysis is 0.1796 inches 
while it is 0.1541 inches when geometric nonlinearity is included. The minimum normal 
displacement occurs at the stiffener intersection, and its value is 0.1392 inches in the 
linear analysis and 0.1490 inches in the nonlinear analysis. Normal displacements along 
the stiffeners vary only slightly from their values at the intersection for both analyses. 
Thus, including geometric nonlinearity in the analysis increases the minimum value of the 
normal displacement of the shell and decreases its maximum value, which is an indication 
that pillowing is reduced in the nonlinear response. 
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The circumferential and axial normal strain distributions on the inner and outer 
surfaces of the shell also show the reduced pillowing effect in the nonlinear analysis. See 
Figs. 5.7 and 5.8. The circumferential bending strain (difference in e## between the inner 
and outer surfaces) is maximum at the stringer midway between the rings (Figs. 5.3 
and 5.7), and axial bending strain is maximum at the ring midway between the stringers 
(Figs. 5.4 and 5.8). These maximum bending strains are substantially reduced in the 
geometrically nonlinear response. 

5.4.2 BENDING BOUNDARY LAYER 

Compare the circumferential normal strain distributions midway between the ring 
stiffeners ( x = ±/) from the linear analysis (Fig. 5.3) to the nonlinear analysis (Fig. 5.7). 
The bending strain magnitudes, which are differences between the outer and inner normal 
strain values, are less in the nonlinear response than in the linear response. In the linear 
analysis the shell exhibits bending over its entire circumference. In nonlinear analysis the 
shell behaves like a membrane in the central portion with bending confined to a narrow 
zone, or boundary layer, adjacent to the stringer. Since the bending is confined to a narrow 
zone near the stringer, the strain gradients are larger in the nonlinear response than in the 
linear response. These observations are confirmed by plotting the circumferential bending 
moment Moo and the circumferential transverse shear resultant Q$ versus 0/Q at x = — / 
as is done in Figs. 5.9 and 5.10, respectively. The bending moment magnitude is less in the 
nonlinear analysis compared to the linear analysis. Moreover, for the nonlinear analysis 
the bending moment is significantly different from zero only in the boundary layer. This 
is similar to the distributions of the circumferential strains. By the shell equilibrium, 
the transverse shear resultant Qe is proportional to the derivative (gradient) of bending 
moment, dM$$fd6. Because the bending moment gradients, or the strain gradients, in 
the nonlinear analysis are larger than in the linear analysis, the transverse shear resultant 
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- 0.2 


i on the inner and 
r analysis at 10 psi. 












has a larger magnitude in the nonlinear analysis than in the linear analysis as is shown 
in Fig. 5.10. A larger transverse shear resultant in the nonlinear analysis means that the 
interlaminar shear stresses are larger in the nonlinear analysis than in the linear analysis. 
Thus, while pillowing is reduced in the nonlinear response, the confinement of bending 
to a boundary layer near the stringer results in larger interlaminar shear stresses near 
the stringer in the nonlinear response than in the linear response. These features of the 
nonlinear response are consistent with the results found by Boitnott 20 . 

5.4.3 INTERACTING LOAD DISTRIBUTIONS 

The distributions of the interacting line loads between the ring and the shell are 
shown in Figs. 5.11 and 5.12. The distributions of the circumferential component, A Sr , are 
antisymmetric about the origin, and A# r has reduced magnitudes due to the geometrically 
nonlinear effect. As shown in Fig. 5.12, the distributions of the normal component of the 
interacting load, A zr , are symmetric about the origin, attain extremum at the origin, and 
exhibit severe gradients at the origin. The negative value of A zr at the origin indicates 
that the action of the ring is to pull the shell radially inward against the action of the 
pressure to expand the shell outward. The peak normal load intensity is changed from 
-1,674 lb/in. in the linear analysis to -1,045 lb/in. in the nonlinear analysis. 

The distributions of the interacting line loads between the stringer and the shell are 
shown in Figs. 5.13 and 5.14. The distributions of the tangential component, X xs , are 
antisymmetric about the origin and A IS has reduced magnitudes due to the geometrically 
nonlinear effect as shown in Fig. 5.13. The distributions of the normal component, \ zs , 
are symmetric about the origin as shown in Fig. 5.14. The normal component A. s is 
maximum at the origin and has a. steep gradient there. The maximum value of normal 
component is reduced from 484,7 lb /in . in the linear analysis to 320.3 lb/in. in the 
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Fig. 5.14 Normal component of the line load acting 
on the shell due to the stringer at 10 psi. 
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Fig. 5.15 Stringer axial force and bending moment 
distributions at 10 psi. 
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magnitude of the interacting normal load intensity at the stiffener intersection as M and 
N increase. The total normal load intensity acting on the shell at the stiffener intersection 
is denoted by A z , where X z — A rr (0) + A r5 (0). This value is -1,190 lb/in. from the linear 
analysis and is changed to -725 lb/in. in the nonlinear analysis for M = N = 16. A 
normalized value of X z is plotted versus increasing M and N values, with M = N , in Fig. 
5.17 for the linear analysis, and in Fig. 5.18 for the nonlinear analysis. The normalization 
factor, X zmax , is simply the value of X z for the largest values of M and N considered in 
each analysis. As shown in figures, X z is steadily increasing with an increasing number of 
terms in the truncated Fourier Series. Consequently, the series for X z does not exhibit, 
in the range of M and N considered, a convergent behavior. In spite of the fact that the 
normal load intensity at the stiffener intersection is exhibiting singular behavior, the total 
radial resultant load at the stiffener intersection converges rapidly with increasing values 
of M and N as shown in Figs. 5.17 and 5.18. The total radial resultant plotted in these 
figures is defined by 

0 

F z = J [X zr Cos6 - X Q r SinO ] (R 0 + e T ) d6 , (5.1) 

-© 

since the resultant from the stringer vanishes by Eq. (4.61). (A more general approach to 
the resultants at the stiffener intersection is discussed in subsection 6.3.2.) From the linear 
analysis F z — —357.5 lbs, and from the nonlinear analysis F z — —382.8 lbs. Since the 
applied radial load due to internal pressure acting on the repeating unit is 1,160 lbs (= 10 
psix20 in.x5.8 in.), the ring resists about 30.8% of this applied load in the linear response, 
and this percentage is increased to 33% in the nonlinear response. The remaining portion 
of the applied radial pressure load is carried by the shell. 
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magnitude of the interacting normal load intensity at the stiffener intersection as M and 
N increase. The total normal load intensity acting on the shell at the stiffener intersection 
is denoted by A,, where A 2 = A, r (0) + A- s (0). This value is -1,190 lb/in. from the linear 
analysis and is changed to -725 lb/in. in the nonlinear analysis for M = N = 16 . A 
normalized value of A, is plotted versus increasing M and N values, with M = N, in Fig. 
5.17 for the linear analysis, and in Fig. 5.18 for the nonlinear analysis. The normalization 
factor, \, ma x, is simply the value of A 2 for the largest values of M and N considered in 
each analysis. As shown in figures, A 2 is steadily increasing with an increasing number of 
terms in the truncated Fourier Series. Consequently, the series for A 2 does not exhibit, 
in the range of M and N considered, a convergent behavior. In spite of the fact that the 
normal load intensity at the stiffener intersection is exhibiting singular behavior, the total 
radial resultant load at the stiffener intersection converges rapidly with increasing values 
of M and N as shown in Figs. 5.17 and 5.18. The total radial resultant plotted in these 
figures is defined by 

0 

F, = J [A 2r Cos9 - A g r Sind ] (Rq + e T ) d6, (5.1) 

-e 

since the resultant from the stringer vanishes by Eq. (4.61). (A more general approach to 
the resultants at the stiffener intersection is discussed in subsection 6.3.2.) From the linear 
analysis F z = —357.5 lbs, and from the nonlinear analysis F z = -382.8 lbs. Since the 
applied radial load due to internal pressure acting on the repeating unit is 1,160 lbs (= 10 
psix20 in. X 5.8 in.), the ring resists about 30.8% of this applied load in the linear response, 
and this percentage is increased to 33% in the nonlinear response. The remaining portion 
of the applied radial pressure load is carried by the shell. 





Fig. 5.18 Normal load intensity k, and total normal load F z at the stiffener inter- 
section for increasing number of harmonics in the nonlinear analysis at 10 psi. 
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of geometric nonlinearity into the analysis. The axial force carried by the stringer due 
to the closed-end pressure vessel effect is increased in the nonlinear analysis with respect 
to its value in the linear analysis. Also, the circumferential force carried by the ring is 
increased in the nonlinear analysis with respect to its value in the linear analysis. Thus, 
the stiffeners resist an increased portion of the internal pressure load, accompanied by a 
commensurate decrease in the load carried by the shell, when geometric nonlinearity is 
included into the analysis. 
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of geometric nonlinearity into the analysis. The axial force carried by the stiinger due 
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The elements of the transverse shear stiffness matrix in Eq. (2.34) are 

A 44 = A 55 = 0.69264 x 10 5 /6/m., .4 45 = 0 

The bending and stretching-bending coupling submatrices for classical model are given by 

'474.937 256.071 0 

D = 256.071 615.194 0 lb in. B = 0 

0 0 276.965 _ 

The extensional stiffness submatrix A is essentially the same for classical theory and the 
transverse shear deformation theory. 

Cross sections of the stiffeners and their dimensions are shown in Fig. 6.1. The 
stringer is an inverted hat section laminated from twelve plies of AS4/938 graphite-epoxy 
tow prepreg with a [±45. 0 2 , 90, ±15. 90, 0 2 , ±45] r lay up and total thickness of 0.0888 in. 
The stiffnesses in Hooke’s law for the stringer in Eq. (2.52) are 

(EA) S = 0.6675 x 10 7 /6. ( EI) S = 0.2141 x 10 7 /6 in. 2 , (GA) S = 0.843 x 10 6 /6 

The frame, or ring, is a 2-D braided graphite-epoxy J-section consisting of 0° and 90° 
tows. The wall thickness is 0.141 inches, and the elastic modulii are assumed to be 
Ex = 7.76 x 10 6 lb/ in 2 , E 2 = 8.02 x 10 6 /6/m. 2 , G'i 2 = G'i 3 = G 23 = 1.99 x 10 6 /6/m. 2 , and 
v \2 = 0.187. Using the ring material properties and the cross-sectional dimensions, the 
stiffness matrix for the ring in Eq. (2.58) is computed from the computer code developed 
by Woodson 29 . The non-zero stiffnesses are 

EA = 0.9088 x 10 7 /6, EI xx = 3.915 x 10 7 /6 in. 2 , EI z: = 0.1867 x 10 7 /6 in. 2 

EIzx = 0.2993 x 10 7 /6 in. 2 , EI UX = -1.322 x 10 7 /6 in 3 , GJ = 0.1346 x 10 5 /6 in. 2 
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The stiffnesses in Hooke’s law for the stringer in Eq. (2.52) are 

( EA) S = 0.6675 x 1076, (El), = 0.2141 X 10 7 /6 in. 2 , {GA) S = 0.843 x 10 6 /6 

The frame, or ring, is a 2-D braided graphite-epoxy 7-section consisting of 0° and 90° 
tows. The wall thickness is 0.141 inches, and the elastic modulii are assumed to be 
Ei = 7.76 X 10 6 lb/in. 2 , E 2 = 8.02 X 10 6 lb/in. 2 , G u = G 13 = G 23 = 1.99 X 10 6 lb/in. 2 , and 
V\ 2 = 0.187. Using the ring material properties and the cross-sectional dimensions, the 
stiffness matrix for the ring in Eq. (2.58) is computed from the computer code developed 
by Woodson 29 . The non-zero stiffnesses are 

EA = 0.9088 x 10 7 /h, EI XX = 3.915 X 1076 in. 2 , EI ZZ = 0.1867 x 1076 in 2 
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EI W „ = 1.705 x 10 7 lb in. 4 , EI U . = -0.1865 x 10 6 /6 in. 3 
GA xe = GA, e = 0.2396 x 10 7 lb 

All the results presented are for an internal pressure p = 10 psi, and the Fourier Series 
are truncated at twenty-four terms in the x- and ^-directions (M = N = 24). Based on 
M = N = 24, the transverse shear deformation model consists of a total of 6414 degress 
of freedom, and classical model consists of 3966 degress of freedom. 

6.3 INFLUENCE OF AN ASYMMETRICAL SECTION RING 

6.3.1 INTERACTING LOAD DISTRIBUTIONS 

The distributions of the interacting line load intensities between the stiffeners and 
the shell are shown in Figs. 6.2 through 6.8. The effects of transverse shear deformations 
and of warping deformation of the ring's cross section due to torsion on the magnitudes 
of the interacting line loads are summarized in Table 6.1. For the component X xs tangent 
to the stringer (Fig. 6.2), there are only small differences in the distributions as predicted 
b.V the four structural models. However, the peak value of the component normal to the 
stringer, A. 5 , is reduced in the transverse shear deformation models with respect to its 
peak value in the classical models (Fig. 6.3 and Table 6.1). 

The distributions of axial force intensity, A xr , between the ring and shell predicted 
bv the classical and shear deformation models with warping are nearly the same (Fig. 
6.4). However, the distributions of this force intensity predicted by the classical and shear 
deformation models without warping have significant differences. Thus, this interacting 
load intensity is more sensitive to the inclusion or exclusion of warping of the ring cross 
section into the structural model. As shown in Fig. 6.5, the differences in the results for 
circumferential force intensity, A^ r . between the ring and shell from the four models are 
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small, except in the vicinity of the stiffener intersection where the effects of including the 
transverse shear deformation into the models are manifested. However, the differences in 
A^ r occur over one wave length of the highest frequency i.e., A0/0 = 2/24. Differences 
occuring over the shortest wavelength may not be significant; more terms in the Fourier 
series are required to verify this. The distributions of the normal force intensity, A~ r , 
between the ring and shell predicted by the four models are essentially the same (Fig. 
6.6). The distributions of the circumferential moment component, A^ r . predicted by the 
classical models have higher magnitudes as compared to shear deformation models (Fig. 
6.7 and Table 6.1). Also note the change in sign of Ag r distributions in the vicinity of 
the joint as a result of inclusion of warping into the models. The classical theory predicts 
much larger magnitudes of normal moment component, i \ zr , compared to the transverse 
shear deformation theory for the models in which warping is included (Fig. 6.8 and Table 
6.1). However, the reverse is true for the structural models with no warping. Also, there is 
a change in sign in the distributions of A zr for classical models with and without warping 
effects. 

The distribution of the normal component of the traction across the width of the 
attachment flange of the ring is represented by line force intensity A zv and line moment 
intensity A# r . The values of A zr are nearly the same in the classical and transverse shear 
deformation models (Fig. 6.6), but magnitudes of A# r are substantially decreased in the 
transverse shear deformation models with respect to the classical models (Fig. 6.7). Thus, 
the asymmetry of the normal traction across the flange width of the ring is decreased in 
the transverse shear deformation models with respect to the classical models. 

The distribution of the circumferential component of the traction across the width of 
the attachment flange of the ring is represented by line force intensity A^ r and line moment 
intensity A- r . The values of A^ r are nearly the same in the classical and transverse shear 


126 



small, except in the vicinity of the stiffener intersection where the effects of including the 
transverse shear deformation into the models are manifested. However, the differences in 
\ 9r occur over one wave length of the highest frequency i.e., A 0/0 = 2/24. Differences 
occuring over the shortest wavelength may not be significant; more terms in the Fourier 
series are required to verify this. The distributions of the normal force intensity, X zr , 
between the ring and shell predicted by the four models are essentially the same (Fig. 
6.6). The distributions of the circumferential moment component, A $ r , predicted by the 
classical models have higher magnitudes as compared to shear deformation models (Fig. 
6.7 and Table 6.1). Also note the change in sign of Ag r distributions in the vicinity of 
the joint as a result of inclusion of warping into the models. The classical theory predicts 
much larger magnitudes of normal moment component, A-. r , compared to the transverse 
shear deformation theory for the models in which warping is included (Fig. 6.8 and Table 
6.1). However, the reverse is true for the structural models with no warping. Also, there is 
a change in sign in the distributions of A zr for classical models with and without warping 
effects. 

The distribution of the normal component of the traction across the width of the 
attachment flange of the ring is represented by line force intensity X zr and line moment 
intensity A #r . The values of X~ r are nearly the same in the classical and transverse shear 
deformation models (Fig. 6.6), but magnitudes of Ag r are substantially decreased in the 
transverse shear deformation models with respect to the classical models (Fig. 6.7). Thus, 
the asymmetry of the normal traction across the flange width of the ring is decreased in 
the transverse shear deformation models with respect to the classical models. 

The distribution of the circumferential component of the traction across the width of 
the attachment flange of the ring is represented by line force intensity X$ T and line moment 
intensity A zr . The values of X 6r are nearly the same in the classical and transverse shear 






deformation models (Fig. 6.5). However, the magnitude of A zr is substantially increased 
in the transverse shear deformation model with respect to the classical model with warping 
excluded, and is substantially decreased in the transverse shear deformation model with 
respect to the classical model with warping included (Fig. 6.8). Thus, the asymmetry 
of the circumferential traction across the flange width of the ring is increased in the 
transverse shear deformation model with respect to the classical model without warping, 
and is decreased in the transverse shear deformation model with respect to the classical 
model with warping. 

The inclusion of transverse shear deformation and warping of ring's cross section 
into the analyses influences the distributions and magnitudes of interacting line load com- 
ponents A, s , X xr , A# r , A Sr* and A zr . The distributions of interacting line load components 
A X5 and A~ r remain essentially the same. The cause of sensitivity to transverse shear de- 
formations is two- fold: First, the tangential displacements of the shell along the contact 
lines are de-coupled from the out-of-plane rotations of the reference surface of the shell, 
and for the stiffeners the longitudinal displacements along the contact lines are de-coitpled 
from the rotations of the longitudinal reference axes. Second, in the transverse shear 
deformation model, the torsional rotation of the ring at the shell-stringer-ring joint is de- 
coupled from the in-plane bending rotation of the stringer at the joint, thereby allowing 
for increased joint flexibility. In the classical model, the torsional rotation of the ring at 
the joint is constrained to be the same as the bending rotation of the stringer (see Fig. 
6.9). The values of these joint rotations for the four structural models are given in Table 
6.2. Notice from Table 6.2 that the sense of the rotation changes if warping is included, 
and that the transverse shear deformation results in a torsional rotation of the ring that 
is about twice as much as the bending rotation of the stringer. 
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Table 6.2 Rotations About The Circumferential Axis At The Stiffener Intersection. 


Rotations in 10' 5 radians 


Description of the Rotation of 
the Structural Component 

Classical Theory 

Transverse Shear Theory 

No Warping 

Warping 

No Warping 

Warping 

Shell normal (0, 0) 

-2.56 

2.58 

-1.06 

2.65 

Ring twist <|> er (0) 

-2.56 

2.58 

-2.67 

3.64 

Stringer normal <J> 0 ^ ( 0) 

-2.56 

2.58 

-0.29 

1.85 
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Table 6.3 Resultants At Stiffener Intersection. 


Components of the 

Classical Theory 

Transverse Shear Theory 

Resultant 

No Warping 

Warping 

No Warping 

Warping 

Cq from stringer, 
lb-in. 

- 1.1696 

-0.0921 

- 0.2953 

- 0.7797 

C e from ring, lb-in. 

1.627 

5.645 

1.363 

6.0192 

Cq total, lb-in. 

0.457 

5.5526 

1.0676 

5.2396 

F z , lb. 

- 564.06 

- 564.56 

- 562.27 

-563.15 


136 



Table 6.3 Resultants At Stiffener Intersection. 


Components of the 
Resultant 

Classical Theory 

Transverse Shear Theory 

No Warping 

Warping 

No Warping 

Warping 

C 0 from stringer, 

- 1.1696 

-0.0921 

- 0.2953 

- 0.7797 

lb-in. 





C 0 from ring, lb-in. 

1.627 

5.645 

1.363 

6.0192 

C 0 total, lb-in. 

0.457 

5.5526 

1.0676 

5.2396 

F v lb. 

- 564.06 

- 564.56 

- 562.27 

-563.15 


136 



V-Mu :i L '-'' " 


x 10' 1 (lb) or M (lb in) 

0S 



x/l 


Fig. 6.11 Stringer axial force and bending moment distribution 
at 10 psi. 
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Fig. 6.11 Stringer axial force and bending moment distribution 
at 10 psi. 



138 


asymmetric about the origin. Only small differences are predicted by the four structural 
models in the distribution of V zs . 

The distributions of the circumferential force and in-plane bending moment in the 
ring are shown in Fig. 6.13. The differences in these distributions predicted by the four 
models are very small. The distributions of the in-plane shear force, V- r , in the ring 
predicted by the four structural models have negligible differences, as shown in Fig. 6.14. 
The out-of~pla.ne bending moment M zr and torque T r in the ring are more sensitive to 
the change in models as shown Fig. 6.15. The distributions of the out-of-plane bending 
moment are symmetric about the origin, and their magnitudes predicted by the models 
with warping included are substantially larger as compared to their magnitudes predicted 
by the models without warping. The distributions of total torque, T r ( = T sr — M^r/Ro ), 
are antisymmetric about the origin. As shown in Fig. 6.15, the torque has reduced 
magnitudes in the transverse shear deformation model compared to the classical model 
when warping is included. The torque predicted by the models without warping is St. 
Yen ant’s torque T sr , and this is negligible as shown in Fig. 6.15. The distributions of out- 
of-plane shear force, V xr , in the ring are shown in Fig. 6.16, and these distributions are 
antisymmetric about the origin. The magnitudes of V xr predicted by the transverse shear 
deformation model are larger compared to the classical model when warping is included. 
However the reverse is true for the V xr distributions without warping. The distributions 
for M# s , A/ cr , T r and V xr shown in Figs. 6.11, 6.15 and 6.16, respectively, indicate that 
these stiffener actions are sensitive to both transverse shear deformations and warping 
deformations. 

6.3.5 SHELL RESPONSE 

The distribution of the normal displacement along a-curve midway between the 
stringers (8 — —0), and along the 0-curve midway between the rings (a- = -/), are 
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Fig. 6.14 Ring in-plane shear force distribution at 10 psi. 
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Fig. 6.16 Ring out-of-plane shear force distribution at 10 psi. 
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shown in Fig. 6.17. As depicted in this figure, there is a negligible difference between 
the results from the transverse shear deformation model and classical model (warping of 
the ring is included in both models). Also, there is negligible difference in the axial and 
circumferential normal strain distributions between the two models as shown in Figs. 6.18 
and 6.19. Thus, the normal displacement and in-plane normal strains for the shell are not 
significantly influenced by the inclusion of either transverse shear deformations or warping 
deformation of the ring into the structural models, in part because the shell is very thin 
for the example studied {R/t = 1268). 

6.4 A RING WITH SYMMETRICAL CROSS SECTION 

As a benchmark for comparing the transverse shear deformation model with the 
classical model, analyses were performed for a ring with symmetric cross section. In this 
case the changes made to the numerical example under discussion are to set the bending- 
coupling stiffeness EI, X , the out-of-plane bending-to- warping coupling stiffness EI UX , and 
the contour warping function parameter u>o of the ring, all to zero. Consequently, the 
0-axis, as well as the z-axis, are axes of symmetry for the repeating unit in terms of 
geometry, load, and material properties. Symmetry about the 0-axis implies there is no 

out-of-plane bending and torsion of the ring; i.e., u r (0) = <?W(0) = <Pzr{0) = A xr (0) = 
\ 6r (0) - A zr (0) = 0 for -0 < 0 < 0. Thus, for the symmetrical section stiffeners only 

the interacting line load components tangent and normal to the stiffeners are non-zero. 
Since there is no torsion, warping of the ring cross section does not play any role in the 
analyses. 

The distributions of the tangential interacting load intensity between the shell and 
ring are shown in Fig. 6.20. The differences in the results from the two models are small 
except in the vicinity of the stiffener intersection. The peak magnitude of \ 6r in the 
transverse shear deformation model is smaller than the peak value for X 6r in the classical 
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Fig. 6.19 Distributions of the circumferential normal strain on the inner 
and outer shell surfaces midway between the rings (x=-l) at 10 psl 





model (50.8 lb /in versus 64.5 lb/in.). However, this difference occurs over one wave 
length of the highest harmonic retained in the analysis, and may not be significant. The 
distributions of the tangential and normal interacting load intensities between the shell 
and stringer, and the normal load intensity between the shell and ring are not significantly 
different in the two models. 

For a symmetrical cross section ring, in Eqs. (6.1) through (6.6) F x = F e = C x = 
Cfi = C- = 0. The only non-zero component of the force resultant is the radial force F z . 
The values of F z computed from the classical and transverse shear deformation models 
are -563.72 lb. and -561.89 lb., respectively. 

6.5 SUMMARY OF RESULTS 

The asymmetrical section ring complicated the analysis of the unit cell, since sym- 
metry about the plane of the ring is lost. Out-of-plane bending and torsion of the ring 
occur as well as a rotation of the shell-stringer-ring joint, about the circumferential axis. 
Inclusion of transverse shear deformations, and warping deformation of the ring s cross 
section due to torsion, into the mathematical model significantly influenced several aspects 
of the response. 

The sense of the rotations of the structural elements at the joint is changed with the 
inclusion of warping deformation in the ring, and the twist rotation of the ling at the joint 
increases by 40% with the inclusion of transverse shear deformation, as is shown in Table 
6.2. That is, joint flexibility increases since element rotations at the joint are de-coupled 
by using transverse shear deformation models. 

The interacting loads that are strongly affected by the structural modeling are those 
components associated with the asymmetric response. These are the axial force intensity 
A rr , the normal moment intensity A zr , and the tangential moment intensity A$ r > between 
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the ring and shell. At the joint the magnitude of X xr is increased by the inclusion of 
warping and reduced by the inclusion of transverse shear (Fig. 6.4). The normal moment 
intensity A zr is a measure of the asymmetry in the distribution of the circumferential 
traction across the width of the ring’s attachment flange. The sense of A zr is changed 
when both transverse shear and warping are included, and its magnitudes are reduced 
by the transverse shear effect (Fig. 6.8). However, the magnitude of the resultant of 
the circumferential traction across the flange width of the ring, as measured by the line 
load intensity A# r , is only moderately affected by the changes in the structural models 
(Fig. 6.5). The tangential moment intensity A^ r is a measure of the asymmetry in the 
distribution of the normal traction across the width of the ring’s attachment flange. At the 
joint, the sense of A# r is changed by the inclusion of warping deformation and additionally 
its magnitude is reduced by the inclusion of transverse shear deformation (Fig. 6.7). 
However, the magnitude of the resultant of the normal traction across the flange width, 
as measured by the line load intensity \ zr , is essentially unaffected by changes in the 
structural models (Fig. 6.6). 

The distributions of the normal actions between the shell and stiffeners (A Z5 , A 3r , and 
A# r ) all show significant magnitudes only in the vicinity of the joint, with much smaller 
magnitudes away from the joint. In fact, they all appear to exhibit singular behavior at the 
joint, but only finite magnitudes are represented by the truncated series approximations. 
The distributions of the actions tangent to the stiffeners (A xs , A^ r , and A 3r ), on the other 
hand, have small magnitudes in the vicinity of the joint and larger magnitudes away from 
the joint. These tangential actions do not exhibit singular behavior. (These results are 
similar to those found for the symmetric stiffeners problem discussed in Chapter 5.) 

In spite of the singular behavior of the line load intensities associated with the normal 
actions, the resultant of these distributions resolved at the joint converge relatively quickly 
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with the number of terms retained in the series approximations. The resultant consists 
of a radial force F z and a moment Co about the circumeferential axis. F- represents the 
action of the stiffeners to pull the shell radially inward against the pressure load, and 
Co is primarily due to asymmetry in the actions between the ring and shell. Foice F z is 
essentially unaffected by the structural modeling (Table 6.3). and its magnitude for the 
example studied is about 17% of the total pressure load carried by the unit cell. The 
remaining pressure load is carried by the shell itself. The moment Co is very sensitive to 
the structural modeling, in particular to the effect of warping as sliowui in Table 6.3. This 
moment Co vanishes for a completely symmetric problem. 

The out-of-plane bending moment and torque in the ring are very sensitive to the 
structural modeling, as might be expected. The magnitudes of the both the out-of-plane 
bending moment and torque increase with the inclusion of warping and transverse shear 
into the mathematical model (Fig. 6.15). However, the shell’s normal displacement and 
strains are insensitive to the changes in the structural models for the very thin shell 
(R/t = 1268) used in the numerical example. 
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CHAPTER 7 

CONCLUDING REMARKS 


7.1 SUMMARY 

Structural analyses are developed to determine the linear elastic and geometrically 
nonlinear elastic response of an internally pressurized, orthogonally stiffened cylindrical 
shell. The structural configuration is of a long circular cylindrical shell stiffened on the 
inside by a regular arrangement of identical stringers and identical rings. Periodicity of 
this configuration permits the analysis of a unit cell model consisting of a portion of the 
shell wall centered over one stringer- ring joint; i.e., deformation of a structural unit cell 
determines the deformation of the entire stiffened shell. See Fig. 1.1. The stringer-ring- 
shell joint is modeled in an idealized manner; the stiffeners are mathematically permitted 
to pass through one another without contact, but do interact indirectly through their 
mutual contact with the shell at the joint. The stiffeners are modeled as discrete beams. 
The stringer is assumed to have a symmetrical cross section and the ring either a symmetric 
or an asymmetric open section. 

The formulations presented for the linear elastic response include the effect of trans- 
verse shear deformations and the effect of warping of the ring’s cross section due to torsion. 
These effects are important when the ring has an asymmetrical cross section, because the 
loss of symmetry in the problem results in torsion of the ring, as well as out-of-plane 
bending, and a concomitant rotation of the joint at the stiffener intersection about the 
circumferential axis. Restraint of cross-sectional warping, as occurs here in the ring due 
to contact with the shell, is an important contributor to the normal stresses in thin- walled 
open section bars. Based on transverse shear deformations and cross-sectional warping 
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of the ring, four structural models are defined. The simplest model uses non-transverse- 
shear- deformable theory, or classical theory, and neglects warping due to torsion. The 
most complex model includes both effects. Models of intermediate complexity occur for 
inclusion of one effect without the other. The formulations presented for the geometrically 
nonlinear response take into consideration only the symmetrical cross section of stiffeners, 
and are based on classical structural theories. 

For all the structural models, the response of the unit cell under internal pressure is 
obtained by utilizing the Ritz method. Displacements are assumed as truncated Fourier 
Series plus simple terms in the axial coordinate to account for the closed-end pressure 
vessel effect. Equilibrium is imposed by virtual work. Pointwise displacement continuity 
between the shell and stiffeners is achieved by Lagrange multipliers which represent the 
interacting line load distributions between the stiffeners and the inside shell wall (see 
Fig. 1.2). Data from a composite material crown panel typical of a large transport 
fuselage are used for two numerical examples. The first numerical example is used to 
validate the structural model, and also to compare the linear response and geometrically 
nonlinear response of the unit cell model with symmetrical section stiffeners. In the second 
numerical example the linear elastic response of the unit cell model with an asymmetrical 
cross section ring is analyzed. 

7.2 CONCLUDING REMARKS 

7.2.1 EFFECT OF GEOMETRIC NONLINEARITY 

It is found that the spatial distribution of the normal displacements of the cylindrical 
shell are more uniform, and the bending strains are reduced, in the geometrically nonlinear 
elastic analysis with respect to what is predicted by the linear elastic analysis. That is, 
pillowing of the skin is reduced by the inclusion of geometric nonlinearity into the analysis. 
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However, in the nonlinear analysis the most severe circumferential bending is confined to a 
narrow boundary layer adjacent to the stringer midway between the rings, and the interior 
portion of the shell behaves as a membrane. In the linear analysis the bending occurs over 
the entire circumference of the shell. The development of the bending boundary layer 
due to the inclusion of geometric nonlinearity into the analysis causes an increase in the 
circumferential transverse shear resultant in the shell adjacent to the stringer compared 
to the linear analysis. Increased interlaminar shear stresses can be expected as a result of 
the increased transverse shear resultant. 

The axial force carried by the stringer due to the closed-end pressure vessel effect is 
increased in the nonlinear analysis with respect to its value in the linear analysis. Also, 
the circumferential force carried by the ring is increased in the nonlinear analysis with 
respect to its value in the linear analysis. Thus, the stiffeners resist an increased portion 
of the internal pressure load, accompanied by a commensurate decrease in the load carried 
by the shell, when geometric nonlinearity is included into the analysis. 

7.2.2 INFLUENCE OF AN ASYMMETRICAL SECTION RING 

The asymmetrical section ring complicated the analysis of the unit cell, since sym- 
metry about the plane of the ring is lost. Inclusion of transverse shear deformations, 
and warping deformation of the ring’s cross section due to torsion, into the mathematical 
model significantly influenced several aspects of the response. The interacting loads that 
are strongly affected by the structural modeling are those components associated with 
the asymmetric response. These are the axial force intensity A rr , the normal moment 
intensity A, r , and the tangential moment intensity Ag r between the ring and shell. At the 
joint the magnitude of A xr is increased by the inclusion of warping and reduced by the 
inclusion of transverse shear (Fig. 6.4). The normal moment intensity A sr is a measuie 
of the asymmetry in the distribution of the circumferential traction across the width of 
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the ring’s attachment flange. The sense of A zr is changed when both transverse shear and 
warping are included, and its magnitudes are reduced by the transverse shear effect (Fig. 
6.8). However, the magnitude of the resultant of the circumferential traction across the 
flange width of the ring, as measured by the line load intensity \» r , is only moderately 
affected by the changes in the structural models (Fig. 6.5). The tangential moment inten- 
sity A (j r is a measure of the asymmetry in the distribution of the normal traction across 
the width of the ring’s attachment flange. At the joint, the sense of A» r is changed by the 
inclusion of warping deformation and additionally its magnitude is reduced by the inclu- 
sion of transverse shear deformation (Fig. 6.7). However, the magnitude of the resultant 
of the normal traction across the flange width, as measured by the line load intensity A jr , 
is essentially unaffected by changes in the structural models (Fig. 6.6). 

The sense of the rotations of the structural elements at the joint is changed with the 
inclusion of warping deformation in the ring, and the twist rotation of the ring at the joint 
increases by 40% with the inclusion of transverse shear deformation, as is shown in Table 
6.2. That is, joint flexibility increases since element rotations at the joint are de-coupled 
by using transverse shear deformation models. 

The out-of-plane bending moment and torque in the ring are very sensitive to the 
structural modeling, as might be expected. The magnitudes of the both the out-of- plane 
bending moment and torque increase with the inclusion of warping and transverse shear 
into the mathematical model (Fig. 6.15). However, the distributions and magnitudes 
of the normal displacement and strains of the shell midway between the stiffeners are 
unaffected by the change in the structural models. 
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7.2.3 SINGULARITY AT THE SHELL-STRINGER-RING JOINT 


The distributions of the normal actions between the shell and stiffeners (Aj S , A- r , and 
hor) all show significant magnitudes only in the vicinity of the joint, with much smaller 
magnitudes away from the joint. In fact, they all appear to exhibit singular behavior at the 
joint, but only finite magnitudes are represented by the truncated series approximations. 
The distributions of the actions tangent to the stiffeners (A xs , \e r , and A zr ), on the other 
hand, have small magnitudes in the vicinity of the joint and larger magnitudes away 
from the joint. These tangential actions do not exhibit singular behavior. In spite of 
the singular behavior of the line load intensities associated with the normal actions, the 
resultant of these distributions resolved at the joint converge relatively quickly with the 
number of terms retained in the series approximations (e.g., see Fig. 5.17). The resultant 
consists of a radial force F z and a moment Ce about the circumferential axis (see Fig. 
6.10). F z represents the action of the stiffeners to pull the shell radially inward against 
the pressure load, and Ce is primarily due to asymmetry in the actions between the ring 
and shell. The magnitude of force F z represents the portion of the total pressure load 
carried by the stiffeners. The remaining pressure load is carried by the shell itself. Force 
F z is essentially unaffected by the structural modeling (Table 6.3). The moment C$ is 
very sensitive to the structural modeling, in particular to the effect of warping as shown 
in Table 6.3. This moment Ce vanishes for a completely symmetric problem. 

The series for the interacting normal load intensity X zr at the stiffener intersection 
does not appear to converge even in the geometrically nonlinear analysis. However, the 
Fourier Series for the total radial resultant load carried by the stiffeners, which is resolved 
at the intersection, exhibits rapid convergence in the geometrically nonlinear analysis (see 
Fig. 5.18). The total radial resultant load carried by the stiffeners is slightly increased in 
the geometrically nonlinear analysis with respect to its value in the linear analysis. 
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7.3 RECOMMENDATIONS FOR FUTURE WORK 

The analysis could be extended to include the important load cases of torsion and 
bending in addition to internal pressure. Torsion of the stiffened shell is equivalent to 
a shear load of the unit cell, and this loading case corresponds to an antisymmetric de- 
formation pattern of the unit cell model. Bending is somewhat more complex since the 
displacement field has period of 2tt rather than the periodicity of the stringer spacing. 

Singular solutions can be investigated to improve the solution methodology. To 
begin with, the simplest configuration can be analyzed; i.e., consider the linear response 
of a classical structural model with symmetric section stiffeners. This problem can be 
reformulated in terms of four integral equations in which the unknowns are the magnitudes 
of the interacting load components A IS (z), A- S (ar), \g r (8), and A zr (8). These integral 
equations are associated with either the tangential or normal displacement constraints 
between the shell and stiffeners. The kernels to these equations are Green’s functions 
(or influence functions), which are displacement solutions to the shell, stringer, and ring 
problems under the actions of concentrated forces. This approach is appealing because 
the solution reduces to solving for the interacting load components in integral expressions 
along the contact lines rather than solving for the displacements and interacting loads 
over the entire solution domain. The difficulty with this approach is in obtaining the 
Green's functions, and solving singular integral equations that may result from the contact 
problems. 
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NOMENCLATURE 
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Laminate extensional stiffness 

Laminate extension-bending coupling stiffness 

Laminate bending stiffness 

Eccentricity of ring 

Eccentricity of stringer 

Three-dimensional engineering normal strains 

Three-dimensional engineering shear strains 

Harmonic in axial direction 

Number of harmonics in axial direction 

Ring bimoment 

Stringer bending moment 

Ring in-plane bending moment 

Stress couples of shell theory 

Stress couple of Sanders’ theory of shells 

Modified stress couples of transverse shear 

deformation shell theory 

Ring out-of-plane bending moment 

Harmonic in circumferential direction 

Number of harmonics in circumferential direction 

Ring circumferential force 

Stringer axial force 

Stress resultants of shell theory 

Stress resultant of Sanders’ theory of shells 

Internal pressure load 
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Stringer axial load share 
Transformed reduced stiffness of lamina 
Transverse shear resultants for shell 
Radius to middle surface of shell 
Radius to the reference arc of ring 
Area of the shell’s reference surface 
Thickness of shell 
Stringer spacing 
Frame, or ring, spacing 
Torque in the ring 

Axial displacement of middle surface of shell 
Axial displacement of reference arc of ring 
Axial displacement of centroidal axis of stringer 
Rigid body displacements 

Circumferential displacement of middle surface of shell 

Circumferential displacement of reference arc of ring 

Stringer shear force 

Ring in-plane shear force 

Ring out-of- plane shear force 

Virtual work 

Normal displacement of middle surface of shell 

Normal displacement of reference arc of ring 

Normal displacement of centroidal axis of stringer 

Axial frequency 

Circumferential frequency 

Normal strains of ring’s reference axis 
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Normal stress components 
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Ring’s transverse shear strains 

Stringer’s transverse shear strain 
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APPENDIX A 


ELEMENTS OF MATRICES FOR LINEAR ANALYSIS USING 
TRANSVERSE SHEAR DEFORMATION MODEL 

The non-zero elements of the submatrices B{j and An in Eq. (4.68) for the trans- 
verse shear deformation model are listed below. The parameter is Kronecker delta assuming 
the values zero for i ^ j, and one for i = j, respectively. 

Elements of (10 MN + 3M + 3A + 2) x (10 MN + 3M + 3N + 2) submatrix [A'u] 


A'n (1,1) 
A'u (1,2) 
A'n (2,2) 
A'u (3,3) 
An (3,4) 
An (3,5) 
An (4,4) 
An (4,5) 
An (5,5) 

A n (6,6) 
An (6,7) 
An (6,8) 
An (7,7) 
An (7,8) 


An Q R 
'' / 

= 2 .4 1 2 0 A 
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=2.4 n a m lQ Rb rn p 
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--2Biior m lQR6 mp 
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^22 n2 i ^55 \](^\ £ 

.- r /?n + - s -V 0 «- 


=2( ^ + ^ )/?n/0 ^ 

=2(^0l ~ Ass)lQ6 n q 

=2(^ + ^0l)lQS nq 

= 2 (%- - A^)(5 n m nq 
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An (8,8) =2(^/3* + A^R)lQb nq 

A'n (9,9) =[A n a^A+ ^/?j[]/0tfmp^», = A'n ( 10 , 10 ) 

A u (9,12) =-[Ai 2 + A 66 + ^ 66 — ^ 66) ]a m /? w /0^ mp ^nq = - A'n (10,11) 


A 1 1 (9,14) — — <4-12 O^m&^mp^nq — 011 (10,13) 

A'n (9,15) =[5n« J m + ( ^ 66 = A'n (10,16) 

An (9, 18) = — (.012 + -066 4~ 066 )^m/3n^0^mp^ng — “011 (10, 17) 

r - in f ^22 ^2 , ^55 . r , . 2(^66 - 5| 6 ) ( #66 ~ 2 # 66 + AkOi,^ ~1 ; 
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Elements of (6 M + 1) X (6 M + 1) submatrix [Ii' 22 ] 
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Elements of (lOMA + 3M + 31V + 2) X 4M submatrix [An] 
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Elements of (10 MN + 3 M + 3N + 2) x (57V + 1) submatrix [Bn] 
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fii2 (8,3) —~ciG6 nq 
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B \ o (9,2) — —aQfinq 
B \ 2 (11,3) — — 

Bi 2 (12,6) =-a p a0<5 n9 
5i2 (13,4) =-a0^ ng 
5)2 (14,5) = CVp(Z0^ng 
fii2 (15, 2) = -a06 n g 
B n (17,3) =-aQ6 nq 

B\2 ( 18, 6) = —Otp(lQ& n q 

Elements of (10 MN + 3M + 3N + 2) x 1 submatrix [£ 13 ] 


5 13 (U) = 1 


Elements of (6 M + l)x4M submatrix [B 2 i] 


(- 1 ) 


m +1 


£21 (M)= a 
^21 ( 2 , 1 ) = lfimp 
B 2 I (3,2) = l&mp 
B2\ (4,3 ) — 
i?21 (3, 4) —l^mp 
B 21 (6,1 )=le s 6 mp 
B 2 1 (7,2) =/e,* mp 
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!! 


Elements of (6 M + 1) x 1 submatrix [B 03 ] 


B23 ( 1 , 1 ) = -1 


Elements of (67V + 1) x (57V + 1) submatrix [$ 32 ] 


B32 (1,1) =2a0 
B 32 (2,2) =a06 nq 
B 32 (3,3) =aQS nq 
B 32 (4,4) =uQ6 nq 


B 32 (5,2) — c 7*fl0 S ji q 


i (Ja 

B 32 (5,3) = -^/3 n aQ6 nq 

no 


B 32 (5,5) ~aQS 


nq 


B 32 (5,6) =—(3 n a06 nq 
-< 10 


B 32 (6,3) =e r aQ8 


nq 


B 32 (7,3) =~aQS nq 
no 

B 32 (7,6) =(1- ^-)a0tf n , 
no 


Elements of (10M7V + 3 M + 37V + 2) x 1 submatrix [F\\] 


F u (1 )=pR 2 Q 
F n (2) =4plQR 
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APPENDIX B 


ELEMENTS OF MATRICES FOR LINEAR ANALYSIS 
USING CLASSICAL MODEL 

The non-zero elements of the submatrices Kij, B{j and F\\ in Eq. (4.68) for the classical 
model are listed below. The parameter 6ij is Kronecker delta assuming the values zero for 
i ^ j, and one for i = j, respectively. 

Elements of (6MN + 2 M + 2N + 2) x (6 MN + 2 M + 2N + 2) submatrix [A'n] 


An (1,1) = 


A U QR 

l 


A'n (1,2) =2j4i2©A 
4 Aoo Ql 


An (2,2) = 


R 


A'n (3,3) —2A\\a 2 m lQRb mp 
An (3,4) =2( - + Bno:^ n )o' m lQ RS m p 

A'n (4,4)=2(^+'^a 2 m + D 11 a i m )lQR6 mp 
A'n (5,5) =2(^ + ^|f- + ^)!3 2 n lQ6 nq 
An (5,6)=2[(^ + ^) + (^ + ^]^0A 


R 2 


nq 


A'n (6,6) = 2 {~ + 

All (7,7) =[^4ll a mA + (-R p 2 ~ + )/^n] = An (8,8) 


A 2 4A 3 


A'n (7,10) =-[^ + + (^ + §r - ^)]«m/3n/e/M«,tfn, = - A'n (8,9) 


A 2 


A 


A 2 4 A 3 


A'n (7,12) =-[A 12 + B\\or m + ^ - ^)/£]a m 0/Wn, = - A'n (8,11 
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nq 


A-„ (9,9) =[(^ + Hf* + ^)/£ + (A M + fliew 

=A'n (10,10) 

r ^22 Bio , Bio D 22 \ /i2 / E\2 . *2 .Bqq D]2 . 3Z1$$ , 2 1 
A'n(9,ii)=[-f + i + (i + ir)^ + (T + X + ¥ + ^ Kl * 

p n lQR6 mp 6 nq = A'n (10,12) 

A'n (11,11) =[Al + ^5 + + 2(Pl ^ 2Cw) <# + (2B« + °n<R> 

a 2 m ]m mp S nq = A'n (12,12) 


Elements of (AM + 1) x (4 M + 1) submatrix [A' 22 j 


A 22 (1,1) = 


(AA) S 

21 


A22 (2,2) =(£4) s aJ,M mp 
A 22 (3,3) =(EA) s a m l8 mp 
A 22 (4,4) ={EIee)sO A m l8 mv 
A 22 (5,5) = {EIgg) s ot m l6 m p 

Elements of (4JV + 1) X (41V + 1) submatrix [A' 33 ] 


A' 33 (1,1) = 
A' 33 (2,2) = 


2EAQ 

R 0 

GJ + (£4.+ ^ i + E! ““'“ 2 


Ro 


Rl 


’)/*» 


^es nq 


A' 33 (2,3) =-(£/„ + ^ £ )4©<5n, ) 

ito ^0 

A' 33 (2,4) = -(£/« + 

A' 33 (2,5) =-f(GJ + £/„ + 


^ ) + ( 


Ao 


EIu/z EIuu 2 


R 0 


Rl 


)Pl 


a 2 

%Q8 nq 
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A _ 33 (3,3) -(EA + ^ 2 '")^® ^”9 
A '33 (3,4) =(EA H — ©^ng 


Ri 

EL 


A ’33 (3,5) =(EI ZX + /?n)-£ 2 @^nij 


A '33 (4,4) ={EA 4 


A 33 (4,5)=(A/„ + ^^)§0^ 


Ao 
EA 
R 0 


&_ 
R 0 


A 33 (5,5) = 


E/» + (GJ + ^^+ f:/ “" rf '' a2 


Ao 


A 5 


PM 


©c 

A 0 ” 9 


Elements of ( 6 MiV + 2 M + 2N + 2) X 4M submatrix [A„] 


An 


( 1 , 1 ) 


(-l) m 


B n (3,1) =— W, 


(4,1) 


-Bn (4,4) — • mmm l6 r np 

B\l (7, 2 ) — — 

B\\ (8,1) = “^mp 
x t 

B\\ (11,1) — — Q m t A 1 m p 
5n (11,4) =-W m p 

5n (12,2) =-"^rri^mp 

i?ll (12,3) — “ l&mp 
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Elements of (6 MN + 2 M + 2N + 2) x (5JV + 1) submatrix [B n] 

5i2 (2,1) =-2a0 
5i2 (4,1) =— 2a0 
Bn (5,3) =— (1 — — )a0^ n <j 
Bn (6,3) 

Bn (6,4) = uQ6 n q 
Bn (7,2) =—aQ6 nq 
Bn (9,3)=-(l-^)a0«n, 

B u (10,6) =-(l - ^)a p a©6 n? 

Bn (11,3) =^f3 n aQ6 nq 
Bn (11,4) ——(iQS n q 
B \ 2 (12,2) — — ttpfl0^n^ 

(12,5) = dpdQf>nq 
Bn ( 12 , 6 ) = ^J^Pn a p a Qfinq 

Elements of (6M1V + 2M + 21V + 2) X 1 submatrix [5 i 3 ] 

Bn (1,1) = 1 

Elements of (4 M + 1) X 1 submatrix [5 23 ] 

5 23 (1,1)= -1 
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Elements of (AM + 1) x AM submatrix [B 2 1 ] 


Bn ( 1 , 1 ) = 


(- 1 ) 


m+l 


Ckr 


f?21 (2,1) =lftmp 

Bn (3,2) =lS mp 
B 21 (4, 2) = o TO /e s $ mp 
£21 (4,3) =l6 mp 
B 21 (5, 1) = o: m l6 s S m p 
B 2 \ (5,4) =/^ mp 


Elements of (4iV + 1) x (5 N + 1) submatrix [B 32 ] 


B 32 (1,1) =2 a& 
B 32 (2,2) =aQ 6 


r nq 


B 3 2 (2,3 )=-^(3 n aQ6 nq 


R\ 

B 32 (2,6) = (1 

Ko 

B 32 (3,3) =(1 + %-)aG 6 ng 

JtO 

B 31 (4,3) =— Pn a Qftnq 
Ko 

B 32 (4,4) =aQ6 nq 

B32 ( 5 , 2 ) = € r dQ6 n q 

B 32 (5,3) =^-(1- ^-)f3 n a06 nq 
Ko Ko 

B 32 (5,5) =a& 6 nq 

B 32 (5,6) =[-^ + (1 -J-)YV naQb 


nq 
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Elements of (6MN + 2 M + 2 N + 2) X 1 submatrix [Fn] 


F n (1) =pR 2 Q 
Fn (2) =4 plQR 
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APPENDIX C 


ELEMENTS OF TANGENT STIFFNESS AND LOAD STIFFNESS 
MATRICES FOR NONLINEAR ANALYSIS 


The non-zero elements of the tangent stiffness submatrices Ajj and load stiffness sub- 
matrix L in Eq. (4.82) for nonlinear analysis are listed below. The parameter 6{j is Kronecker 
delta assuming the values zero for i ± j, and one for i = j, respectively. 

Elements of (3 MN + 2 M + 2 N + 2) x (3 MN + 2M + 2N + 2) submatrix [A'n (««&«!<)] 


A'n (1,1) 


A U QR 

l 


A'n (1,2) =2 Ai 2 0 
An (1,3) =0 


An 

A'n 

A'n 

An 

A'n 


A'n 


A'n 


(1.4) — A\\Q Rot‘ m w m £) jm 

(1.5) =^~(v n + I3 n w n )6 k n 

K 

(1.6) = ^ p — Pn(v n + fin W n )hn 

it 


: 1,7) =o 


(1,8) = (Vmn + fin w mn)fijmf>kn 

IK 


(1,9) = 

( 2 , 2 ) = 


A «W -0n(Vmn + fin w mn)^jm^kn + §.5Ai\QRafi m W mn l)j m f)kn 

2R 

4A 22 Ql 

R 


A'n (2,3) =0 

An (2,4) — ZA\2Qlot^ m w m b j m 

A'n (2,5) = 2y4 p, — (Pn + finWnfihn 

A'n (2,6) = 2A ^ 2 Q -fi n (v n + fi n w n )hn 
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An (2,7) 
A'li ( 2 , 8 ) 
An (2,9) 
An (3,3) 
An (3,4) 

An (3,5) 
An (3,6) 
An (3,7) 
An (3,8) 

An (3,9) 
An (4,4) 


=0 

j 4 o 2 0 / 

R 2 

A '22 0 ^ 


(^mn ^ mn )Sjm^kr 


" Pn(Vmn + flnWmn)^ jm^ kn + ^12 0^ a m w mn^jm^kn 

R z 


- 2 Aiia m QlR 6 m p 

j 

z 2A\20( m Ql$mp 4- 2A\\OL m OLpOtjQRwjl2 
j=i 

A\2 0 / 
i? ~ 

^ 12 ©/ 


R 


^p(^pn 4" jp6 kn 

a p 0n( Vpn “f” /3n^ ? pn 


= 0 


A 12 0/ 


a m (v n + /3n^n)^mp<5jtn + ^ ~^~ a P^ V J n + Pn™jn)fanll 

j= i 


^ 12 ^-a m /J n (v n + finWn)&mpbkn + ^ ^“ a p/^n(^jn + fln Wjn)bknh 

j=l 


i! 2 ' 


u mp 


+ ^ Ana m a p aj 0 .ftWjn< 5 /cn 7 2 

j=l 

S 2 ( 4 rr + Ona 4 ra ) 6 i^m P + 2 ^i 2 a 2 m 0 /wo^m P + Ano^©^, 

R~ 

K J 

4- -|- ^ / ^ > 66 )Q/a^(ufc + flk w k)~ &mp + 2 An 0 Ra m aja p UjIa 

2 " R >= l 

+ £ 2A 12 0a m W / 2 + £ Dnf + ^) 0ama > fc + ^w k )v jk h 


3 - 1 
J 


J = 1 *=1 


j=l 
J K 


J K 2 
+ ^ 3Aii0i?a 2 a p a m tu 2 / 2 3 + ^ 1.5AnQRotjOt m a p Wj k l23 

j = i *=i 


+ V V + ^)0a m a p (t; ifc + P k w jk ?Iu + £ 2A a2 0a m a p ii» i /« 

— J £ J 7 R R i i 


j = 1 *=1 
J I< 


j= 1 
J 


+ V £(^i + ^^)0a m a p /3fc(uA: + dkWk)Wjkh + £ 2 A 12 © QpQj 
j=lfc=l " R 1 =1 
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A’ll (4,5) = (4ll + — ^—)9l a2 m( V q + Pq'W q )w m 6 kq 6j m + AxoQloilnPqWmq&kqb jm 

R R 

- ^j^-OtmPqUmqbkqfijrn H + Pq w mq)b kq b jm 

- A^Qlcr m v mq b kq b jm + + ^r) /a mK + PkWk)WmkIntj 


jm 


k= 1 


+ Y^^~W + ^^-) 0a »" Q i( V i9 + PqWjq)WjIufoq 


3 = 1 
J /\ 




j=l 

r . An , 2Ae6 


A'n (4,6) = (-^~ + — J^)0top0n(v n + ^ntWn)Wptffcn^jp + ^12©^p W pn^Arn^p 

it it 


- -ap^lXpn^fen^ip 4 ^2“/5n(^pn + finWpn)^ kn^jp 

-A e 6 Qlall 3 n v pn b kn S JP + + ^)l<* 2 p Pn(v k + /3 fc ti>/> pfc / 3 o*j P 


*=1 


+ ]T(^. + ^l)Qa p a j f3 n (v jn + PnWjJwjlnSkn 


i=i 

J /< 


+ ^ jP(^ + ^¥L) aj a p Pn(Vj k + P k Wj k )w jk hoh3 


j- 1 *:=! 


I/ i /\ 

A'n (4,7) __ Ae 6 e f a m/3 n (v n + (3 n w n )6 mp 6 kn - -^-Ot m p n {Vj n + PnV>inWkn 


R 

J 


j=l 


j=l 


^22 


A'n (4,8) =~A660/« 2 m ( U n + /?nWn)<5mp<5fcn + ^ ” («n + ^n^n^mp^n 


- ]T ^~a v p n u ]n hb kn + + ^) 0 <W V n + 0nW n ) Wj h6 kn 


r ■ 

j=l j=l 


- A & eOa p atjVj n I 5 6 kn + ^(— p 4 ^-) 0a p a j( v j" + fi n ™i n)™jhif>k 

1=1 j=i 


j=i 

j 


+ (Vjn + PnWjn)hhn ~ ^ A66©tt p am(*>jn + PnWjn)hbkn 

3 = 1 " j = 1 
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A'n (4,8) 
conVd. 


K n (4,9) 


J K 


+ EE< + )ajOt p (v jk 4- PkWj k )wjkhoh4 + E AnOaj 


j - 1 fc = i 


A A 


j=i 


J I< 


*a p f3 n w jn Is6 kn + y E(^- + — + PkWk)wjkhho 


j = l fc=i 


= j4l2©/0! m (/?n^n 4" W n )6 m pfo kn 4" fin (®n 4" fin^n)^mp^kn 

+ y p 2 " fin( Vj n 4 ” fin^jn )^1 4 * ^ ^ - 4 . 12 H" ^jn)-^ 6 ^/cn 


i=l 

J 


i=i 

J 


+ ^2 3AnQRa 2 j& P a m WjW jn 


■^23^kn “h ^ ^ A\\ Q R(X. j&pOL m U j n I§f) k n 


j= 1 
J 


j = l 


I< 


+ ^ AuQajO!pWj n I 5 6kn + y + -y)«rai(^ + fikW>k)~ &mphl 


j=l 

J 


k~ 1 


^ ^ ^66 Oc^pQl jfin Vj ji I§ ^ kn “f - ^ ^ -^4l2 ^jn-^2 ^kn y^ ^ R ^ 


j=i j=i i=i 

*/^n 4 " y* n~ n )®®p®jftn{ Vn H - fin^n 


R 


i=i 


H“ ^ jp 4* jj )®®p&jfin{Vjn “1“ fin ^jn )^j^3 J i^kn 


24 


j = 1 

+ EE ^ ( 4f + + fik^ik?h2k 

j - 1 *=1 Z 
J K * ^ . 

+ y y (-F- + + fikw jk )hhi 

j=i k= i 

+ ^2 + —jp-)<*j<*pfin(Vk + fik*>k)v>jkhlM 


j - 1 /c=l 
J K 


+ ]T ]T(^ + + ^^^i fc/3o/34 

j=i fc=i 
J I< 

+ y y 1.5^ni?a^a p a m «;] fc /23/2S 

j=i fc=i 
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An (5,5) 


An (5,6) 


An (5,7) 


An (5,8) 


(^ + 

v f + »,)/, + £ 

fc=l K ‘ k = 1 

h v + &«>t)/7 + E ^r'K + /W/. 

fc=l K k=l 

y E + e e Hr'^ + 

j-l z i=l fc=l 

1 EE W* nrrtW* + thf+ 

j = i i-=i z j =1 * =1 

( Ai + &^)PWn, + ^0/VtMn, + 

t E + «-»)/, + E iTrA-M”* + 

fc=l " k=l 

f ]T -^T l ( V k + 0kWk)h + -JZT@ nl ( Vk + ^Wk?h 

km 1 K ~ k=l 

+ E 5(nr + + EE + 

j-i ^ K n i=i fc=i 

+ S ? ( ^f + + SS (_ jT + )a 2 j f3 n lw j w jk I 4 

j=i t=i ^ J =1 fc=1 

_ A660 Ct m i5 n W m f>nqf>jm ~ ^ —^~lot m (i q Wmkh^jm 
A fc=l 


-Otm^iymk + 0kWmk)hf>jm 

A K A i2 

ggO!^j)0/ui m Sjm^nq + ~^Qlo: m U m 6j m 6 nq + ^ -j^lct m U m kh*> jm 

U k= 1 


/3 . 1 /3. 
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I! 


A' ii (5,8) 

coni' d. 


A’n (5,9) 


+ E -firKPkVmk + WmkWjm + E + ^T^) a i 0U, i /l2<5 


Jt=l 

A 




i=i 


A’ 3A22 A 

+ ^ ~K V k “I" Pk^k){ v mk Pk^mk )^9^jm ~ ^ ^ jm 


k= 1 

J /\ 


/c=l 


+ £ £ W + '^y^oin + EBt + 


j=l /t — 1 

J A 


2 V i? 


j=l fc=l 


+ 55 E + PkWjk?hh\ 


2A 3 

j = l k — 1 

* 4,, A a 22 

+ / ] —n2~K v mk + Pk W m k )h^jm + “^/A^n ( tW + Pk^mk)!^ jm 

k — 1 k — 1 


A 


+ E ^MkVmk + U> m *)/4^m + E + ^)a^n0^/l2^n, 


fc= 1 
A' 


3 = 1 


A 3^ A 

4“ ^ ^ ^ ^ ^ Pk^mk)^9^jm ~~ 


k=l 

A 


k—1 


A 


' jm 


An 2 A 6 6 

/? + £ 

E ^7r<XmPklUmkl7Sjm + E^^ + “4p) /a mK + PkW k )w mk I 16 6 


+ E Ai 2 la 2 m l3 q w mk Is6 jm + E(— + —pT~) la2 m( v k + P kW k )w m I 7 b 

k - 1 

A 


k= 1 

A 


A?=l 

A 


Jt=l 
J A 


+ ^ ^/a m /? n U m */4^m + ^ + /W^ll 


fc=l 
J A 


j=l k=l 


v-E EE 1 A\2 2/166 x 2/0 2 r r ■ V^V^/^ 12 . 2^66 v 

+ E E 2^~R~ + _ R-} a jP nW j kIl ° l12 + EE^ - ^ - + R ) a m(*j 

j = 1 fc=l z J=1 fc=l 


J A 


x — ^ \ ^ A-i'y 2^4.66 x 9 n T T 

*{Vjk + PkWjk)WjI 7 In + 2^ E(-R- + -£-)<* jPnWjWjkUI 12 


j=l /c=l 


J A 


+ E E(^ + -~j^} ama j( V j k + PkWjk)WjkIuh 6 


j= 1 fc=l 


A A 
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An (6,6) 


An (6,7) 


A'n (6,8) 


= 2(^ + + ^QPlqoKg 

+ ^2 ^-^nWkVk + W k )h + i2 ^jrPnKVk + PkWk)h 


k= 1 
K 


+ E + E ^iA.W»* + A»* ) 2 A 

1 fc=l 


fc=l 

J 


+ V i(- Ai + ^r )a W mw 1 s " + EE Hr /(Vjfc + PkWjk) 2 PnPqh 

_Z K ri ? _i t-i 


j= 1 
J K 


J=1 /:=! 


+EE5(^ + H i) “' ?w '" ? * /, ° 


j = l fc=l 
J K 


2 V R R 

A\2 2 Agg 


+ E Dt + 

i=l fc=l 
^66©^ 


i? 


/< ^ 

®mfiq Wmfinqfijm “ ^ -^/a m /3 n /3g W m A;/4^jm 

/c=l 


+ ^ ^ — ^4“ & m fin l{Vmk 4 f^k^mk )^3^jm 
fc=l 

— A66«m)®^^ m ^' m ^ ni ? “I —QlamPnUmfijTnfinq 

R “ ft 

K A A 2 <> 

+ J] -^//3 9 (/3^mfc + Wmfc)/4«jm + 

*=1 Rm *=* 

/\ . ft 

+ T -J^f3nl3ql( v mk + (^k^mk)h^jm “ ^ A^a 2 m j5 q lw m kh^ jm 

fc= 1 ^ *=1 


fc=l 

K 


+ E + E i(^ + 

/c=i i =1 


/e=l 

K 


R 
3 /1 2 2 


3 4 

4- ^ -j2^-lflq(Vk + 0kWk){v m k + PkWmk)hbjm 
k— 1 

+ EE^ + H i) °^“’5^o/. ! 

j=l /c=l 


+ EE<^ + H^ ftwW " 

j = l fc=l 

+ ^ ^ + /?fc™jk) 2 A)/ll 

J=1 fc=l ^ 
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A' ii (6,9) =(^T0n + M-iot^QlwmbjmKq + ~-Qla m 0lu m S jm 6 nq 

K A ^ ty 

+ ^2 ~E5~iPq( V mk + PkW m k)h^jm + ^ -jglPn(Vmk + fi k Wmk )h^ji 


^fiPnflqiPkVmk + W m k)hf>jm + Y 2^~]f + ~ 


K 

+ Y, ^T l 0n0 q (Vk + PkW k ){Vmk + PkWmk)h^m ~ £ A 66 a 2 m Pqlv mk I 7 6 5 
fc=l " *=1 

+ ^ “* J^y a mPq( v k + PkWk)w m h&j m 

k = 1 fc=l 

+ —^lOi m fi n fiqU mk Ii&j m — ^2 — Jjr a mPqPkl u mkl7ftjm 


+ + —jr~) la mP<l( V k + PkWk)w m kh^3 


R R 


+ EE 


j=l *=1 
J /v' 


+ E E ^ 


j= 1 *=1 




j=l fc=l 


A A 


+ EE( ^ + + PkWjk)Wjhhi 


j = l k= 1 


4- ]T Y(~ir + ‘~^’^ lama ^^ v i k + ^wjk)wj k lnhf, 


j= 1 ^ = 1 


A"n (7,7) = (Ana^ + -^2‘/3n)^©A<5 mp (5 n9 

JC ^12 

An (7,8) =(Aj2 “l" ^66 )^ro/^n^® ^mp^ng 4* E -^-la m (v k + fi k Wk)h&r 


-E^“A0^„ + EE^ P (ujfc + P k wjk)hh 


j = i *=i 


_ EET aj ^ fc/4/5 


j=i ^=1 
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A A ^ A(y(y 

A' n (7,9) =A l2 a m Ql6 mp S nq + 53 -jf- la m (3 n {v k + f3 k w k )hS mp ~ 53 -jf la ™M v i' 


+f 3 k w k )lT 6 mp + ^2 A\\QRa m a v ajWjl 2 f>nq ^ ^ aj( 3 n 0 wjls, 6 r 
i=l i =1 


“ S S ^-(XjPnPqWjkhh -lb!! ^- a mP q {vjk + /W)/e/T 

j=l fc=l J =1 fc=1 

J K J A ^ 

+ ^ 5 ^ AuRa m a p ajWj k l 2 h + 53 5 Z —j^ a P Pn{ v jk + PkWjk)hh 
j= l t=i J - 1 fc=1 

+ (A* + ^-)alieR 6 mr K, + ^ 6 («> 0 «, 

AT . A' , 


i=i fc=i 


A'n (8,8) =(^ + ^)0le I6 mp 6 nq + (^66 + ^)allQRS mp S nq + -^Qlw 0 6 mp 6 t 


/\ A 

+ ^^QqotimpKq + 53 ~ET^ kVk + W k)l4$m P + 53 ~# 2 ~^Q^ Vk + PkW k ) 
2 K , „ a-=i 


*/ 3 <5mp + ^2 15rP nl ( Vk + PkWk)l7bmp + ^ + U?/, )“ / 9 <5 r 


+ ^ A 660 a j( a m ^2 + a p ^5 

j=l i =1 

+ t + t blf + T? 1 )"?®"?'-*" 

J = 1 J = 1 

+ 53 E + 53 S Hf^ fc + P kw i k ? hh9 


j= i fc=i 
J A' 


j=l A:=l 


+ EE EE 'fPq{vjk + fhwjk)hh 


j = l k~ 1 
J K 


j=l k - 1 


+ 53^3 jjrPn(v k +PkWk)hh - 5353 /l66 ^ (am/2 + 


j=i jt=i 


j=i *=i 


+ E E ;<tt + 2 J ^)°W,^ho + EE<ir + )«’ 

j = 1 fc=l z -*- 1 fc - 1 


+ EE^ (Ufc + P k Wk)(vj k + P k Wjk)hh 

j= i fc=i 
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K n (8,9) 


A 22 , o2 D 22 \ , ,D\2 + 2-0(56 , 2 


[(^ + ^^r) + ( 


^22 


R ' iZ 3 7 ' v R 


)al n \P n Ql6 rnp 6 nq + —^Qpn^ w 0 ^mp 


’nq 


+ -^-^-Q^nQo^rnp^nq + A\2 jOtm&n w j^^nq 4* ^ ®fi n Wj I\ S nq 


2R 

J 


i = 1 


i=i 


+ E - E ^600 .a.DnWjhK, 


3 - 1 
J 


j=l 


A' 


+ ^ -j^-Qotjf3 n v,jIi6 nq — ^2 


i=i 

/\ 


*=i 


+ ^33p n 8 q l(Vk + 8kWk)h&mp + X] -Jfi-WniPkVk + Wk)U&r 

fc=l k=l 

+ ^2 + 8kWk)hbmp + ^2 ^W'^ Vk + ^ fcWfc ) 2 ^ 9 ^ m P 

/c=l ^ fc=l 

+ EE ~jD~® j$ n ^ jk^l ^4 4 EE + PkWjkfhhd 

j = 1 *=1 ^ i=l fc=l 

+ E E #A.(A”i‘ + <■>*)'.* + E E TfrM-te + Pku>ik)hh 


J mp 


j = 1 fc=l 
J A 


j=l fc=l 
J A 


+ S E 4^-( u jfc + 8kWjk)hh ~ ^2 S A 66am(«i^2 + a p I e )v jk I 7 


j=l k=l 
J A 


EE 


j = l fc=l 
J A 


j = l &=1 

J A 


£ 


j = l /c=l 
J A 


E E M^ja r /3,w it i,h + e E + 

j - 1 fc=i i= 1 *=i 

-EE - 4 e + £Et + ^ )<*?&> >w^» 

j = l fc-1 j=l * = 1 


J A 


+ EE^r /?n(»Jfc + 8kWk)(Vjk + 8kWjk)I\h 


j = 1 fc=l 
J K 


-R 3 


+ E E<^r + + fo w k)w 3 hh 


i=i fc=i 
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K u (8,9) 

coni' d . 


A'n (9,9) 


J K 


+ Y2(~jr + -jr-) Q j a m(vjk + PkWjk)wjhh\ 


j - 1 k=l 
J K 


+ Yj + -^) Q j a m( u * + PkWk)wjkhIi6 


j=l /c — 1 
J A 


+ ^ + flk w jk) w jkhsh\ 

j = i fe=i 

:[(^jy + Pn~jgr) + 2 ( 1- —) a m0l + A>11 -Ra^] 0M mp ^n 9 

+ (^12C*m + ^p”/?n)9^0^mp^ng + (-y-# a m + -^/?n)®90«5mp^n</ 
j J ^ 

-|- ^ ^ A\ 2 Q[oi jOt m l2 + CXjCXpIs + OtpCt m I§ )'WjS n q A ^ ^ ^ ^ 


nq 


j = 1 

J 


J = 1 


+ H ^)a]0/?^U?j/ 2 O^ng + Ai i 0 Ra j OCpCXm u j h & 

U 2 R R i=i 

j * 

+ ^ L5AneRa)a p a m wp 23 S nq + ^ Ai 2 /a^(/3^i;fc + w k )h^mp 


nq 


i=i 

A 


fc=l 


A 


+ Y1 jjrWiVk + 0kWk)h Sm P + Y2 + W k )I 4 6 


mp 


k= 1 
A 


fc= 1 
A 


+ Y2 ~^r^P<i( Vk + PkWkjllbmp + Y2 21 p>(W Vk fi kWk ^ 


k= 1 

K 


3A22 


fc=l 


+ V \{-JT + -jp-) a m l ( v k + PkWkf hs&mp + ^ ~^-Q a jl^l. u jh^nq 
t= 1 J- 1 


fc=l 
7 A' 


+ E E + E E + /W) 2 -w. 9 

j=i fc=i i=i fc=i 


J I< 


J K 


+ Yj -pfflnPqiPkVjk + W jk )hh + EE ^-/3n(Vjk + /3kW jk )hh 

3 = 1 k=l K ‘ j=lfc=l 

+ EE ^-/) 9 (ujfc + fi k Wjk)I\h - EE 

i? ^ ' 


j=l /c — 1 
J A 


j = l /c=l 

“ - 4 j * 

- ]T J! Y1 A M a m<XtPqVjkhh 

j= 1 fc=i i =1 fc=1 

J A 


- ££ A maj a p 0 n v jk hh + EE + ^)a?A./»X*JWjo 

j= 1 ^=1 


1/^12 , 2^66 


j — 1 A:=l 
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A'i i (9,9) 
cont'd. 


J I< 


- + + —jr ] 13 n ^ Wj w j k /4 /2 ° 


j= 1 k = 1 
J K , 


3^2° 

+ 2^ -nTPnfiqiVk + PkWkHvjk + PkWjk)I\h 


j=l k=l 
J K 


v -^ { A\2 2/156. 

+ 2^ ZJ~/T + ~^) a i a 'mPg(^ + Pkmfajhl 7 


i=l /:=! 
J 7\T 


y""v ^12 2.466 

+ 2^ 2-J^T + —p-) a j a m/3q(Vjk + PkWjkjWjhhl 


j= 1 k=l 
J K 


Y^, -4 12 2^66 


j=l k = 1 
J 7\ 


Y~'V _ '/^12 2^66. 

+ Z. 2/ “ft" + —ft-) a j a mPq(v jk + /3 k W jk ) Wj I 16 I 2 i 
j= 1 fc=l 

• ; ^ J I< 

^ ' .4] j /(’ajQpQ m Uj k I 8 I 8 'y ^ ^ " \.bA\\Ra‘j(Xpa m w i jf.l22l2 

j=l k=l J = lfc=l 

j I< J K 

+ YL YL SAuRajapCimWjWjklsIoa + Z 5Z A i2a p a m f3 k v jk I 6 I s 


j=l k = 1 


j - 1 *=1 
J A' 

+ Z Z ^i2(oy«m/2 + Oja p I 5 + a p a m I 6 )w jk I 8 

j= 1 A= 1 

Y- Y 1 4j2 2/166 o 

^ z_ ^ 2^ 2 ~ W Jj~' a ™ a p(Vjk + PkWjkY I 241 28 

j= 1 fc=l 
J K 

Y^ y^/ ^l 9 z/i66 

+ Z ( “^" + — ^)«p«m(W/t + PkW k )(v jk + 0 k w jk )l 6 l 28 

j=l Ar=l 

J K 

Y^ y^ /*i 9 z/i66 

+ 2^ Z^ (- jr + + PkV>k)V)jhh 


.7=1 Ar=l 

J 7\ 


Y^ y^/ 4j9 2/166 

+ 2^ (_ ft~ + — )«j« P y97i(^A: + 0kW jk )WjI 3 I M 


i=l Ar=l 
J 7\ 


Y^ y^/ -4i° 2/166 

+ Z-, 2J "^r + -^)aja p /? n (v* + Pk™k )™jkh he 


j = 1 A:=l 
J 7C 


V“^Y^ /^1 2 24.66 x n 

+ Z^ ZZ "IT + "^"J a i a pPn(Vjfc + ^fcWjit)Wjfc/26/34 


j=l fc=l 
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Elements of {3UN + 2M + 2N + 2) x (3MN + 2M + 2JV + 2) submatrix W«.wi)l 


1 ( 1 , 1) =0 

1(1,2) =2 Q{R + wo) 

L(l, 3) =0 
L(l,4) = &W m fijm 
1(1,5) =(»„ + /3n^ n )0^fcn 
Z/(l,6) ={fl n V n + W n )0Skn 

1(1,7) =0 

1(1,8) =^(t> mn + /*n ^mn 

1(1,9) =i(/?nUmn + W mn )®bjm^kn 

1(2,2) =20(2/ + go) 

1(2,3) —2a m Qlw m fijm 
1(2,4) =2 a m Qlu m S jm 
1(2,5) =0 
L(2,6) =0 

L(2, 7) —Q. m Qlw mn b jm^kn 

1 ( 2 , 8 ) =0 

7,(2, 9) =a m <dlu mn &jmhn 

L( 3,3) =0 

j 

L{ 3,4) =2 a m Ql(R + wo)6 mp + ^2a p 0Wj/i 

i=i 

7 ,( 3 , 5 ) =a p Ql(v p n + 0n w pn)f>jptikn 
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£(3,6) 

£(3,7) 

£(3,8) 

£(3,9) 

£(4,4) 

£(4,5) 

£(4,6) 

£(4,7) 

£(4,8) 

£(4,9) 

£(5,5) 

£(5,6) 

£(5,7) 

£(5,8) 


--a p Ql{Hn ^pn + Wpn kn 
= 0 


-a, 


u 2 

i&l(v n + fln w n)hmpf } kn + “<*p/3n © ^kn 


J = 1 


j J i 


J = 1 


j=i 


=a, 


_i 2 

7 0/(/? n V n + Wn)<5mp^&n + — Q' p/3 n 0 V j n ^1 <5 fc n 


i=i 


J J 1 

+ ^ ^ ^pQ^jnA^/cn + ~^&n®{° L jIh “I" Ck m^6) v jn^kn 

j = 1 J=1 

J 

= 0(2/ + 9o)^mp + 2aj0?Xj/i 

j=l 

0/lXpn^ jp^fcn 

"Clp0/ XXpn ^jp^/cn 

^ ^ 2 

= Q? m 0/(/?n^n “f" )^mp^n “1" ^ ^ ^Q(p/3 n 0^jri-^6^A:n 

i=i 

J J 1 

+ a m 0Wj n /l<5fc n + -^n0( o m^l ~ a .?^2 ) v jn^kn 


j = 1 


j=l 

J 


2 1 

-a p /? n 0(ij n /5^/cn + -/? n 0(aj/i — Oirn^UjnSkr 

J=1 ^ j=l 


i-1 3 

J 

' Otj 0 Ujn^l $kn 
j=l 

= 0(2/ + tfoMng 
= /?„©( 2/ + ^o)^ng 


* 1 

tt m/^n0^m^ng “I" ^ ^ “1“ )^7^jm 


fc=l 


A' 2 ^ 1 

+ -r&mPqlWmkh m + -OmPn^mkh^ jm 

tl 3 k= 1 3 


K 


-Q-m 0 lu m ^nq + y ^ 


*=1 
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1(5,9) 

1 ( 6 , 6 ) 

£( 6 , 7 ) 

£(6,8) 

£(6,9) 

£(7,7) 

£(7,8) 


£(7,9) 


£(8,8) 


K 2 * j 

- ttrM 0n 0/^tti &nq "i” 77 ^Tn/(A^7z£l "f* 


J TYl 


k = 1 


/c=l 


= 0(2/ + 9o )£ 


n 9 


K 1 

= O m 0/i£- m ^ ^ — Z?n^3 

*=1 ^ 

A' 2 * 

+ Es“ m^qlVmkhSjm + a ml w mkh&j 


jm 

A 


2 1 
= ^T7i/^7i 4“ ^ X^77i/^g^^77i/c^4^jim “I” ^ ^ 

k= 1 15 Ar= 1 * 


/v 


“h ^ ^ ^ m^^m k ^8 &jm 


k=l 


K 1 K 2 K 

—® m l{flnl 8 ~ ^A^mp H~ ^ Oi m lv^I^6 m p 

k= 1 ^ fc=l ^ *=1 


J A 


2 1 

+ ^ -OijPnQwjhbnq + ^ o^n®( a P^l _ Otmh)WjS nq + ^ ° p t; j A; / 1 /3 


J = l 
J A 


i= i 


j = l A: = l 


J A J A 1 

+ EE ^/^(^pA H~ ®ml(> )^j "f* EE ^/^n^pA 4~ & jl$)m jkl 8 


j = l A = 1 
J A* 


j=l Ar= 1 


+ EE 


j = l k— 1 


A 


= a m 0/(tf + ^o)^mp<5n ? + ^a p 0^j/i(5 n9 + ^ - a m lv k I 8 6 


j = 1 *=1 

A 


2 7U * ~ K ~ o ~ trip 
J K 


1< 1 K J l< 

+ £ + /WM,., + E + E E CXpWjkl 1 A 

k= 1 A:=l j = l A: = l 


Ar=l 
J A 


J A i J A ^ 

+ EE — /A(a p A + a mA) v jkA + ^ g/?n( tt pA tt jA) v jfcA 


i=i fc=i 
J A 


j=l /c=l 


+ X] O Pq( a mh - Otjh )Vjkh 


j~ 1 fc=l 


J A 


j j j\ 

= — 0(2 / + ^O^mp^ng + a j© w jA^ng + ^ a j u jkAA 

1 j=l j=l A:=l 
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Z(8, 9) — ( ^(3 n Q{ { ll + Qo)^mp^nq “I" ^ ^ ^ Otjf3 n QUjI\b n q + ^ ^ 

J-l j=l 

J K i J /\ 

“I" ^ ^ ^ Qplb )^j/c^7 "I" ^ ^ ^ ^ )^ j Ar-^4 

j=l Ar=l j = l *=l 6 

J K 

+ EE 3~ QjIl)UjkIg 

j - 1 A: = l ** 


X(9,9) =| 


J J /v 

0 ( 2 / + qo)6mp6nq + ajQujI\6 nq + ajuj^hh 

j=i j=i fc=i 


Elements of (2 M + 1) X ( 2M + 1) submatrix [A' 22 (Wstr)] 


A 22 (!■>!) = 


(EA), 

21 


A '22 ( 1 , 2 ) =0 


A '22 (1,3) =0.5(EA) s a 2 m w sm 
A '22 (2,2 )=(EA) s a 2 m l6 mp 


j 

A 22 (2,3) — / ^ ( E A) gQj-^cx jOipW s jlq 
i=\ 


j 


A '22 (3,3) =[0.5 (EA) sqi + ( EI) s a 2 J } 


+ ^ 1 . 5 (£ M ) s a m Q ^ a p w ^/ 2 3 

j=l 


■^2 )^j^nq 
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Elements of (2 N + 1) X (2 N + 1) submatrix [/i33(«r)] 


A’33 (1,1) 

A' 33 (1,2) 
A'33 (1,3) 
A '33 (2,2) 


A 33 (2,3) 


Am (3,3) 


EA)r Q 

Rq 

— ®( v rn + fintUrn) 

v K , „ 


' Ri yn Rq 

' (~7? T ~ 0<}( Vrk + 0kWrk)h + X - ffT fin ( v rk + fik™r k)h 

\ R o k ~i o 

' *££ (■» + PkW-rk) 2 h + X ^ 7jT~(& kV rk + w rk)h 

2 R 0 k=l 

+ fil PnQS nq + ^~^r~ fin&^rO^nq 


■+^ { ^}f3nQS nq+ { -^(3 n Qw r oS. 

J' 

[^L/} n p q (v rk + ( 3 k w rk )h + XI + ^^) /7 

Ao *:=i 0 

^^p3~ Pn( v rk + fikWrkf h + X ^"d2 ~fin(fikVrk + V>rk)h 
2 A 0 jt=l 0 

- + fit } ^3 ~]Q^»q + ~ fit® w rof>nq 

( -^T L fin(Vrk + fikU>Tk)h + X + PkWrk)h 

R o k = i " 0 

/if ( EA)r n n < a „ . ... 
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li 


The integrals I\ to I34 in the elements of the submatrices K\j and L are given by 


i 

I\ = J Cos(ajx)C os(a m x)C os(a p x)dx 

-1 

1 

1 2 = J Sin(ajx)Sin(a m x)Cos(a p x)dx 

-l 

© 

1 3 = J Sin{Pk9)Sin(f3 n 9)Cos(f3 q 0)d9 

-© 

© 

1 4 = j Cos((3 k 8)Sin(j3 n 8)Sin(l3 q 9)d0 

-0 

1 

Is = J Cos(a m x)Sin(ajx)Sin(a p x)dx 

-1 

1 

-i 

© 

/? = J Cos(f3 n d)Svi(f3i c 6)Sin(f3 q 8)d9 

-0 

0 

/ 8 = J Cos{P n 0)Cos((3 k 0)Cos(/3 q 0)d9 

~© 

© 

h = j Sin 2 (/3 k 9)Sin(/3 n 9)S in(fi q 0)d0 

_© 

© 

/ 10 = j Cos 2 (P k 8)Sin(M)Sin{/3 q 0)d0 

-© 

1 

In = / Co5 2 (ajx)Co5(a m ar)da: 

-I 

l 

I 12 -I Sin 2 (ajx)Cos(a m x)dx 


195 



© 

/ 13 = J Sin(Pk 9 )Sin(/ 3 n 0 )Sin{( 3 q 0 )dd 

-© 

i 

iu = J Cos(ajx)Sin(ajx)Sin(a m x )dx 

-l 

© 

/ 15 = J Cos(M)Sin(PkO)Sin(( 3 g 6 )dd 

-0 

0 

/ 16 = j Cos{/ 3 k 0 )Cos(l 3 n 6 )Sin{l 3 k 6 )Sin{( 3 q 6 )d 9 

-© 

/ 

/ 17 = J Cos 2 (a jx)Cos(a p x)dx 

-i 

i 

hs = I Sin 2 (ajx)Cos(a p x)dx 
i 

7^9 — C 0 S“ {QLjX^C OS^QtjnX'jC OS^QipX'jdx 
l 

ho — / 5 i ?i 2 ( q j a* ) C o 5 ( a m x ) C o s ( a p x ) d x 

-i 

/ 

hi = / S in(a jx)C os(a jx)S in(a m x)C os(a p x )dx 

-/ 

0 

I 22 = j Sin 2 (p k 9 )Cos{ 0 n 0 )dO 

-0 

l 

hs — / 5in 2 (ajx)5trc(cK m x )St7i(a p 2)dx 
/ 

I ‘. 24 = / Co5 2 (ajX)5m(a m x)5in(a p x)cfx 

-i 

© 

/ 25 = J Cos 2 (Pk 0 )Cos(/ 3 n 0 )d 0 

-© 
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!i 


© 

h§ = j Cos(j 3 i c O)Sin(f 3 k@)Sin(l 3 n 6 )Cos(l 3 qd)dO 

-© 

© 

J27 = J Cos 2 (p k 0)Cos(l3 n 0)Cos(l3 q 6)de 

-© 

© 

/28 = J Sin 2 (j3 k 0)Cos(j3 n 0)Cos(j3 q 0)d0 

-© 

/ 

/29 = J Sin(ajx)Sin(a m x)Sin(a p x)dx 

-i 

© 

ho = J Co$(fi k 0)Sin(l3j c 0)Sin(f3 n 0)d0 

-© 

© 

hi = J Cos 2 (f3 k 0)Cos{/3 q 0)d0 

-© 

© 

/ 32 = J Sin 2 ((3 k 0)Cos(l3 g 0)d0 

-© 

/ 

I„=JcMc S x)SM^)SMo r ^ 

-l 

l 

h 4 — J Cos(ajx)Sin(ajx)Cos(a m x)Sin(a p x)dx 

-i 
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APPENDIX D 


ELEMENTS OF RESIDUAL FORCE VECTORS 
FOR NONLINEAR ANALYSIS 


The non-zero elements of the residual force subvectors R s heii, Rstr and Rr in Eq. (4.81) 
for nonlinear analysis are listed below. The parameter S tj is Kronecker delta assuming the 
values zero for i ^ j, and one for i = j, respectively. 

Elements of (3 MN + 2 M +2iV + 2) X 1 residual force subvector [R s heti( u s h e ii\ p)] 

For the sake of simplicity the residual force subvector [ Rsheii ] is written as 

n ^ next. pint . 

ftskell — P * shell “ shell 


in which the elements of are 


pint. / i \ _ 

* shell V 1 / — i 


+ 


M 


/luOR 9 •> , , . a 

E E -1— 1 “»«*•» + E E -j r (v ” n + 

m=ln=l m = ln=l 


N , M 

F%i, ( 2 ) = 2 / 1,2090 + + E + E A..ei« 


2 2 


m^m 


+ 


n=l " rn=l 

/ , 9 ^^A 22 0/ * )2 

55 — 2 — amW; ” n + 2^ 2R~ t7mn + ^ mn 

m = 1 n=l m=l n=l 


771 — 1 ft — A 

Ai 2 a m ©/. 

R + 

n=l 


* (t>n + fin w n)bmp + 


+ 


m = 1 

M N . 

^(^n+^ mn )2 /l 

m = l n=l 


M iV j 

55 ^11 QR a rn a P w mIl 4" 55 55 ” 2 

m= 1 m=ln=l 


^ A A u 0R 2 

55 55 
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^she'll ( 4 ) — 4 “ 2 ( 4 ” -^11 ) 0 /-fttiJ m <$ m p “I" 2 i 4 j 2 ^rn, 0 ^M^O^m 

4 " >li 1 0 ft< 7 o ^m^mp “I" ^ ^ g ( ^>~ 4 ~ _ ft )0^J^(l> n 4 ” fin^n )** ^m^ 7 ?i 


AftaOl 


Qmfin(Vn “I" ^n^nj^mn^mp “I" ^ ^ A 12 0 /ctJ^ (/?n^n 4 ” ^n)^mn^r 


+ E- ^q" ( V 71 + /^n^n)(^mn 4“ fin ^mn )$mp 4" ^ ^ 2Aj 1 0/^0.'^ CXpU m t^m /4 


~ X! ^6® /a m( v n +^n%Kn<5mp + ^ i4 U J13 + ^ 4l2<*m 

n=l m = 1 m = 1 

M M N 

* © ^ra 1 2 4 ~ ^ ^ 2 , 4 i 2 0 Ckm^p ^ m ^4 “ ^ ^ ^ " < 4 - 66 0 ^p°m ( ^ron 4 “ / 3 n W m n )^mn^ 


m = l n= 1 


JVJ /V ^ 2 ^ 

"H ^ ^ ( j^~ 4 “ ~j£ ) 0 &'m< 2 p(^n 4 " / 3 „ W n )w m V mn l 4 


I V' V" l ( ^ 12 I 2 ^ 66 m / * \2 r , ^22© 

+ 2 ^ 2 ^ 2 ”ft~ + -^) 0 a mM w mn +/ 3 n^mn) tO m /l 5 + Z^ Z^ ^ 2 ~ 


m = l n=rl 


*(^’mn 4 “ fin'WmnY I\ 4 “ ^ ^ ^ ^ 4 “ ^ ^ )®&mQtpfin(Vn 4 ~ ^n^n)%^mn ^4 


+ EE /ll20a m ap(/?nVmn + Wmn)w m nh + E E ~^~ Q m u; mn^ 


+ E^ EZ l-54ll0/tamapW m lU^ n Ji3 + E^ El ^110^ol«p“mn»mn/4 


-EE- _ " <^p/ 3 72(^77171 4 " fin^mn) UmnI 4 


^22 , D 22 


A 22 .-,') D‘22 


FI hell ( 5 ) = 2 (-^ + -^-)p^Qlv n 6 nq + 2 f -gi 4 - ft j)/Jn 0 (%^, + -^©<Zo(t’n + 

2^22 ^12 
* W n )finq ~i" ^2 ®^o(^n “I" Pn^n^nq ~\~ / , ^ “I" 


+ Z |(“ir + T^ 6 ^" +/ J nWn)roi^r 1 , + E^ ^ 12 a 2 m P n Qlw m w mn 6 ri 


2 y R R 


M A 

^ QlQm&n^m ^mn 


& nq ^ ' ■ | 4 ( 560 /<l jn ti> rn t) mn < 5 r ig -(- E ' pV Olw ri 
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Khell^) _ + (. Umn + I3 n w mn )6 nq + -SrPA v n + ^nWnfh + E + It 

cont'd. “I -« n=i 

JV M N A l 

+ E ~lj£~ ( V n + /?nW„)(/? n V n + W„)/s + E E + ^ nU?mn ^ 

i ^ m — 1 n=l 


M jV ^ ^ ^4 12 2^66 

-f- ^ ^ ^ ^ - ^~ - lQ m (v mn + (dn^mn) u mn AQ6{X m lv mn W mn + ( ^ )l 


m = 1 n=l 


*a^(v n + /3 n W n )w m W mn ~~~Ol rn f3 n lu rnn W mn + -^-{v mn + PnWmnKPr 


1 E 1.5A 2 2^ , „ w , „ \2r , V' V' 

+ W mn ) h + X/ E R3 ^ nWn ^ Vmn (^n w mn) h + E E 2' # 

m = 1 n.= 1 


m=l n=l 


M N 


m= 1 n=l 


M N 


1V1 i ▼ 

+ ^)a 2 m /(t> n + /?„u; n )«i n J8 + E E 0 - 5 ^ 2 ^^ 

^ m = ln=l 

Fjjft, (6) = 2{^f + Pl^-WnQlvJ nq + 2(^1 + ti^ji-)QhOnt> nq + ^0/3n</o(/3n^n 

M ^ 

+ U n ) 6 n? + ^QW nW 0 (v n + PnW n )tn q ~ E 

m=l 

Af , 

+ — a m /3n0i(^mn + Pn^mn )^mKq “ ^66©^ a m^n^m u mn^nq 

" ™=i 

M ^ M 

+ E “(“=■ + + PnWn)Wmf>n q + E ^ w m^mn^n q 

^l 2 R R 

W /t _ A 22 I 3 v ' ^22 ^ 

+ E + PnWmn)bn q + E ~^ 3 ~^^ Vn + ^ n ^ n ) + E #2 

771 = 1 K " = 1 n = 1 

N M N . 

*(v n + PnWnfl 3 + E + @ n w n)(Pn v n + W n )F + E E [ # Z 

^Tj " m = l n=l 

1 , A\2 2^66 \ 2 1 ( 

*a m (»m» + finWmn)Umn - AeeO^^nWmn + (“g" + fl - ) a m lw mWmn\Vn 

+ f3 n W n ) - ^-a m p n lu mn w mn + -^pr(Umn + PnWmn)(PnVmn + Wmn) 


M N 


+ E E 1,5 pa ~&( l?n + ^W n ){v mn + PnWynnfh + E X 

il m — 1 n — 1 


m= 1 n=l 


m=l n=l 


Af TV 


+ ^)Q 2 m /?,/(t; Tl + /J„«>n)w ^„/8 + E E 0 - 5Alja "' W ™ /6 

^ m=l n=l 
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Fshi'il (') - + Pn~^r)®lR u mn + (A J2 + A 6 6 )®m(3nQlv mn + Ai 2 0: m Qlw ri 


N 


Ot m f3 n Ql(v n + &n'W n )'Wm ^mp^nq 4 ^ ^ p 4* PnW n )(v mn 4* Pn 


n=l 


*^mn)^mp^3 ^ ] ~~p~ ®mPql{Vn 4 Pn^n)^mn^mp^ 5 — ^ ^ R /^ti 0 


n= 1 


m=l 


A/ 


A/ TV 


*(^mn 4 Pn w mn )^m ^ng^4 4* ^ ^ 4" ^ ^ ^ 

m=l m = l n=l 

A/ AT ^ 

* a p(^mn + /3n^mn)“A^3 “ ^ ~j^&mPq{ v mn + Pn w mn) w mnhh 

m = l n=l 

A/ AT 


A 12 

2R 


m = 1 n = l 


F shell (8) — |(^i2 + Aee)a m /3 n Qlu mn + (— -jp + -^~)/?n + (A 66 + —^-)cr m R 


, \ M2 q2 D 22 
+ 1 R +i3n R3 +a 

^22 


R R* 

2 (^12 + ^D^e) 


R 2 


Qlv rj 


R 

, 0.5^12 


/3 n 0/l^ mn 4" ( — — Q m U m — .4 6 6 Qt 


R 


m W m 


4 2 Wm )( v n 4 ftnWn)Ql 4 ( — <70 4 ■^ 2 '^o)(%n 4 P n w mn )© j 

TV 

+ ^2 {^~^ Wmn ~ -^66 OtmW mn )la m (v n + & n Wn) + —jgf- [(1 + 0l){w n v 

72= 1 

,1 " Ao 2 

4 Vn'^mn) 4 2/? n (l? n V mn + Wn w mn) J j^mpP 5 4 / ~^T PqK v n 4 /? n tU n ) 




A/ 


* ( ^mn 4 Pn^mn)^mp-^3 4 ^ ^ 4 Pn^n) i^mn 4 Pn^mn)&mpl 7 


n — 1 


M A A M 

+ Yl^~^ amUm + -j?T W rn)(Vmn + PnWmn)&6 nq I\ + ^ [Ai 2 a^/3„ 


w. 


m = 1 


m = l 


. I/A 12 2A 6 6s. 2 / , Q s , 2 A 66 

4 2 V 4 ^ 4 /^n^n) ^ 66 ^ 771 ^^^ ^mPn^Ti 


0^ ? m ^2 


^ yl|o 2ylg6 ^ 

4 ^ ^ o ( ~cT 4 0 “ )®mQ( V rnn 4 Pn^mn) w m^nql 11 ~ ^ ^ 


m = l 


m=l 


Af A/ 


m = 1 72 = 1 


E r— ^ 4^2 ^2^ 

/ „ ®m“mn + (finVmn + W mn ) 


(v* 
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03 



Fikti i ( 8 ) 

conVd . 


0.5A 2 2 a , , a t . V' V s L5 ^ 22 

= +/3n^mn)/l^5 + XI X/ P2 ^( t7 ” ln + Perrin) h 3 + /_/ 2-^ #3 

^ m = l n=l 


m= 1 n=l 


M N 


*{v n + PnW n )(v m n + PnWmnfhh + YY Jj 3 "( t ’ ran + Pn™mn) hhi 

m = 1 n=l 

M N 

+ EE 


m = l n=l 


~ Ae60‘m v mn —a m 0 n Umn 

R 1 1 


M N — - 

*W mn hh + XI E°- 5Al2 ^^ /2/6 ~ X/ z2 A ^mapPn w mnhh 


M N 


m= 1 n=l 


771 = 1 n=l 


M N 


+ E E 5 <tt + 


m— 1 n=l 

M N 


2 y R R 


-f ^ Y/(~JT + -^■ ) a m( v mn + ^n w mn) w m w mnhhl 


m — 1 n=l 

M N 


+ E E + XT + 


77i= 1 n=l 


Fikiii (9)= {.-t,, + i^f + ""X? + 7; 

4h + + 2 (P.i + ^ ggl a j^ + + ( 


,0.5Ai2 


R 


-qo 


+ 4 lr^o)(Vmn + PnWmn)Pn + ( 0 . 5 Au<? 0 # + Ai 2 «)oO Q m w mn 
it* 

+ ( — QmUm + ^U> m )(Vn + (3nWn)Pnl + A \ 2<*m / (/ 3 n*>n + W n )w m }e6 mp 6 

R R 


- V A 6e a 2 m l(v „ + t3 n W n )v mn h6 mp - X/ + ^nW„)u mn /3«mp 

_“A n=l 


n= 1 
/V 


+ ^ |^-a m -li mn (Un + /3n^n) + + u n™mn) + 2/?n 

n= 1 

N , , 

*K*W + «; n tff mn )] } PqlSmpIs + Y ~W^ n + + ^mn^mph 


N N 

+ V l ' bA ?- {v n + PnWn) 2 (Vmn + 0n™mn)P<,6mph + Y A U a m lw ™n{PnV n 
“ R J n=l 

71=1 N 

+w n )K r h + E i(^ + ^K'htW’%1.,6 + E j(^x 

n=l 
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^convP = + + PnW n ) 2 w mn 6 mp Ii7 4- ^ 4^a m /3„0u m ( v mn + f3 n w mn )6 


m = 1 


m \ A 2A 

* h -f Yh [^ 1 2 a 2 m Wmn + 2 ^~W~ + —jr~) a m 0 n W m( v n + fin™ n) - A 6 G a 2 m ( 3 n V T 


m — 1 


^66 

R 


M 

Q^m^nql 2 4" ^ ^ j4l2 0Om®p(^n®mn 4" 2u!mn)®m^n5A 


m=l 


M . M 

E /\. 9 9 ■> o 9 

~nY /^n®%(^mn 4" /^n^mn)^n^l 4" / ^ 1.5i4n0i2o; m Q!pW7 m iyf7in^ng-^13 


m = l 
M 


m=l 


A/ 


+ E + ^) a mPn®{v mn + p n w mn )w 2 m 6 nq I u 4- A„0£a 2 m a p 


m = l 


rn = 1 


AT N 


iVi A A 

*(^m w mn 4* 4" ^ ^ ^ ^mn 4* ( fin ^mn “I” ^mfi) (^ri 

m=l n=l 


Af N 


M N 


1.5.4 4 


fin^mn )fiq^l ^5 4“ ^ ^ ^ ^ ^4l 2 ^p(/^n ^mn 4" ^mn)^mn ^6 4” ^ ^ ^ ^ ^>3 

m= 1 n=l m = l 77=1 

A/ N nc , 

\3 


0 5A 

*(^n 4" fin^n )(^mn 4“ fin ^ mn )~ /5g/l ^7 4" ^ ^ ^ ^ ^>3 * (^mn 4" /^n ^mn ) fiqll -^12 


777 = 1 n=l 


A/ JV 


\ — ^ ^ T/ - 4 i 2 2466 0 / o x . 9 466 a 

+ 2 ^ 2 ^ [(“cT + “]p) a m<>n + fin^n)Wm “ 4 6 6 « m *W ~ — a m/ / n U 


466 


m = l n= 1 


A/ N 


M N 


* fiqWmnfals 4" ^ ^ 2 y 1 ‘54u (X m <XpRw m W mn I^I\^ 4" ^ ^ ^ ^ 0.54 lo^m. ^mn^6 


m = l n= 1 


777= 1 n=l 


A/ TV 


+ EE- ^>o~ ( ^mn 4“ fin^mn) I\I^ 4“ EE° .5Aua^ n a p Rw^ nn Ii3l\s 

777 = 1 77= 1 777 = 1 71 = 1 

M N M N ^ 

4“ ^ 2 v 4u Oj^Qfpi?lZ mn 1L? mn /4/6 ” EE Ckp/^n(^m7l 4“ /^n ^777 71 777 77/3 


A/ N 


777=1 71=1 

M N 


777=1 71= 1 


M N 


4“ 2 ^ y ^ h ^ )^m ( ^mn 4“ fin'^mn)^rn^mnfiqlhl\\ ^ ^ ^ 4g6^m^j 


777 = 1 71= 1 


m = l 77=1 


A/ N 


*( ^777 77 4~ fin ^777 71 ) ^m 7 l/ 3/4 4“ ^ ^ ^ l(Ap + ^^-)a 2 m (v n 4- () n w n )w 2 mn P q hh 


777 = 1 71 = 1 


2 V R R 


M N 


4- 2 D^ir + ^r )ama P( Vn + ^^(Vmn + PnWmn)Wmhh 


777= 1 71=1 


i? R 
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Hheii ( 9 ) = + y y I(ii 1 + ^l)a m a p (v mn + l3 n w mn ) 2 Wmhhs 
Hr! ^ ^ 2 i? R 

m = l n = 1 


coni' d. 


.)_ y Y{'—— + -—jP~) a m a p( V n + Pn w n){ v mn + Pn w mn)Wrnnhh’ 


77i — 1 n=l 

M N 


Ai2 2A& 6 
~R + R 


i V' — 1 — — — - -f" — 77 ^- ) a m<X P ( v mn + /3 n W mn ) 2 W m „fl5ll7 

9 R R 


m=l n=l 
M N 


+ YYh^. + H^)(t, mn + f3 n w mn )a 2 m w 2 rn J q hhu 
2 R R 


m = l n=l 

and the elements of [F e s ^ ell ] are given by 

M , N 

' 2 0t<; ™ + X 2 




m vv 1 

F^it (1) =(/? + li>o) 2 © + X O® 10 ” 1 + 5^ 2 0 ^" + ^ + X OPnVnWn 

m = 1 »=1 n«l 

M N , M N i 

+ E E + «’-») + E E 2 eft> ” 

m = ln=l m-ln=l 

M 

F she ii (2) =4 (fl + u> o )0/ + 2(i? + w o )09o + X 2 a m 0 /M m w m + X X Qm0/ 

m = 1 m=l n=l 

M N 

F shell (3) =2(tf + MJo)«m0^m^mp + XI "p 0 ™™^ 1 + X a m&l{v n + /5n»n)«mn^ 


M N 


V'mn'Wmn 


m— 1 n — 1 

M N 


+ X (x m Ql((3nVn + w n )w mn 6 mp + X X 2 ° p0 ^ mTl Wmn ^ 1 


n= 1 

M N 


rn= 1 n=l 

M iV 


aj i_ o \ 2 

+ X X o Qm ^ n0umnU ’ mn ^ 4 + X X Y p ^ nQvmnWmnl1 


m = l n= 1 


m = l n=l 


M 


F s e ^j, (4) =2( i? + U?o)<*m0^m^mp + ( 2? + 9o )&w m 6 mp + X 2a m 0u m w m /i 

m=l 

N M N 

+ X (XmQliflnV n + W n )u mn 6 mp + X X a rn®U mn W m nh 

n=l m=l 71=1 

M N , 

+ EE -0/? n [a m (Ji - h) + 2a p h] ^mn^rnn 


m = 1 n=l 
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M 


M TV 


Khali ( 5 ) =( 2/ + 9o)(«n + 0 n W n )Qb nq + OtmQlfinWm ^mn^nq + EE«-'* 


m = 1 


m = 1 n=l 


M 


M TV 


2 


771 n + PnW mn) Um^nq + EE 5 "I" fiqls)^ 

m= 1 m=l n=l 

M 


77171 ^ 771 71 




^ng 


^ she// (®) — (2f + < 7o)(/^n^ , n “I" ^n)®^ng "1“ ^ ^ 

m=l 

M M TV 

+ E ^m®^(/^Ji^mn "1“ %n) + EE Q J m /iZmn^mn^6 


771 = 1 

M AT 


771 = 1 71=1 


+ EEr m^[/^n(-^6 ~ A ) + 2/?g/5] ^77171^77171 

771=1 71=1 

^ shell ( J ) = ( ^ “1“ ^0 )®m QlWm n^mp^nq 4" {fin^n “I" ^ti ) ct m ©^^Tn^mp^ng 
A/ M 

4" ^ ^ ^p©^m^mn-^l^nij “1“ ^ ^ ^ /^ti © ( £* 771 ^4 4“ ^p-^l)^ 771^77171^ 

m = l m— 1 

TV TV ^ 

4* ^ ^ ^7n^^n^mn/6^mp 4~ ^ ^ m ^ ( /^ti -^6 — 

71 = 1 71 = 1 

TV N 

+ E Qm^n^mn^mp 4* ^ ^ ^ ^ m/^n ^(^71 ^mn -^3 4" 


ng 


J mp 


mp 


71= 1 

TV 


71=1 


M TV 


v— > 1 ^ ^ ^ 

+ 22 ^ a ml(finh + (iqh)v n W m n^mp + 22 ^2 ^ a p( v mnh + W;„/s)A 


g -»ft- \/- »t-u > i-y-v / - t* — n*t* - tny - / v / y 2 

71=1 771=1 71=1 

M TV 


+ E E W ^771 A 4" &pll )(l3 “h -^6)^77171^77171 

771=1 71=1 

Fsheil ( 8 ) = ( v n + 0 n W n )oc m Qlu m 6 m p8 n q 4" (^77171 4~ fin w mn ) ©^tti p 8 n g 

TV 

T ~ Af mp 


4“ cy © ( ^77i 7i 4“ fin Wmn )<70^mp &n g 4" ^ ^ 

71=1 

A/ M 

4” ^ ^ 77i © ^ m^mn h^nq + 22 3^ n0 ( ^ 77i 1 2 4“ ^ng 

771= 1 771 = 1 

" 2 1 
4“ ^ ^ ^Z^-mfin^^m'^mn^2^nq 4“ ^ ^ ^ /^n © ( ^ m -^1 ^p^)^m^mn^7 

m = l * 77i=l * 

A/ 2 ^ 1 

4“ y — Q m fin © ^ in ^ 77i 7i 1 1 ^ nq 4“ ^ ^ r^m^n^5 4" /^g-^6)^n^ 77171 ^mp 
1 ., *5 

771=1 71=1 

M TV M TV 

+ EE ^771^77171^77171-^1-^5 4“ EE -OL m (finh + fiqh){h 4- h) u 


77171 ^771 n 


m = 1 n = l 


771=1 71 = 1 
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F ( J L u (9) =(/? n Vn + W n )a m Qlu m 6 mp 6nq + {PnVmn + W mn )Ql6 mp 6 nq 


+ 


i0(/? n V m „ + w mn )qo^mpKq + {R + wo)a m 0lu mn 6 mp {>nq 


Aj ^ 

+ ^2 amQUmWmnhbnq + ^ + a p U)u m V mn 6 nq 

m = 1 m=1 

M 0 N 1 

+ ^ ~rOl rn ji n Qu m V mn I\(>nq + ^ -a m l{Pnh + Pqh)VnUmnt>mp 

m = l n= ^ 

M JV 2 

+ V amQw m Umnh6nq + Y X a m/3n^n«mn ^ mp + Y a m lw n U mn I 6 d mp 

tlx "=1 3 "=» 

M A' ™ ™ 

+ V' -Pn(0!mh + Ot p U)UmnVmnh + 2^ / , a rnU r nn w mnh h 

tlx tlx 3 »=1 

M AT A A 1 

+ EEW«- / >+“ p l4)u mn Vmnh + /_/ ^ ^ Q 'm/?g(-^l + h)u m n v mn 5 


i = l n=l 


m = l n—1 


Elements of (2 M + 1) X 1 residual force subvector [RstA^str )] 


M 


Rstr ( 1 ) = — ^ Q'% 5 (EA) s a m W 8 m 

~ m= 1 

M 

R sir (2) = (EA) s a^ n lu srn 6 rnp + ^ QMEA) s a m cx p w srn l 2 

771 = 1 

M 

(3) = [Q.5(EA) s qi + (E I) s ot 2 m l} a 2 m w 8m 6m P + ^ 0.5(EA) s Q^Q>u; 5m /i3 


771= 1 


Af 

+ Y^(EA) s a 2 m a p u sm U)smh 

771 = 1 
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Elements of ( 2 N + 1 ) x 1 residual force subvector [i? r (u r )] 


Rr (1) =2 { -^ew r0 + £ i ~lZ§ L Q(v r n + Pn™rn) 2 


\{EA)r , (El) 


(EA) r , 2 (£J)r 


tfr ( 2 ) = p— + i-= 5 ^])S* 0 t> r „« n , + [ L — ^ 

/to tto ^0 

N 

4 “ ~ ^v> ®(^Vrc 4 “ fin'Wrn )^r0^nq 4 " ^ ~Jyi figi^rn 4 " Pn^rn )~-^3 


+ E L f^'»'» + ^“™) 3 - , ’ + E i ^r ! 

n=l ® n= 1 ® 


(^rn “I" ftn ^rn )(/^n ^rn 4" 1V rn )I\ 5 


d /0 \ _\{EA)r , ■} j: , \{EA) r { ^(EI ) r - x 

i? r (3) — [ n 4" ft n j Pn^^rn^nq 4” [ D "I" /9 n ^3 j 0UVn*r 

/to /Iq /tQ /to 


(£A) t 


/^n©(^rn 4 ” /^n^rn )^r 0 ^riqr 4 " ^ 4 ” /^n^rn) ^7 


“I" ^ " 2R~ ^ rn 4" ^3 4" ^ ^ ^2 Pgi^rn 4“ ^n^rn)(/^n^rn 4" ^rn)^ 5 


(£A) r 


The integrals Ji to I\$ in the elements of the residual force subvectors are given by 


I\ = J C os 2 (a m x)C os(a p x)dx 

-1 

1 

I2 = J Sin 2 (a m x)Cos(a p x)dx 

-1 

0 

h =J Sin 2 (P n 0)Cos(p q 0)dO 

-© 

1 

I 4 = J Cos(a m x)Sin(a m x)Sin(a p x)dx 

-/ 

0 

h=J Cos(M)Sin(M)Sin(f3 q 0)d0 
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h 

h 

Is 

h 

ho 

hi 

In 

1 13 

hi 



Cos 2 (i3 n 0)Cos(l3 q 6)d6 
Sin 3 (fi n 6)Sin(fi q d)dO 
Cos 2 (fi n 6)Sin(fi n 0)Sin(fi q 8)d0 
Cos 3 {a m x)dx 
Sin 2 (a m x )Cos(a m x )dx 
Sin ( CV771 X }C Os( QrnX 'jC OS(dpX 'jdx 
Cos 3 (a m x)Cos(a p x)dx 


I 

■ j Sin 3 (a m x)Sin(a p x)dx 

-l 


+e 

: J Sin 2 (fi n 0)Cos(fi„0)dO 
-e 


1 

j Cos 2 (a m x)Sin(a m x)Sin(a p x )dx 


© 

/ 16 = j Cos*(PnO)M 
-© 


© 



Sin 2 (P n 9)Cos{j3 n 0)C o$(3q0)d0 


hs 


© 

: J Cos*(M)Cos((3 q 0)de 

-© 
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